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Compressibility Effects in Circumferential 
Inlet Distortion in Axial Compressors, III 


(Application of the Actuator Disc Theory)* 
By M. Z. v. Krzywoblocki, Michigan State University, U.S. A. 


Summary. This is the third paper of a series dealing with the compressibility effects in circum- 
ferential inlet distortion in axial compressors. The first paper contained the approach related to 
the method of Ehrich (valid in the incompressible fluid flow) and the second to the proposition 
of Y eh (valid also in the incompressible flow). The present approach refers to the philosophy based 
on the actuator dise theory, developed by Bragg and Hawthorne in the case of a compress- 
ible fluid flow. 


Introduction 


The present paper is the third one in the series of papers dealing with the compress- 
ibility effects in circumferential inlet distortion in axial compressors. In the first paper** 
of the series the author generalized to compressible flow regime the method of Ehrich® 
of analysing the circumferential inlet distortion phenomena in axial compressors in 
the incompressible fluid flow domain. In it the author derived certain similarity rules 
which enable one to transfer the results from the incompressible flow domain to com- 
pressible flow regime. Besides that, a straight-forward generalization to the compress- 
ible flow was worked out. The same technique was applied in the second paper*** in 
the series to the method of solving the circumferential inlet distortion as proposed by 
Yeh’. In the present, the technique of the generalization is applied to the philosophy 
of the actuator disc, proposed by Bragg and Hawthorne}. 


1. Fundamental Equations 
1.1. Basic Notions 


Assume a curvilinear orthogonal coordinate system {s, n, 6} with ds denoting an 
element of the arc along the space streamline, o or s; dm an element of the arc along 
the normal to the streamline (curve tangent to normal to streamline at a point P), 
and db along the binormal. The fundamental Cartesian coordinate system {2, y, 2} 
will be denoted by the symbols {y', y?, y?} = {y'} (i = 1, 2, 3) and the system {s, n, b} 
by the symbols {z'} (i = 1, 2, 3). Let a set of equations of transformation: 


RY, on Ee) (1205) 

* The results, presented in the present work were obtained during the author’s association 
with the Compressor Aerodynamic Design and Development Unit, General Electric Company, 
Cincinnati, Ohio, U. S. A., with Mr. R. A. Novak as its Supervisor. The author wishes to express 
his deep appreciation to Mr. R. A. Novak for the numerous discussions on the subject in question. 
For the heavy task of typing the present work, thanks are due to Mrs. Maxine Spillman, the 
secretary of the unit. 

** (sterr. Ing.-Archiv XIII, 214 (1959). 

*** Ogterr. Ing.-Archiv XIV, 79 (1960). 
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define the curvilinear orthogonal system {2}. Since the Cartesian system is determined 


by a set of orthonormal base vectors ‘ba (i = 1, 2, 3), the position vector r® of any 
point P (y}, y?, y3) can be represented in the form: 


sf eae AU Gail ee MeL ONS 2 
p= DY. Oa == Sar aaa (1.1.2) 
(ds)? = gi? dai dx. (1.1.3) 


Assume an (orthogonal) cylindrical polar coordinate system {6, 7,2}. Instead of the 
coordinate “‘7’’, we shall introduce the reduced radius ‘“‘R”’ referring to the unblocked 


circumference: 
R=r—(2a)s?*tM=rK: K=1— Qar)*tMm, (1.1.4) 


were the symbols denote: 


¢ — thickness of the blade at certain r and z (coordinate measured along the axis of 
the compressor) ; 
M — number of blades. 
The quantity K is assumed to be a parameter, not affected by the differentiation 
or integration operations. We shall use the system of coordinates: 
Ss Kise Nin 6 a6 Oe aera (1.1.5) 
(ds)? = gi dX dXs; gf = K+ gi. (1.1.6) 
We shall also use cylindrical polar coordinates {6,7,z} or {O, R, Z} where R = Kr; 
Oi On Ate= K 2. Krom: 


T: R=Kr= K (a yy)"; 6 = tan 2); Ze (1.1.7) 
TT: ¢ = 700s 0—= BK" e630. ih i Sic == a (1.1.8) 
(ds)* = R? K-? (d0)? + K-2 (dk)? + K- (dZ)?, (1.1.9) 
we have: 
Pe Rea! iT elery 
g = [gj] =| 0 BOO aa Ree (1.1.10) 
0 OF ee) 
a: {(K-SRK, 020) sisi) (OKO) 5 Gas (OL 0, ke): (lta 
LRA 0 
(ga = | 0 TOA mE ea See (1.1.12) 
pao Oy rr? 
In this system, the velocity and force components are: 
Velocity vector: C = a; A'= ai As; (1.1.13) 
C: Physical components: C,; C,; C:; (1.1.14) 
Contravariant components: 4° = K R- 0, ; 
A®= K C,; (1.1.15) 
AX s2K-GC,* 
Covariant components: Ager RC, 
Ar = K-1 Gs (1.1.16) 


Az= Ce: 
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Force components: 


Physical: R, ; R,; .R-: ; (1.1.17) 
Contravariant: P’ = K R-! R,; P7 = K R,; P?=K R,; (1.1.18) 
Covariant: Pe = K-1 R R,; Pe = K- R,; Pz = K” R: (1.1.19) 


Force potential @: 
Contravariant components of the gradient: 


Vo= ih 0D /0x; = igh OD oa" ; (1.1.20) 
a = KR: @° = K? R 06/00 = K* RG, o- 
pie BEL ie ow eae de. (1.1.21) 
a, = KA: R= KG, » 


Concerning the parameter K, the following philosophy is accepted: Let us follow the 
principle that one should deal only with a three-dimensional coordinate system, like 
Cartesian {2x, y, z\, cylindrical polar {rs 6, z}, etc. In some cases it is more convenient 
to use another coordinate system like R = Kr, 6 = 0, Z = Kz, where the parameter K 
is actually a scalar point function of {r,z}. But when calculating the square of the 
element of arc, (ds)?, the assumption that the function K should be subject to ordinary 
functional variation at a point (dK = K,,dr + K,,dz), leads to unsolvable complica- 
tions: Either one is compelled to preserve the expressions dK and (dK)? in addition to 
(dR)?, (d@)?, (dZ)? which would mean that the transformation (7, 0, z) > (R, O, Z, K) 
involves a change in dimensions of spaces in question (from 3 to 4) and consequently 
is not a one-to-one transformation; or one is compelled to use the total differential 
dK = K,,dr+ K,,dz, to attempt to express it in terms of dR and dZ (if possible) 
and to insert it into the general equation for (ds)?. The partial derivatives K,,; K,;; 
must be assumed to be known. Actually, they depend upon the mechanical design in 
question. But this will complicate enormously the entire mathematical procedure beyond 
any practical use. The problem is, that in Kq. (1.1.4) the functional relation ¢ = ¢ (r, z) 
depends upon the mechanical design and is unknown in its most general form. Hence, 
the exact explicit forms of the relations: 


Poet ay Koa Kah, Z); (1.1.22) 


are unknown. In practical applications, it is sufficient to assume certain mean values 
for the function K. To avoid all the complications, mentioned above, the assumption 
is made that the function K behaves like a point-wise constant parameter, i. e., instead 
of a point function K, one assumes a point step-wise function K. Thus, when the opera- 
tions of differentiation and integration are applied, the parametric function, K, is not 
affected by them. ; 


1.2. Fundamental Equations in Steady State Conditions 


Equation of motion (Euler’s equation): 


(C-V)C+o01Vp=P=V72. (1.2.1) 
We assume that the external force has a potential ; 


continuity and state: 


J 


div (0 C) = 0; p= ko T; R, = gas constant; (1.2.2) 
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pressure-density relation: 
P Pot = (0 oo)” exp [1 (S — Sp)];_ S = entropy, (1.2.3) 


where the subscript zero denotes some uniquely defined initial conditions (either some 
stagnation conditions or the conditions in the incoming uniform flow) ; energy equation : 
in steady state conditions a particle moves along a streamline; assume some initial 
conditions, denoted by the subscript zero; other points on a streamline are denoted 
by the symbols S,) S, etc. Assume that on the length of an arc [0, S,] the velocity is 


denoted by Cy, pressure by py, density by Qo, velocity of sound by co, and there are no 
external forces and no heat addition from outside; on the arc [S,, S] the used symbols 


= > 


are: C, p, 0, c, there are acting some external forces of the vector sum P and some heat 
from outside is added to the particle of the total amount measured from the point S, to 
the point S equal to Q (S). Then the equation of energy along a streamline has the form: 


£02 ery (y — 1)? py Op7) = 8 CO? ey Ona CS) 
S Ss 
—|P-do=G(y,9); | P-do = (8), (1.2.4) 


& S, 


where a streamline is characterized by y = const., and ? = const. 


1.3. Cylindrical Polar Coordinates Forms 


The equations, given above, are expressed below in cylindrical polar coordinates. 


(a) Continuity (with the use of Eq. [A-3]): 


(oCy),e + (RoC) 2 + (hoC,), z = 90; (1.3.1) 
(b) Motion (using equations given in the Appendix): 
C,Co,, + 7 OC, Ca rr OC, Cy CLC, =o 47 pin 1 as (1.3.2) 
Cr Cnr + 1 Oy Ong = 7 tC + C,C2 = —0" p,, -- Be; (1.3.3) 
C, Cr + 7? Cy Cag + C02 = — op, . + Bz. (1.3.4) 


In (O, R, Z) coordinate system: 
C, Car + B+ CO, Oo,5 + BAC, Cy + 0,0, z= — 01 R41 p,, + R3@,,; (0.35) 
C, Cir + R+"C, C9 — R4C,? + C,C,, z= — 0 1p,r+D,p; (1.3.6) 
C, Cr+ RC, Cu9 + C.O,2= — 01 7p,z4+ G,z. (13.7) 
Equations of state and of the pressure-density-entropy relation remain the same. 


(c) Vorticity components (referred to the a;-system): 


K(Cor— C2) = (1.3.8) 
KR (0 = 0, =a BC, = OM; (1.3.9) 
BRO C Oo ee (1.3.10) 


(d) Energy equation: 


3 (Coo + Ceo 4 Cro) rey 1)-* pp Op? = 4 (C2 + CP + C2)+y (Pee Coe 
— Q(8) — ® (8S) = y (vy — 1) 1 py eg ' = 4 CPF — Q (Si) — G (S) = G4 (y, 9), (1.3.11) 
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where the subscript “‘st’’ refers to the properly chosen stagnation conditions where 
Q (Sx) = 0 and © (S,) = 0, and the subscript “7’’ refers to the limiting conditions 
with p (S;) = T (S;) = 0. The velocity C; may-be assumed to be constant on a particular 
streamline but may vary from streamline to streamline. The symbol S, denotes the 


value of the running coordinate at = on and @ (S;) the total amount of heat added 
on the distance [S, = S/], i. e., [S, — S,]. 
1.4. Generalized Crocco Equation 
From the first law of thermodynamics 
dQ =«c,dT + pdv; T- dQ = dS, (1.4.1) 


after some elementary transformations, one obtains: 


odp = dH — TdS — 4d (C2), (1.4.2) 
with: 
h=oT; H=h+}4C?; @=C024 O24 C2, (1.4.3) 
or: 
01V,=VH — TVS —4P (C2). (1.4.4) 


Using Eq. (1.4.4) in the equation of motion (use also Eq. [A-7]) furnishes: 
VH —TVS=VO+ Cx curl C, (1.4.5) 


with the components (referred to the a;-system) : 


Kenai. = 1S, ih ho O®,, 42°C. £=C,¢; (1.4.6) 
Ke Sn De te Cy. = On: (1.4.7) 
K(H,z— 78,2) =K@,z+C,4—C, 6. (1.4.8) 


Another equation to be considered is: 

C-(VH — TVS — 7) = 0. (1.4.9) 
Since the force acting on a blade is not perpendicular to the direction of the incoming 
flow, one gets: 


== SS 


C-(VH — TVS) =C-V@ =|C|-| P|: cosa. (1.4.10) 


1.5. Streamfunctions 


Equation of continuity is satisfied by means of two functions, p and #, fulfilling the 
equations: 


Oo C, iy D(R,Z) : (135,12) 
D (y, 9) > 
Rig C, NET oy (1.5.2) 
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1.6 Transformation of Equations 
Kq. (1.3.6) gives: 
R[C,(RCy), r+ C,(RCy),z + Cy R>(RCy).o] = —O7R* pet RD, 4; (1.6.1) 
or 


= d > 
CV (RO) = (BOs) = — 07 Py + Do (1.6.2) 


Similarly, Eqs. (1.3.6) and (1.3.7) furnish: 


aC, 
gp BOP = 0 a Oe (1.6.3) 
10; , 
GE SO Re Oe (1.6.4) 
S = C,0/0R + RC, 0/00 + C, 6/eZ. (1.6.5) 


Calculate the ratio (p o-1) from the pressure-density-entropy relation and insert it into 
the energy equation, thus obtaining the expression for the density 0: 


e = {[3 (C? — C?) + Q (S) — @ (Si) + @ (S) — @ (Si) - 


“(y= 1) poe po exp (Oe se ee es (1.6.6) 
The velocity of sound is defined in the usual way: 
C? = (@p/e)5~ const. = Y PQ =y RT. (1.6.7) 


Using energy equation (1.3.11) one gets: 
02 = 4 (y —1)(— 0%) + (y — 1) [© (8) + Q(8) + @]=y RT =yper. (1.6.8) 


Let us take the logarithm of Eq. (1.6.8) and divide it by Eq. (1.6.8) itself, thus obtaining 
after some operations: 


— 4C,-* grad C? = (y — 1)-1 grad (In 7') — C,-? grad [@ (S) + Q (S) + G]. (1.6.9) 
Insert Eq. (A-7) into Eq. (1.2.1) and divide the result by C,?, thus obtaining: 


=> 


(o x C) C,-* = — $ C,-* grad C? — y— grad In p + C,-* grad ©, (1.6.10) 
or using Kq. (1.6.9): 
(w x C)C,-* = (y — 1)-! grad In T — y-! grad In p — C,-* grad [Q (S) + @]. (1.6.11) 
But the first two terms on the right hand side of Eq. (1.6.11) give: 
grad In [pi/y (Ra @)—1] 


(vy — 1)-' grad [In T — (y — 1) yn p] = per : (1.6.12) 
or with 
p = Co’ exp (S c,-1); p (CC 9) = exp (Se, 4); (1.6.13) 
(y — 1) grad In T — y~! grad In p = (yR,)- grad 8. (1.6.14) 
But with 7Q = TVS (dQ = T dS): 
(y Ra) VS —c,2 TVS = 0, (1.6.15) 


and Kq. (1.6.11) furnishes: 
(wo x C) = — VG. (1.6.16) 
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Combining Kgs. (1.6.6) and (1.6.8) one gets: 
@ = {v1 C2 9’ py exp [(S) — 8) eA] 0-9, (1.6.17) 


Introduce dimensionless quantities, using as the reference velocity the velocity C), 
assumed to be constant on each particular streamline: 


G C7} = W (Cy; Cry ¢.); C, C;-! =C;. (1.6.18) 

Then: 
c? = $ (y — 1) (— w®) + (y — 1) Cr? [@ (S) + Q (8) + GI]; (1.6.19) 
@ dimensionsless = 97 = {y-1c.?_ 0," py“ - exp [8S — 8) e1]}/¢—-, (1.6.20) 


Let us multiply scalarly Eq. (1.6.16) by C: 


> > > >_> > > 
(ax C)-C = [oC Cl]=—Veé-C=—0, (1.6.21) 
since in the present case, the triple scalar product is equal to value of the determinant 
in which two rows are equal (such a determinant is equal to zero) and, on the other 


. . a . . . 
side, along a streamline /G = 0, since G is constant along a streamline. 


2. Running Coordinate System 
2.1. Fundamental Notions 


Let the Cartesian coordinate system be denoted by y' (¢ = 1, 2, 3). Assume at each 
point of a curve (streamline) an orthogonal running coordinate system which princi- 
pally is a cylindrical polar coordinate system consisting of two coordinates located 
in the osculating plane and of the third one perpendicular to it. In the osculating plane 
there are located: 

(i) the radius of curvature, r, of the curve, K, at the point P; ds is the length of the 
are of K between two points, P and P;, on the curve K; dO is the angle between 
the tangents at P and P,; the principal normal at P is denoted by the letter p; 

(ii) the binormal, b, at P is perpendicular to both the tangent and the principal normal 
at P. 

One has: 


==i¢ == 7-!; » = Curvature of K at P. (Zale h) 


We choose the system of coordinates x’ = (6, p,b). This is an orthogonal system, 


hence, the quadratic form is: 
- Oy dy 
(ds)? = Gjk agidz; Gb Ty ak 2 (2.1.2) 


with gj, = 0(j + k). In the present case one has: 


(ds)? = (x-1 d6)? + (dp)? + (db)?. (2.1.3) 
This furnishes: 

ee ) 0 
g = [gi] =| 9 1 O|=x*; (2.1.4) 

0 0 1 

Pe 0 0 

g = [g]=| 0 1 0; =x? (2.1.5) 

0 1 
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Since p and r are located on one and the same line, with p pointing towards the center 


of curvature and r from the center of curvature, one has: 


ap = or; = . (P3) =r t= we. (2.1.6) 
ag 4a ji0y 0): bar iO, Oey ged) (2.1.7) 
Velocity vector components: 
Physical: C = (C7, 0, 0); Cy = | C : (2.1.8) 
Contravariant, Atl Ces 03 0; (2.1.9) 
Covariant, AR GE O55 (2.1.10) 
Curl C: fae te DOG 
ox x e-LiCe spas 
Vee — “(47 Ce),p = 4 OC, — Cry». (2201) 
The continuity equation reduces to: 
(0 C5),5 = 9; @e Cy = (p, 6) only. (2.1.12) 


The function g can be only a function of p and b (not of 6). The stream functions, y and #, 
are constant along a streamline. From Kq. (1.5.1) one has: 


(2.1.13) 


eG = 


which may be interpreted in the sense, that the Jacobian of the transformation 
{p, b} <> {y, 3} is equal to (9 C,;). Let us introduce new coordinates: 


O@=60;P=Kp; B=Kb; R= Kr, (2.1.14) 
then one has: 
(ds)? = K~-? R? (d0)? + K~ (dP)? + K- (dB)?; (2.1.15) 
| KR? 0 Ota rk 
g=([gil=| 0 kK? 0 |= K-®R; (2.1.16) 
0 0 Hite 

a: KR; b: K- a,: K-}; (2.1.17) 
C: Physical: Cee 080; (2.1.18) 
Contravariant: KRG; 020; (2.1.19) 
Covariant: K-1 RO;; 0; 0; (2.1.20) 
Force P: Physical: Ry Rips Rss (2.1.28) 
Contravariant <A’: K R- R;; K Rp; K Rp; (2.1.22) 
Covariant A;: K- R-R;; K-1 Rp; K-* Rp: (2.1.23) 
Curl C: ay: 0; a3: K* Og, 9; ag: K7R2O; —e CG, 2: (2.1.24) 

Continuity equation: 

D (y, 9) 


e Cy = K* Dy; Dy = Dy (y, 8) = arp gy Ks (2.1.25) 
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Gradient: 
ay: K? R- 2 618 6; dg: K?0/0 P: a9: °K? alo B: (2.1.26) 
Cx curl C= p(n’): 
Ga py =O: (2.1.27) 
ie peek, Onn — KOR O22: (2.1.28) 
Os: p> — K2 C. Ce ip. (2.1.29) 
Also from Kq. (2. 1. 13): 
q! D (y, 9) 


2.2 Equations of Motion and Energy 


Components of Crocco’s equation: 


Hy; —T 8,5 = ®,;; (2.2.1) 
H,p — T S,p = @D, p + Cy Cr, P — R> Cres (2.2.2) 
Hys— TS, = ©, 5 + Os Cy, 2. (2.2.3) 


The acceleration vector is (see Eq. [A-9]): 


A=(C-V7)C =a, K* RC; O;,; + a, K? R- C52, (2.2.4) 


with a,, a, given by Eq. (2.1.17). With Rd = dt,t = tapegnt to the curve, Eq. (2.2.4) 
is in a full agreement with Eq. (15) in 4. The vector A lies in the oscullating plane. 


Euler’s Equations: a, (tangent or streamline direction): 


KR* 05045 + KR 0 pj = Ri; (2.2.5) 
a, (principal normal direction): 
K Ro C52 + Ko-! p,p = Rp. (22.6) 
Energy equation: 
PC ty 1) pes — (8) ~O() = Gy, %); (2.2.7) 
entropy: ve 
S (y, 9) = Ba (y — 1)* In (pe); (2.2.8) 
or: - . = 
SY — 8 (y R.)-; S® (y, 8) = y-* (vy — 1) * In (po). (2.2.9) 


Let us construct the gradient component of Kgs. (2.2.7) and (2.2.9) in a@,-direction 
(or P-direction): 


OF Cy,ep + v (vy —1)7 po! (pp, ep — @ 1 0, P) — (Q + ©), Pp = Gyr; 
(y — 1) (y p> pp — 07 @,P) = 8%; (2.211) 
this enables us to construct the difference (with the return back to S): 


G,p — T S,p = Co Cy, ep + 0 pp — (Q4+ ®),P, (2.2. 2) 


(2.2.10) 
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where Q, p = 7’ S,p; ©,p = K- Rp. Insert the value for e-! p, p from Kq. (2.2.12) into 
Eq. (2.2.6), thus obtaining [see Eq. (2.1.24)]: 


G,p = — K-* Cy £5; és = K* (RC, — Cy. 2). (2.2.13) 


Apply the same technique to the gradient of the energy equation in a, (— B)-direction, 
thus, obtaining (in this direction the acceleration vector does not possess any component) : 


Ko" p,a = Ra; (2.2.14) 
GB = A Cy Ep; fp = 1 .& Ch, B, (2.2.15) 
with &p [see Eq. (2.1.24)] denoting the component of the curl C-vector in a, (P)- 
direction. Let it write the equation of motion in the form: 
A+eo3Vp=VO; A=34V7(C-C)—CxeurlZ, (2.2.16) 
and let it apply the gradient operator to Eqs. (2.2.7) and (2.2.9), thus obtaining again 
Crocco’s equation in the form: 


VG= Cx eurl C, (2.2.17) 
Using Eq. (A-1), one obtains from Eq. (2.2.17): 
a, K*G, p + a, K*G, p= — A, Cz én + as Ca Ep. (2.2.18) 
2.3. Transformation of Coordinates 


There are no components of curl C or of (w x C) in C-direction [see Eqs. (2.1.24) and 


(2.1.27)] (briefly, ae C, wx CC); this implies that one can consider a one-to-one 
transformation {P, B} <~ {y, 0}. The technique is a well-known one; from: 


dG = G,,dyp + G,» dd, (2.3.1) 
one gets: 
G,p=G,yyie+ GoD, p; (2.3.2) 
Gjp=G,yya + Go, .B; (2.3.3) 
solving for the unknowns G, , and G, 5, furnishes: 
D Y 
G, y = (4,72 8,2 — G2 8,7) J, = DEE Jats (2.3.4) 
D (G, ) D (y, 9) 
Ge=— ppp: J. = D(P.B)- (2.3.5) 
Insert into Eqs. (2.3.4) and (2.3.5) Eqs. (2.1.30), (2.2.13), and (2.2.15), thus obtaining: 
G,y = — 01 (EP B,p + &3 9B, a); (2.3.6) 
Gio = 01 (Ery,p + sy, zB). (22327) 
Similarly, from G = G (P, B) and: 
dG = G,pdP + G, 2, dB, (2.3.8) 
Gs i Ge Pie + GB By yi Gy See 
one obtains: 
D(G,B DG 
fe Dante I Ga = — STE n JA, (2.3.10) 
D(P, B) 
T= Dye Ii = KE (9 OG). (2.3.11) 
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Insert the values for G,p andG, g, Eqs. (2.3.10), into Eqs. (2.2.13), (2.215), and calculate 
two derivatives: 


Gag Gr Ba Sa P,y) (0d) Gee = (Er Bin = En Pye) (0). (2.3.12) 


Using Eq. (2.1.30) in conjuction with the expression for J,, Eq. (2.3.5), one can cal- 


culate the derivatives C3, p and Cj, ,; these inserted into Eqs. (2.2.13) and (2.2.15) 
furnish : 


Cr ho (Rea ots ee Pe 15: (2.3 fa)e 
G,p= — K*o? (0,087) — J, Jy, 2). (2.3.14) 


This approach can be generalized. Consider Eq. (2.2.17) in the {0, P, B}-coordinate 
system. Using the expression (2.1.30) in conjunction with the second Eq. (2.3.5) for 
Cy in Kas. (2.1.18), (2.1.19) and (2.1.20) one constructs the derivatives Cy. p and C7, z. 
These are inserted into Eq. (2.2.17) [see Eqs. (2.1.27) to (2.1.29)], furnishing after 
various operations: 


VG = K*91J.2 R37 (Ro) +4 K4 02 7 (J), (2.3.14) 
or 


VG = K*(0 RB) J, V (Re J;), (2.3.15) 


where the gradient V has only two components in P (= a,) and B (= as) directions. 
Another useful form of Kq. (2.3.14) is: 


VG = — K*[o3J,2 Vo — 0? R2J3,V (J, R)). (2.3.16) 


One can construct various other formulas of value in some future applications; using 
Kgs. (2.2.7) and (2.2.8) or (2.2.9) one constructs the expression: 


VG—yTVS=V (Cs) + ypo*Ve -—V(Q+ ©); (2.3.17) 
write Eq. (2.3.15) in the form: 
VG = —K*to3J2Vo0 + KJ, 02 R37 (RJ), (2.3.18) 


calculate Vo from Eq. (2.3.18) and insert it into Eq. (2.3.17), thus obtaining after a 
series of operations the formula: 


(20,2 — 1) K*R*S,V (RI,) — Cs KI PV Cpe Of) = 
= ot (— VE + (yy — 1) yy B.1 0)27 VS — Ga? C2 FO (S)). . (2.8.19) 
Using Eggs. (2.1.13) and (2.3.5), the second term in Eq. (2.3.19) is equal to: 
— OF Ot J Ki yS 5, (2.3.20) 
which furnishes another form of Eq. (2.3.19); 
K* R23. J, (02,2 J, VR —V(RJ,)] = &([—VG+ 
+ (y= 1) pt Ret G9 BS 305? C,-* VS (8)), (2.3.21) 


where the gradient / has the components in two directions (P and B) only. These com- 
ponents are: 


EGR 3 (O10) Sd Ja = LS Gp + 
+ (y Ba)" (y — 1) O72 S,p — CG? C,-? ®, P); (2.3.22) 
an — K4J, Js = oa {[— Gye + (y dig) (y —1)- Cy? a Cra ®D, 2). (2.3.23) 
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Eqs. (2.3.22) and (2.3.23) can be transformed. To this end let us begin with the func- 
tions P = Pty, 3) andi =" B97 0): 


dP = P,ydy + P,»d0; dB = B,y dy + B,o ad; (2.3.24) 
L=P,yyp+ P,od,p; O= Py ya + Pye FB; (2.3.25) 
i Veg Bra yy Sl ba wp by yk oe ee (2.3.26) 


these equations jointly with Eq. (2.3.11) furnish: 
P,y=8,03; P.g==v.eJ; By= 26255 Bao=prd. (2327) 


Returning to Eqs. (2.3.22) and (2.3.23), let us express on the right hand sides of these 
equations the derivatives 0/6 P and é/é B by means of Eqs. (2.3.10); next, let us group 
together the terms containing the derivatives B, 3; B,y; P,9; P,; it appears that the 
coefficients standing in front of P,y and B,y and P,,,and B,,, are identical, respectively ; 
let us calculate these two coefficients from the two equations in question; and in the 
final step let us express the derivatives P,,; P, 9; B,y; B,»; in terms of Eqs. (2.3.27). 
The resultant equations are: 


07 [Gy — (y Ba) (y —1I) Gi? 8,y + OF C4 G4) = 
= [(Cs2.C,-2 — 1) K* R= J, 4 K4*d yp] Op — A JS Oe; (2.3.28) 


0 1G,9 — (y Ba) (y — 1) Cy? 8,9 + Cy? C,* ©, 9] = 
= — [(Cy2C,-? — 1) K* R2OJ, + K*J,,ply,e + K*J,,By,p. (2.3.29) 
These equations have the property that on their left-hand sides they contain partial 


derivatives with respect to {y, #}-coordinates, and on the right-hand sides the partial 
derivatives with respect to {P, B}-coordinates, respectively. 


3. Cylindrical Polar Coordinates 
3.1. Fundamental Equations 


In the previous Chapter we have derived the equations referring to the {y, 0}- and 
{6, P, B}-coordinate system, respectively. 

In the present Chapter we shall use the {0, R, Z}-coordinate system. From 
Eqs. (1.5.1) to (1.5.3) one obtains 


CPO alGzsne (3.1.1) 
——. 0 Fie (3.1.2) 
O, = KAR osez (3.1.3) 
CO. = K* (R 9) ae (3.1.4) 


Let us express Kq. (2.2.17) in cylindrical polar coordinates using Eqs. (1.1.20), (1.1.21) 
and the formulas given in the Appendix; hence: 

VG =a, K? R*G,, + a. K?G,r + a, K?G, 2; (3.1.5) 
with 


as 


G,— Kh? hy ag = Ky a Sk. (3.1.6) 


Compressibility Effects in Circumferential Inlet Distortion in Axial Compressors, III 109 


€ Soul | a, K R2(C.t = On) + 
eC. e Ch) 
4s GORUG eer Cre), (3.1.7) 
eee Cre re CGT) — C913 6 7 [C..9 — (Cyr) ret - (3.1.8) 


Let us multiply Eq. (2.3.6) by dy and Eq. (2.3.7) by dd; substitute for the differentials 
0,p; 0,8; v,P; y,e their values calculated from Eqs. (2.3.27); this furnishes: 


dG = J, 01 [EpdB (yp, 0) — Eg dP (y, 0)]. (3.1.9) 


To verify the result (3.1.9), let us reduce the expression (3.1.9) to a two-dimensional 
flow; in this case one has: 


ep = OF ena > 0 == B: (3.1.10) 
Dy, B 
J. = Diep = we = ala, (3.1.11) 
and 
dG = — wo dy, (3.1.12) 


which is in full agreement with the result quoted in °, pp. 87. 
Consider Eq. (2.3.15) in conjunction with Eq. (2.3.11), thus, obtaining (with K 
not being subject to the differentiation operations): 


dG = C; R Rd (RC;) =} Rd [(R G;)?), (3.1.13) 
or 
dG = 4d (C2) + 4 R C52 d (PR). (3.1.14) 
With: 
Oj? = 0? = 0 + OF + C2, (3.1.15) 


one obtains 
dG =} R-*d [(RC,)"] + 3 (RC,)2d (BR) +} R-*0;2d (RY) + 4d (C2 + C2). (3.1.16) 


We now use expressions: 


F=F (y,9); dF =F,,dy + F,» dW; (3.1.17) 
y=y (6, R,Z); = 8 (6, R, Z); (3.1.18) 
Hh = Fry Ys + Fy 9,03 F yr = Fy WR + F 5 OR; 
ete a ze (3.1.19) 
From Eqs. (3.1.19) one can calculate the derivatives: 
D (yw, 3)|-1 : 
Foy = (FP z9.1 — F609, 2) ae a ; (3.1.20) 
D (wp, @))-1 
Fis = (Frey, z— Ff, 24,8) Fara A (3.1.21) 
with the Jacobian D (y,#)/D (Z, 6) given by Kq. (3.1.3); similarly 
D (y, 8)|-1 by 
Fy = (F,¢9,x— F, Rr, 4) Faces > (3.1.22) 
D (y, @))-1 
Fy, 5 =(F,ry,5 —F 64, R) eal 5) (3.1.23) 


110 M. Z. v. Krzywoblocki: 


with the Jacobian D (yp, #)/D (0, R) given by Eq. (3.1.4). Let us perform the following 
re-modeling of Eq. (3.1.16): 

(i) First term: Apply: d = é/dy- dy + o/od - dd; 

(ii) Second term: Left in the same form; 


(iv) Fourth term: C,?: 


) 
) 
(iii) Third term: Left in the same form; 
) 
(a) expressed by Eq. (3.1.3); 

(b) 

(c) 


d = d/éy - dy + a/ed - dd; 
(...),y by Eq. (3.1.20); (C ..)ve. by: Bay (S1aar), 
C2: (a) expressed by Kq. (3.1.4); 
(b) d = d/ép- dy + a/0d - dd; 
(c) (...), y by Eq. (3.1.22); (...),@ by Eg. (3.1.23). 


The results of these operations are: 


UG = R-*(R Cy) (RC), y dy + (RC,), 9 A) + Keo BY {(— Lye + 
a Wz O, 4) dp + (116 Pz Jaye y, 9) dd} ai Kk? (o R)” \(Io,4 9, R a Iz, xB, 9) dy ae 


+ (— Iney r+ Le, ry, 5) 9} + 4 (BO, d (R-*) + 4 R= Cid (R); (3.1.24) 
D ( ma) ry D , 8) 
hs 02K (9 RY) hee BA SR (81.25) 


One may perform some additional re-modeling of Eq. (3.1.24). Namely, one may easily 
derive from Eqs. (3.1.19) a third set of partial derivatives F, ,, and /’, » analogous to 
Eqs. (3.1.20) to (3.1.24): 


D (wy, & 
Fey =f, R82 ize e) Ip eal = ; (3.1.26) 
Foy = (Ff, zy, r —F ry, z) Ip te Al > (3.1.27) 


with the Jacobian D (y, ?)/D (&, Z) given by Kq. (3.1.2). Applying to the term contain- 
ing d (R-?) in Kq. (3.1.24) the operator d = d/éy - dy +2/00-dd, and Eqs. (3.1.26) 
and (3.1.27) one gets: 
} RCP d (R-2) = — R°C,2 R-*(B,y dy + R, dd) = — K*(9 RB) C, (8, 2 dy —y, zd). 
(3.1.28) 
Consider the last term in Kq. (3.1.24): 
R21 0,2dR = R- (0,2 + 0? + C2) dR: (3.1.29) 


the differential dR will be evaluated in three different ways depending upon the velocity 
component it is attached to: 


(i) C,?: we use the operator d = 0/éy - dy + 0/09 - dd, and the partial derivatives 
of the form (3.1.26) and (3.1.27); 


(ii) C2: we use the above operator and Eqs. (3.1.20) and (3.1.21); 
(iii) C2: the above operator and Kgs. (3.1.22) and (3.1.23); the result is: 
} BR O;?d (BR) = 
= K* R> {[0, 04 (B, 28, z— R. 28, e) + 
+ C,(@ R)> (RB, 29,5 — Bio 9, 2) + 
C(t ji (R, oO,r— Rr @, 9) | dy + 


Compressibility Effects in Circumferential Inlet Distortion in Axial Compressors, III 111 


+ [Cho (R, zy,r— R,ry,z) + 
+ C, (eR) (Ry, z2— BR, zy,4) + 


+ C.(@ Rk) (RB, ry,5 — Roy, r)] dd}. (3.1.30) 


Using the concept of 0 = RO,, (1.3.10), and Kgs. (3.1.25), (3.1.28) and (3.1.30), 
one can represent Eq. (3.1.24) in the form: 


dG = RO (60 dy + 03 dd) + 

+ K* (9 B)- {(— Crp 9, 2 + C,, 29,4) dy + (Coy, z— Cr zy): 

dd} + K* (9 RB) {(C., 49, x — Oz. 2 8,9) dy + (— Caspr t+ Cary, ,) dd} — 
— K* 9 R* 6 (8, zdy — y, zd) + 

eK RA 108) 9 (Ry 2 One ks 2 Ok) 

+ 0, 0-1 (R, 29,5 — R,, 9, z) a 

+ C071 (Ri, 9,2 — B,r%,,)] dy + 


+ [0% 97 (BR, zy, 2 — Rey, z) + 
+ C,e> (Ryo y.z— RB, zys0) + 
+ Cro (B,ry.o — Ryoy, r)] dd}. (3.1.31) 
Assume that 
G=G(y,0) and d@=G,,dyp + G, dd; (3.1.32) 
this inserted into Eq. (3.1.31), with the condition that the functions y and # vary 
independently, and that the coefficients of the differentials dy and d# on both sides of 
Kq. (3.1.31) must be equal, furnishes a system of two simultaneous equations: 
Gy = ROM ON) + K? (0 BR) ((—C..5 9,2 + C2954) + 
TAC 00, 2 = Cz, x %)} — K* 9 R72 OM 8, 2+ 
+ K* R* (69 91 (R. 28, z— RB, z9,8) + 
+ C, 07 (R, 29,4 — Ry, 9,2) + 
+ 0,07 (Ryo 9,x— B,r8,s)]; (3.1.33) 
G9 = R76 O% + K? (0 RB) [(C,, oy, 2 — Cn zW0 — 
(Coan — Cz 29,0) + Ke RA OY yz + 
+ K* R* [6 91 (RB, zy,r— Bey, 2) + 


Cro (Riv, 2— BR, ev) + 


+ C,0-1(R,ry,5 — Boy, e)). (3.1.34) 


One can derive another form of equation for dG analogous to Kq. (3.1.13). Namely, 
using Eq. (2.2.7) one has: 


dG = 4d (032) + y (vy — 1). d (p 0) — TS (8S) — d® (8). (3.1.35) 
The use of Eq. (3.1.15) furnishes: 
dG = dH — T dS (8) — d® (8). (3.1.36) 


The differentials in Eq. (3.1.36) may be treated by means of Eqs. (3.1.17) and (3.1.20), 
(3.1.21), or (3.1.22), (3.1.23), or (3.1.26) or (3.1.27) with Jacobians given by Kgs. (3.1.2), 
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or (3.1.3) or (3.1.4), respectively. With Eq. (3.1.36) and with the use of Eqs. (3.1.26), 
(3.1.27) and (3.1.2) one gets: : 


Gy dy + @.9 49 = ((G,29,z—G, 29, 2) dy + 


si (G, 6 Rie G, RY, z) di] Ke (o C4) (3.1.37) 
G,y = (G29, z—G, 8, vn) K? (9 Cy); (3.1.38) 
Gg = (G,.20.8 — Gry ee eee) (3.1.39) 


One may construct a form of a partial differential equation: 


Cy Si G. 9 <= UG, Rr, z ev] GE, R) Os (G, ah ie G, Ry, z)] Ke (o Oy) ® (3.1.40) 


3.2. Another Form 


The forms of Eqs. (3.1.33) and (3.1.34) are inconvenient in that they contain the 
radius of curvature R, which may be an unknown function. For this reason we shall 


derive another system of equations for G,,, and G, 9. 
Starting with Eq. (2.2.17) and using Eqs. (3.1.5) to (3.1.8) one gets: 


R26, = C6, C=O, (3.2.1) 

Goa = CO. Fee (3.2.2) 

G,z= C7 —C, &; (9°278) 

&= (Ci r—C,,2)3 4 = RB [(Cy R),z— Ce]; € = B>[C,,4— (Cy R), er). (3.2.4) 
With C, R = 6, one gets: 


Goo = C, (C4 — OR) — Cz (043 — C255); (3.2.5) 

G, x = Oz (Cyr — O,, z) — Cy R+ (C,, 5 — OR); (3.2.6) 

G, z= CO, R- (6°. — Cz») — C, (Car — C, 2). (3.2.7) 
Let us apply Eq. (3.1.26) and (3.1.27) in conjunction with Eq. (3.1.2): 

Gy = (Girt, 2— Cz RG 0) (3.2.8) 

G,9 = (G, zy, rz — G, ry, z) (Cy 0) K?. (3.2.9) 


Assuming that the Jacobians of the transformations {Z, R}<—>{y, 3} do not vanish, 
we can write with 06 = 6® (y, @): 


ds => 6! “ W,R — 03 0, R; 0, = oe WY, 2 + 0"3 v. ze (3.2.10) 

Let us insert Kqs. (3.2.10) into Eqs. (3.2.6) and (3.2.7) thus obtaining: 
Gir = 0, & — R26 C,,) + BR O® (6% y, p + 0% 9, a); (3.2.11) 
Giz = RAOO OO 7 OY Ge) ROO, Oe (3.2.12) 


Inserting Eqs. (3.2.11) and (3.2.12) into Eqs. (3.2.8) and (3.2.9) furnishes with the 
use of Eq. (3.1.2): 


G,y = BK? (6 9)-1 € (0,0, 2 — C. 4, y— R41 pL 

"(C..6 8, 2 City x) 4a? Gg) ON (3.2.13) 
G9 = — RK? (6 9) (Cy, x + Cry, 2 é +R 04 K?- 

“(C9 Y z— Cony, r) + RO OY. (3.2.14) 
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The system of Eqs. (3.2.13) and (3.2.14) represents the final system of two equations 
in two unknown functions y (0, R, Z) and # (6, R, Z). Having the functions y and # 
one can find the velocity components 103, Of C,}, Eqs. (3.1.2) to (3.1.4), and the vorti- 
city components, Eqs. (3.2.4). The density 9 can be found from the continuity equation 
of the form: 


R-1 (0 C,)) + R-1(RC, 0), r + (0 C,), z= 0. (352.15) 
Another approach to solve the problem in the most general case may be as follows: 


Consider Eq. (3.2.13) as the fundamental one; additionally combine Eqs. (3.1.21) and 
(3.1.23) with Eqs. (3.1.3) and (3.1.4); this furnishes: 


Go (y, zg CP = YR Cs) = (G, R CA ste G, zg OY Ws 6° (3.2.16) 
let us insert Eqs. (3.2.11) and (3.2.12) into Eq. (3.2.16), thus obtaining: 


Gray, zO7* + yr Cz) = R76 (6% y, rk + O89 8, r2) Cr+ Oy, 2+ 
+ 6% 8, z) C7! — (C,,5 C32 + Cae Cp}, (22217) 
or 
G, i) (y, cs Co Bie WY,R Ce == f=? 9© {0% (y, 0) Ca a Oe (y, 0) Ce aay 


—(C,,.07 + C9 C>)}- (3.2.18) 


3.3. Reduction to a Two Dimensional Flow 


Consider an axially symmetric flow with streamlines located in the meridional 
planes; then C, = 0; Eq. (3.1.2) implies that there must be: 3, r = 0, z = 0, and that 
3 = & (8) only, # = 6, say. In Kq. (3.1.30) one has to put down: #,, = 1; all the other 
derivatives with respect to 6 are equal to zero. The right-hand side of Eq. (3.1.30) in- 
volves only the terms (C, R, z — C, R,r). But from the geometry of the dynamical 


system in question one has RC (C,, C,) and: 


C, Oz: = R, r(R, 2), (3.3.1) 
which implies that the term - Rs 

Soha OC, (hy 0; (3.3.2) 
in Eq. (3.1.30) vanishes. Thus, all the three terms containing C, and C5 in Kgs. (3.1.16) 
and (3.1.24) are equal to zero. The remaining terms give: 


dG = K? (9 R)-1{1,, z — Is, r} dy, (3.3.3) 
which with the use of Eqs. (3.1.25) reduces to: 
dG = K? (9 R)- (C,, z — Cz, rv) dy. (3.3.4) 


This is exactly Crocco’s equation °, quoted also in Ref. 1, pp. 244. 


4, Particular Cases 
4.1. Steady Flow With A Rotational Symmetry 


Assume a flow in an axial compressor with an axial symmetry but with C, + 0; 
this implies that all the derivatives with respect to the variable 6 are equal to zero. 

Eq. (3.1.2) shows that for CO, being different from zero the function 3 must be a func- 
tion of R or of Z or both. 

From considerations in Section 3.2, it seems that in the case of an axially symmetric 
flow there should be: & = 6. We choose: 


SG (RY) = (BZ). (4.1.1) 
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The relation 0 — 0 is used in Kgs. (3.1.3) and (3.1.4); the relation d = # (R, Z) is left 
in Kq. (3.1.2). 
With Eq. (4.1.1) and with p,, = 0, Eqs. (3.1.2) to (3.1.4) take the form: 


Cy = K* Qt he A (4.1.2) 


C, = K* (Roe) y,z; C,= — K? (Re) ye. (4.1.3) 


Starting with Eq. (2.2.17) one may use Eqs. (3.1.5), (3.1.6), (3.1.7) and (3.1.8), thus 
obtaining: 


ER 6 Ce = C.7 =; (4.1.4) 

KG, p= Cpe =— Cys (4.1.5) 

KG,z= yy — C, &; (4.1.6) 

E= Cyr — C23 9 = 7" (Cor), «3 

C= = 7G rib (4.1.7) 
Or: 

Gor = Cz (Cz,r — Cr, z) + Cy R (Cy R), ev; (4.1.8) 

G,2= Gk (CC, &), 2 C, (Car — Gp 2): (4.1.9) 

C, (Cy BR), rk = — C, (Cy RB), z. (4.1.10) 
Let us apply Eqs. (3.1.26) and (3.1.27) in conjunction with Eq. (3.1.2): 

Gi, = (CRO, 2 = Gy URI NC Gg Oy Te (4.1.11) 

G9 = (GZ zy.r — Gry, 2) (Cy Q)* K?. (4.1.12) 


One can insert Eqs. (4.1.8) and (4.1.9) into Eqs. (4.1.11) and (4.1.12). 
From Eqs. (1.3.6) one can see that for p,, = 0, ®,, = 0 (no forces acting and no 
pressure variation in the circumferential pees one has: 


é 
(C. se + Or gz) (RO) = F (RO) =0, (4.1.13) 


which is in agreement with Kq. (4.1.10). 

Let us assume that we found the equation of a streamsurface, y (Rk, Z) = const., 
which enables us to calculate the relation Z = g (R) valid on the surface. This allows 
us to express all the functions and all the derivatives in Eq. (4.1.13) in terms of one 
variable only, R, say. This immediately gives us a very important physical result, 
namely, that the angular momentum (RF C;) = 6 is constant on a streamsurface and 
is possibly only a function of (y, 9 = 0). We shall use this result. Namely, from [with 
9D = 9 (y, 9): 


do) = ORAR + 0% dZ, (4.1.14) 
one gets: ‘ 
ON) = OF Rey +02 Zr u3 OG = OF Reg 2 OU Zee, (4.1.15) 
and 
7 D (y, 8 
OR = (0, Z,9 — 0% Z, y) ae Ht (4.1.16) 
D (y, 8 
Oz = (0% Ry — OR, 9) Den af (4.1.17) 


Inserting Eqs. (4.1.16) and (4.1.17) into Eqs. (4.1.8) and (4.1.9) furnishes: 
Gor ae C, (C.R i: Cx, z) SiG bts C? Q herd (OZ, 9 4 Os Z, os (4.1.18) 
Gj2= RC? 0 K-05 Ry 0% Ris) = CC, p= Cy (4.1.19) 
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This inserted into Eqs. (4.1.11) and (4.1.12) furnishes: 
Gry = (C2 — Or,2) (C28, z — O, 9,2) (Cy 0) K? + 


(a) ae cal POO ea | 

+ R-* 9m D (yp, 8) O,z+ Dy, d) D, Ry; (4.1.20) 
Goa = — (Cir — Cz) (Cp, z+ C,p, x) (Co @)-? K*? — 

2 | D (0%, Z) D (0, R) | F 

pee R ee) in (y. 8) Y,2z Deg 3) Y, Ry E (4.1.21) 
Let F = F (y, 0) be any function; then: 

dF = F,, dy + F,y dd; en ee 
Prre=Fjyyrt Fie OR: Fiz=F,yy,z2+ F,9 9,z. (4.1.23) 


This furnishes: 


Fy = (Fr 9,2 — F,z9,r) _~ 3th ; (4.1.24) 
Fy = (Fiz rx — Fr yz) ne men (4.1.25) 
with D (y, 3)/D (R, Z) given by Kq. (3.1.2). 
Let: 
F=R; Ry =0,z K* (Cho); Bis = — y,z K* (CO, @)-*5 (4.1.26) 
== A; = — 0," K2(Cy 0); 2,9 =v, x K* (Cy 0). (4.1.27) 


With these equations one gets: 


D (0, Z) D (6), R) we - | vA 
Diy. 8) e+ Dew PF = [Fy edz — yz bn) 


— 03 (— dn 9,z + 8,2 9,n)] K* (9 Oy)! =-6); (4.1.28) 


D (60), Z) _ D(A), R) 
Dy) ?*T Diy, 


YR = [ay (WR UW aor YZ YR) =F 


— 03 (y, 29, z — y,z 8, r)] K? (Cy g)-? = — 8%. (4.1.29) 
Using Eqs. (4.1.28) and (4.1.29) in Eqs. (4.1.20) and (4.1.21) furnishes: 
G,y = RK? (6® 9) € (C, 8,z — C, 8,r) + ROM 6; (4.1.30) 
G9 = — RE (0 9) +E (C.y.zt Cyr) + ROM OY; (4.1.31) 
€ = "CL p— Cp 27 0 = RO. (4.1.32) 


Kgs. (4.1.30) and (4.1.31) are particular cases of Eqs. (3.2.13) and e 
one has to consider Eq. (4.1.10); with the use of Eqs. (4.1.16), (4. 
(4.1.27), Eq. (4.1.10) takes the form: 


a (Cyr + Oly, 2) + OY (0, 8,2 + 0.8,2) = 0. (4.1.33) 


Using Eqs. (4.1.3) one sees that the first term in this equation is equal identically to 
zero. To make Eq. (4.1.33) satisfied, one has two alternatives: 


14). Moreover, 
); 


3.2. 
1.17), (4.1.26) and 


(i) 6° = (RC;,),9 = 9, which would imply that 0 cannot be a function of #; 
this is true, since 7? = 0; 
(ii) (C, 9,r+C-9,z) =0; using Eqs. (4.1.3) this may be fulfilled with 0 = y; but 
in this case C, = 0 [see Eq. (4.1.21)]. Consequently, this must be excluded. Return- 
Sz 
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ing to Kq. (4.1.31), one sees that with G,» = G,,= 0, and 6% = 0, and with the use 
of Eqs. (4.1.3), this equation is identically satisfied and has to be excluded from the 
system in question. Consequently, 0 = R OC, may be only a function of y. This seems 
to imply that the function C,, Eq. (3.1.2) should be also only a function of y. There 
is also another item, which will be discussed below. Consider the first term in Kq. (4.1.30). 
Apply the following reasoning: Let the coefficient in the first term: 


R K? (6 0)- (C, 8, z — C, B, z), (4.1.34) 
be equal to —K 2(R o)-1. This implies that the following equality should be satisfied: 
Cpt (Cb 2= 0, 6) = (4.1.35) 
Using Eqs. (3.1.2) and (4.1.3) one obtains: 
— K? (Re) (y29,2+ 4,28, x) = — K? (Ro)? (y.r9z—y,z9,n). (4.1.36) 
This implies that there must be: 
—9y,z9,r= +y,z9,r; or: & = 8 (4), (4.1.37) 
or 
o = 0 (0, Z) only. (4.1.38) 


The function # which satisfies Eq. (4.1.38), satisfies the condition that C,, Eq. (3.1.2), 
is only a function of y, and allows one to pass from Eqs. (3.1.3), (3.1.4), to Eqs. (4.1.3), 
respectively, is of the form: 


0=04+Z. (4.1.39) 
Then Eq. (4.1.30) transforms into the form: 
Gy = — K* (Re) E+ ROM OY, (4.1.40) 


This is the only equation left in the dynamic system in question under the assumed 
conditions. It is identical to Eq. (16b) in ! and the entire procedure was constructed 
so as to make this equation identical to that in 1. The velocity component C,, Eq. (4.1.2), 
is now equal to: 

C, = K® 0-1 y,p. (4.1.41) 


Calculating y, x from the second equation (4.1.3) and inserting this result into Eq. (4.1.41) 
furnishes: 
Op SRG. (4.1.42) 


This implies that, if on a streamsurface, y = const., the magnitude RO, = 6 is 
constant, then also the magnitude (— R? C,) is constant on y = const. 


4.2. Particular Solution 


In order to demonstrate a possible method of solution of Eq. (4.1.40), an example 
will be solved. This obviously may be solved by an iteration procedure and in the first 
iterative approximation one may assume: 


Gy = 03. @ = const. (4.2.1) 


In this case, the problem reduces to the problem treated by Bragg and Hawthorne 
in }, They treat the equation (in their notation): 


Gy = 0 =F? (0.0; 9 sae ne i= tn eee (4.2.2) 
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Below, we shall briefly outline the main points of their reasoning : 
(1) they assume: ry =f (p)s P= Yot Dd) yn exp (+ kz); (4.2.3) 
where y,, and y,’ are functions of r alone; &, is assumed positive, and the sign of the 


index ( + k,) is such that the terms in the summation all vanish at great distances from 
the origin. With | z| > oo, the summation terms vanish, leaving: 


EW) = Poorr — 1! Woo,r3 (4.2.4) 


(2) let r, and 7, denote the radius at the tip and at the hub, respectively; let c,, denote 
the mean axial velocity of the flow, which is constant along the annulus in incom- 


pressible flow; then with (7 c,) = — y,,, one gets: 
“tip “tip 
Ga | ORS Or = | 2076, ary (4.2.5) 
hub hub 
tip 
C7 2y) (r/? a reese (4.2.6) 
|hub 
(3) Putting 
L(y) =— wy 17; w= const., (4.2.7) 
into Eq. (4.2.4), furnishes a Bessel equation: 
Woo tee as Woo r + we Lo aa Woo = 0, (4.2.8) 
with the solution: 
Poo ("? Cm)? = — R(ap-1) A, (u R), (4.2.9) 
R= rr; Ap (wR) = Jp (wR) + (a! a) Vy (pe R); (4.2.10) 
a and a’ are constants; a’ is arbitrary and: 
a = (1 — Ry) [Ay (u) — By Ay (ue BI, (4.2.11) 
with Rk, = 77,1. Each of the disturbing functions has a solution: 
? Gar tags == B(B A=") Bi AR); (4.2.12) 
ag (Ade) lg (Ack) (6 O41) Y fA Bh) AM = er? 2 (4.2.13) 


where k = k, (without subscript), and b, b’ are constants of integration. 
(4) From Eq. (4.2.2) one has: 
OD) 0,0 = = er = pn? yp, (4.2.14) - 
and: 
Ee Oe Seay ia he al rca (4.2.15) 
with g being another nondimensional constant of integration; another form is: 
(Cy Cm)? = (wu R-)? (72 Cn)? yp? +g R? = g R24 
+ fa A, (u R)+ D>) b, An B, (A, BR) exp (+ &,2)F. (4.2.16) 


The constants, appearing in the system, are calculated from the matching conditions 
at infinity, at the annular wails, and at the blade rows. 

As the second solution in the case of an incompressible flow, Bragg and Haw- 
thorne solve the case with G,,, + 0, and r 7 being a function of r only. 

Obviously, when g + const., one can apply the iteration process to calculate the 


higher order approximations. 


118 M. Z%. v. Kraywoblocki: Compressibility Effects in Circumferential Inlet Distortion 


4.3. Outline of a General Case 
In a general case we consider the system of equations (3.2.13) and (3.2.14): 


G,, =REE 0)-1 € (C. 9, z — O, 8, r) — 


— R- ae RK? (C,,59,2 + Cao 9, xn) + R= OM 64); (4.3.1) 
G, 9 = — RB K* (6 0) E(C,y,r — Czy, 2) + 
4+ R- 9-1 K? (C,,99,2 — Coops r) + R26 O°. (4.3.2) 


Moreover, one has to consider the following additional equations: 
Eggs. (3.1.2) to (3.1.4) for {Cy, C,, Cz}; 


Eqs. (3.2.4) for &,n, ¢; 
Kgs. (3.2.15) for the density o. 


We consider the system of equations (4.1.3) and (4.3.2) as a system in two unknown 
functions, y (6, R,Z), and 3 (6, Rk, Z) with the function G (y, 3) being known. One may 
attempt to find a solution of this system by some sort of an iteration process. In the 
first approximation one can assume that G, , = G, » = 0, that the flow is axially sym- 
metric and that consequently one can use the technique explained in Section 4.2. 
Next, assuming the knowledge of the function G (y, #) one can propose the following 
iteration procedure (this is an illustrative example, and other iterative procedure can 
be proposed, as well): 


ow pe) a oO”) eee a (R K2)-1 GL) (n) o(”) (&))-1 - [oe a (R o"))-1 K2 (on Bry + 
+ C®, 0) + R-2 9G). @ 90) 7; (4.3.3) 


ow) prt Lise ow yy ) = (R K%)1 99 ow) (E)-1 apat Gin a3 

+ (R09) K2 (Oy — Oy) + R*9O HO], (4.3.4) 
where: 
00 =6-+ Z,y given in Section 4.2. The functions 0”, m > 0, may be assumed in 
the forms analogous to the form y in Section 4.2. It is clear that in each particular case, 


depending upon the form of the function G (y, #), one may apply some particular 
iteration process and particular refinements. 


Appendix 
Invariants 
Gradient: V = a; 0/0 x; (A-1) 
Oa; = g"' djaa*: (A-la) 
aA, 2A aA, Ay . OAW 
iat Ae = 97> “(ZS oe = “s)4 a, +(24— 9) a, + (Be AL (A-2) 
chy me g- 0/da* (g% A*); (A-3) 
Vector Product: D. = Rs x - (A-4 


) 
p= g” (git Jd + gi? J? + gi8 J3); (A-4a) 
p? = g” (9? J} + g% jee g?8 Jy: (A-4b) 
p= g” (931 Jt + g?? J? + g?® J); (A-4c¢) 

) 


1 — 2 3 3 . as . 
Jt = 7? gh — 73 32; J? = 73 gl — yl g8, J3 = pl 52 — 72 gl, (A-5 
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Scalar products: 


A B= 94 Bi =A, Bis (A-6) 
(7) C= 17 0-0) = Ox one (A-7) 
curl C = ee 
PiCPOp= 1G) ao Kt R050; % | 
OK? Of Onih:: (A-8) 
ag: K? Cy Cy, 8; | 
O-Ve: x a K? R-2 0; Ch: 
we 6 48 | (A-9) 
On he Cee a, ==. 0, 
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Die mit Saulen unterstiitzte vertikale Kreiszylinderschale 
unter Windlast 


Von S. Sarkadi Szab6, Budapest 
Mit 9 Textabbildungen 


Zusammentassung. Die vertikale Kreiszylinderschale ist durch Windkrafte belastet. Die Unter- 
kante ist in Aquidistanten Punkten (Abb. 7) untersttiitzt. Die Oberkante ist frei oder ringformig ver- 
steift. Zufolge der unstetigen Stiitzkraftverteilung entstehen Biegespannungen in der Nahe der Unter- 
kante. Die in den Randbedingungen auftretenden Funktionen haben die Periode 22. Sie erleiden 


2 
jedoch Spriinge infolge der Stiitzen, deren Abstande = - sind. Der Verfasser schlagt eine Methode 


vor, womit man die trigonometrische Reihe einer solchen Funktion als die Summe von vier Fourier- 
Reihen anschreiben kann. Damit lassen sich die Biegespannungen in gewohnter Weise berechnen. 

Summary. The vertical cylindrical shell is loaded by windforces. The lower edge is supported 
at equidistant points (Fig. 7). The upper edge is free or carries ringshaped stiffeners. Due to the 
discontinuous distribution of the supporting forces, bending stresses arise near the lower edge. 
The functions occurring in the boundary conditions are periodic with period 22. How- 


20 
ever, they suffer jumps at the supports, whose distances are ae Author suggests a method for 


the determination of the trigonometric series of such a function as the sum of four Fourier-series. 
This makes it possible to compute the bending stresses in the usual manner. 
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1. Einleitung 


Nachfolgend wird die in Abb. 7 dargestellte vertikale Kreiszylinderschale behandelt. 
Die Schale ist an m aquidistanten Stellen unterstiitzt und durch Windlast belastet. 
Die Unterstiitzung erstreckt sich auf eine Lange 2 der Unterkante und ist gelenkig 
drehbar. Die Windkraftverteilung wird folgendermaBen angenommen: 


X=Y=0,2=p,+ pi cosg + p, cos 29. (hay 


2. Ziel der Arbeit 


Die statische Berechnung, die in diesem Artikel vorgelegt wird, bezweckt die Er- 
mittlung der Schnittkrafte. Wegen der unstetigen Unterstiitzung kann man die Schnitt- 
krafte nicht nach der Membrantheorie allein ermitteln, weil Biegespannungen in der 
Nahe der Unterkante entstehen. Zu ihrer Berechnung sollen die Funktionen, die zu- 
folge der Randbedingung unstetig sind, in solche trigonometrische Reihen entwickelt 
werden, welche die Unstetigkeiten darstellen konnen. Eine fiir diesen Zweck geeignete 
Entwicklungsmethode wird im Abschnitt V mitgeteilt. Die Berechnung der Biege- 
spannungen kann dann in gewohnter Weise durchgefiihrt werden. Ein numerisches 
Beispiel wird in Abschnitt VII gegeben. 


3. Vergleich der vorgelegten Methode mit den aus der Literatur 
bekannten Methoden 


2c 2 2 bes , 
aL ant aL tf Mi tt aL Ahnliche Windkraftverteilun- 
cat gen hat Girkmann! behandelt; 
JINN Vow Vg~ er berechnet jedoch nur Membran- 
oo Va Yass spannungen, wie sie bei stetiger 
Unterstiitzung entstehen. Olsen? 
und Ohlig* untersuchten die 
Schnittkrafte der an m Stellen 
unterstiitzten Schale, aber nur fiir 
eine zentralsymmetrische, kon- 


stante Belastung. [hre Verfahren 
sind nur fiir eine solche Last an- 


; aa 
wendbar, deren Periode rere 


da sie unstetige Stiitzkrafte mit 


F 22 
der Periode sn on der Unterkante 


hervorrufen. In unserem Falle ist 
laa die Verteilung der Windkrafte 
(1.1) von der Periode 2, aber die 
Stellen der Unstetigkeiten (Stiit- 


zen) sind durch die Winkel k io 


(« = 1, 2,3, ...) gegeben, die ain 

ae ea der Literatur angegebenen Me- 

or thoden sind daher hier nicht 
Abb. 1. Randkrafte der stetig unterstiitzten Schale brauchbar. 


’Girkmann, K.: Flachentragwerke. Wien: Springer, 1956. 

2 Olsen, Olav: Kontinuerlige skall. Akademisk Trikningsentral. Oslo: 1951. 

3’ Ohlig: Raiumliche Tragwerke des Stahlbetonbehalterbaus. Beton und Stahlbetonban, 1953, 
Oktober. 
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Die vorliegende Arbeit gibt eine Methode an, mit der eine Funktion, welche die 
Periode 2 hat und dabei Unstetigkeiten an den Stellen k = aufweist (Abb. 1 b) 


in eine trigonometrische Reihe entwickelt werden kann. Diese Methode ermoglicht es, 
die Randbedingungen (4.2) so zu formulieren, daB sie mittels der Biegetheorie behandelt 
werden k6nnen. 


? 


Die Entwicklungsmethode ist fiir alle orthogonalen Funktionsreihen brauchbar, 
nicht nur fiir trigonometrische Reihen. Die in Rede stehende Kreiszylinderschale 
stellt nur ein Beispiel dar. 


4. Grundsatzliche Beziehungen 


Die Berechnung der Schnittkrafte der — durch Windkrafte belasteten — Kreis- 
zylinderschale nach der Membrantheorie ist wohl bekannt. Die Schale kann aber nur 
dann als Membranschale betrachtet werden, wenn ihre Stiitzkrafte langs der Unter- 
kante dieselben sind, wie sie sich aus der Membrantheorie ergeben. Diese sind natiirlich 
stetig, daher konnen unstetige Randbedingungen nach der Membrantheorie nicht 
erfiillt werden. Aus diesem Grunde 
muB man zur Berechnung der 
Schnittkrafte der hier unter- 
suchten Kreiszylinderschale (spa- 
ter der Einfachheit halber nur 
,schale‘ genannt) die Biege- 
theorie mit heranziehen. 

Bei der Lésung des Problems . 
wird in gewohnter Weise voraus- 
gesetzt, daB die Schnittkrafte in 
hinreichender Entfernung von der f 
Unterkante als Membrankrafte 
betrachtet werden konnen. Die 
Biegespannungen vermindern sich 
schnell mit wachsendem Abstand 
von der Unterkante. 

Nach der Biegetheorie greifen 
an der Unterkante der unbela- 
steten Schale solche Krafte m,1, 
Ni an, welche mit den Mem- 
brankraften 0, N,qo Zusammen 
die durch (4.2) ausgedriickten 
Randkrafte 7,, 2, ergeben: 


Abb. 2. Membrankrafte der Schale 


Ul Nxo te Nxt 

Nx = Nye cle Meet 

eh ite) soe (4.1) 

Ny = Ny 0 + 2,61 Abb. 3. Sehnittkrafte der Schale 
Ny = Nyo + Ny1.- 


Diese Superposition fiir n, langs der Unterkante ist in Abb. | dargestellt. Hierbei 
ergeben sich betrachtliche Schwierigkeiten, weil n, und n,, die folgenden Randbedin- 


gungen erfiillen sollen: 
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23 {0 im Feld 

aunt | Flache von nm,» an der Strecke der Stiitze 
P 0 im Feld ee) 
ere | Flache von 2,9. an der Strecke der Stiitze 

m, = 0 

= —_ — 1 OMx el 

Ou = 0, qe iG Sia a h chests ’ 


wo die Flachen von no bzw. m,) an der 2 L langen Strecke zu rechnen sind. 


Die Untersuchung der mit beliebigen Kraften belasteten elastischen Kreiszylinder- 
schale wurde schon von mehreren Verfassern durchgefiihrt. Jede Methode ist auf 
dem linearen Elastizitatsgesetz aufgebaut. Alle Verfasser beniitzen dieselben Gleich- 
gewichts- und Deformationsgleichungen. Die genauen Differentialgleichungen sind 
aber kaum lésbar, darum fiihrt fast jeder Verfasser Vernachlassigungen ein. Mac Na- 


e/KYA 


Z 
= 
= 
ER 
= 
= 
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Abb. 4. Die ringformig versteifte Membran- Abb. 5. Die oben freie Membranschale 
schale 


mee‘ verglich die fiinf bekanntesten Methoden und leitete die folgende Differential- 
gleichung fiir die elastische, nur am Rande belastete Kreiszylinderschale ab: 


glm+d(n—7 Swe 


= 0, (4.3) 


wo w die radiale Verzerrungskomponente der Schale bedeutet und die Poissonsche 


1 
Zahl u = | zu Null angenommen ist. Eine wohlbekannte und allgemein gebrauchte 
Léosung von (4.3) ist 


w= a C e*l2 cos p —. (4.4) 
b = 


* * McNamee: Existing Methods for the Analysis of Concrete Shell Roofs. Proc. Symp., London: 
954. ¥ 
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Zu jeder ganzen Zahl p gehért je ein im Gleichgewicht stehendes Kraftsystem lings 
der Unterkante. Die Summe dieser Randkrafte soll zu den Membranstiitzkraften addiert 
werden und soll die Randbedingungen (4.2) erfiillen. In Abb. 1 a—c ist diese Summierung 
dargestellt. Infolge der Randbedingungen miissen die schraffierten Flachen einander 
gleich sein. 

Die Formeln der Membranstiitzkrafte sind bei beiden in Abb. 4 und 5 dargestellten 
Schalen ahnlich aufgebaut. Diese Formeln lauten wie folgt: 


GO, My 0 = A,, sing + A, sin 29 (4.5) 
No = Ag, COSM + Ag. Cos 2 g, 


wo die Werte der Konstanten im Falle Abb. 3 


Ay = pil, A, = 2 p,l 
py 4p, I? 
Ay = Ja : Ay = ve 


5 I? 


gq t2teH 
Ay = p, |, Ay = 4h Pl 
gq t 2 (L+H) 
12 
; rams 
p,/? 3a? Pa he 
Ay = aa Ay = 4h my a 


qi t20+e) 


sind. Die Randbedingungen (4.2) sind durch unstetige Funktionen ausgedriickt, deshalb 
miissen sie durch eine Reihenentwicklung transformiert werden. Die transformierte 
Form soll ebenfalls trigonometrisch sein wie (4.3) selbst. 


5. Eine Art der Reihenentwicklung unstetiger Funktionen 


Die Randbedingungen (4.2) sind durch unstetige Funktionen ausgedriickt. Die 
Periode der Funktionen ist 22, der Abstand der gleichférmigen Unstetigkeiten ist 


aber — . Die gesuchten trigonometrischen Reihen sollen also aus Gliedern der Periode 


von ungefahr Pu bestehen. Der Fouriersche Entwicklungssatz ist zu diesem Zweck 


nicht brauchbar, weil der Abstand der Unstetigkeiten mit der Periode nicht iiber- 
einstimmt. 


Die Beziehung (4.6) kann auch in folgender Form angeschrieben werden: 
Nxo = (Ag, Cos py + Ags cos 2g): 1 (¢), 
wo 1 (m) den sogenannten Einheitssprung der Operatorenrechnung bedeutet 


ee g>9 
eG gy < 0. 


Diese Funktion kann als eine Membranstiitzkraft (entstanden durch den Einfluf 
irgendeiner dimensionslosen Last) aufgefaBt werden und kann in eine solche Funktion 
1* (py) umgestaltet werden, dal} 


a 0 — im Feld 
Bre Wee | Flache von ee an der Strecke der Stiitze 


NxO 


sei, wie es in Abb. 6 dargestellt ist. 
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5 Yh a 
Die Fourier-Reihe der unstetigen Funktion 1* (y) von der Periode — ist bekannt 


und lautet © 9 9 
SIM 7 NE 
1* (9) = » ie OU 
ak (5.1) 


= »’ Ani COS NM YP, Ri Lyn 
n=1 
Die Funktion (5.1) approximiert die in Abb. 6 eingezeichnete Kurve mit Ausnahme 
der Unstetigkeitsstellen zufriedenstellend. An den Unstetigkeitsstellen tiberspringt die 
Fourier-Reihe immer die darzustellende 
Funktion um ungefaéhr 9% (Gibbssches 


Vise 1) Phanomen). Dies beriihrt aber den Flachen- 
ee ae inhalt der Kurven nicht und soll daher hier 
—---- -—--— 4- unbeachtet bleiben. 

Mit Hilfe der Funktion 1* (pg) konnen 


die unstetigen Randfunktionen (4.2) auf 
/*(#) folgende Weise angeschrieben werden 


(C) 


salt+ 12) | B 


20) 2c 
2g as 2L 
| 
Abb. 6. Randkraft der symmetrisch oder durch Abb. 7. Die oben mit Ring versteifte und unten 
Winddruck belasteten Schale auf m Punkten unterstiitzte Kreiszylinderschale 


fia = (AG, cos p + Ag. cos 2p): 1* (vy) = 


= (Az, cosy + A, cos 2 ~) + ¥'dn1 cosnm@p = 


5.2 
= YF {Aa [cos (nm + 1) 9 + cos (nm — 1)¢] + ne 
+ Ags [cos (nm + 2)p + cos (nm — 2)q)}. 
Analog kann abgeleitet werden 
<0 = (A, sing + A, sin 29) - 1 (g) 
Mei = (Ay, sing + Aj, sin 29) -1*(y) = (5.3) 


rant 


= 25 {A,, [sin (n m + 1)» — sin (nm — 1)p] + 


+ Aj, [sin (nm + 2) p — sin (n m — 2) q}}. 
WOR? =e aoe SOs 
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Die trigonometrischen Reihen der unstetigen Funktionen (4.2) und der Abb. 1 b—ce 
konnen also als die Summe von je vier Fourier-Reihen berechnet werden. Diese sind 
schon solche Funktionsreihen, die sich nach der Biegetheorie in bekannter Weise 
behandeln lassen beziehungsweise mit (4.4) verglichen werden kénnen. Da die Fourier- 
Reihen konvergent sind, existieren die Summen (5.2) und (5.3). 

Die hier angefiihrte Entwicklungsmethode ist fiir die trigonometrische Reihen- 


entwicklung solcher Funktionen brauchbar, die eine Periode 2 2 haben und in k = 


(i = 1,273, ... m)-Stellen Unstetigkeiten aufweisen. In ahnlicher Weise ist sie aber 
auch fiir andere orthogonale Funktionen verwendbar. 


6. Die Berechnung der Schnittkriafte 


Nach der Substituierung von (4.4) in (4.3) ergibt sich nach Weglassen von kleinen 
Gliedern die folgende sogenannte charakteristische Gleichung: 


Az ? 


Vee ie) ty jae es 12a?’ (6.1) 
deren Wurzeln 

Api,2 = @ (%1 +2 Boi) = — Apso 

Ays,4 = @ (%p2 tt By2) = — Arg 


sind. Die Form der allgemeinen Lésung von (4.3) wird die folgende sein: 
eS Bell ar cos Boi a+ Woe sin Boi eye ess ate 
p 


+ (Ws Cos Bpo x + Woa sin Byo x) e~ °??*! cos pop. (6.2) 


Dieser Ausdruck enthalt vier Integrationskon- iy 


stanten, welche aus der Randbedingung (4.2) 


berechnet werden k6énnen. 


7. Numerisches Beispiel 


Die Schnittkrafte des in Abb. 8 dargestellten 
Kihlturms sollen berechnet werden. Die Ab- 


messungen sind: 2a = 36,00 m, ¢ = 0,15 m, 

1 =47,00m,m = 30, 2c—1,10m, ¢ = — = 0,292. 

Werkstoff: Stahlbeton mit H = ai t/em?. NV VV 
Belastung: Winddruck mit py) = — 35,5 kg/m’, Ss 


p, = 59,25 kg/m?, p, = 94,8 kg/m?. 
Die Membranspannungen sind folgende: 
Noo = 639 — 1065 cos y — 1705 cos 2p 


Ny~o = — (59,25 sin y + 189,6 sin 2g) x + 
+ 2785 sin gy + 5178 sin 2 


Nxo = (1,6458 cosp + 10,5333 cos 2 py) x? — 
— (154,705 cos y + 575,333 cos 2 ¢) 
x + 3635 cosg + 3773 cos2py (7.1) 


ee ei cosine or Pid Samal oY (7.2) Abb. 8. Der durch Winddruck belastete 
Nxo = 3635 cosp + 3773 cos 2¢—. Kuihlturm 
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Die unstetigen Randkrafte der Schale kénnen aus (5.2) und (5.3) berechnet werden. Ks 
seien nur die ersten drei Glieder n == 1, 2, 3 beriicksichtigt. Die Randbedingungen (4.2) 
lauten jetzt wie folgt: 


Nx1 = 3266 (cos 32m + cos 28 y) + 3146 (cos 31 y + cos 29 y) + 

+ 1985 (cos 62 m + cos 58 py) + 1912 (cos 61m + cos 59 y) + 

+ 520,6 (cos 92g + cos 88 y) + 501,6 (cos 91 y + cos 89 ¢g) (7.3) 
= 4482 (sin 32 y — sin 28 y) + 2411 (sin 31 m — sin 29¢) + 

+ 2724 (sin 62 p — sin 58 y) + 1465 (sin 61 y — sin 59 py) + 

+ 714,5 (sin 92 m — sin 88 p) + 384,3 (sin 91 g — sin 89 ¢) 


mM. = 0 
gh == 0. 


Die Berechnung der Biegespannungen kann in drei Schritten durchgefiihrt werden: 


a) Lésen der charakteristischen Gleichung (6.1), 


b) Befriedigung der Randbedingung (4.2) bzw. (7.3) vermittels 12 algebraischer 
Gleichungen. Die Ergebnisse sind W,, .. 


W pa, 


c) Berechnung der Schnittkrafte nach der Biegetheorie. Die nach a) und b) berech- 
neten Wurzeln sind in Tabelle 1 und 2 angefiihrt. Die nach der Biegetheorie berechneten 
Schnittkrafte haben die folgende Form: 


Nx = »' cos pp { e ™* (€,, cos By & + Cy, sin B, %) + 
+ e~ %* (c,5 cos By & + Cy, sin B, x) } 


Ny = > cos pp {6-"* (a: C08 Py 


+ e~ ™* (¢., cos B, & 
m, = 2 cos pp { e~ °** (cs, cos By x 
+ e~ “* (C33 Cos Be x 


Nig == >) cos pp {e- “* (¢,, cos 6, x 


+ e7 %* (C43 Cos B. x 


M9 = >) sin pp { e-™* (cs, cos B, x 


+ € “(6.4 cos By x - 


Tabelle 1. Wurzeln der charakteristischen Gleichung 


Cy Sin B, x) + 
Coy 8. Ba x), 

C3 Sin B, x) 4- 
Czq Sin By 2) } 


Cag Sin B, 2%) 4 
Cag Sin: By x) } 


Cs Sin B, x) + 
Cra Sin. By Z) } 


P Hy Xo By By 

28 1,959 1,159 0,503 0,298 
29 2,014 1,214 0,499 0,301 
31 2.125 1,324 0,493 0,307 
32 2,180 1,380 0,490 0,310 
58 3,623 2,822 0,450 0,351 
59 3,678 2,878 0,449 0,352 
61 3,789 2,990 0,448 0,353 
62 3,845 3,044 0,447 0,354 
88 5,289 4,489 0,433 0,368 
89 5,345 4,544 0,433 0,368 
91 5,456 4,655 0,432 0,368 
92 5,511 4,71] 0,431 0,368 


p WwW, es W, Ws 
a eee ca ee ec a et 
28 + 2343,5706 + 1556,8623 — 2850,0929 + 2544,0989 
29 + 1606,6736 + 1046,3484 — 1879,4102 + 1935,1035 
31 + 107,35351 — 67,317758 — 13,068955 +711,45226 
32 — 435,96054 —496,07911 + 626,67294 — 263,96795 
58 + 681,67307 + 583,28494 — 698,44874 + 722,65585 
59 + 477,26463 + 399,13842 — 486,82567 + 528,00474 
61 + 52,255945 + 9,7786076 — 48,435901 + 132,42593 
62 —111,32608 —141,86735 +119,71183 — 19,167365 
88 + 117,65922 + 106,91846 — 118,64100 + 122,54810 
89 -+ 82,696659 + 74,061956 — 83,186125 + 89,070499 
91 + 9,9778930 + 4,3232080 — 9,7761890 + 18,758424 
92 — 18,640245 — 23,382085 +.19,095148 — 8,7631250 


Tabelle 2. Wurzeln der Gleichgewichtsgleichungen 
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Die iibrigen Schnittkrafte — namlich n,1, d:, ¢ — sind nicht berechnet. 


Die Schnittkrafte wurden zu den Membrankraften (7.1) addiert und die Summen in 
Abb. 9 eingetragen. Die Membrankrafte sind mit punktierten Linien eingezeichnet. 
Interessant sind die mit den Drillingsmomenten kombinierten Randscherkrafte: 


;00m +4087 


Reet Ot 


ee 1 axe 
Meg = Rae g agr * 
Thre Werte ergeben sich zu 
0 : 1 2 
s 30” 30° 
.oo* O + 1367,68 + 7219,21 4+ 3846,27— 3065,61 0 
Negit 0 — 1308,58 + 22086,02 + 11767,48 — 9422,02 0 


59,10 + 29305,23 + 25613,75 — 12487,63 0 


9 +8196 


D +14 435 


#30070 


Abb. 9. Die Kurven der Schnittkraifte der Schale 


(Bingegangen am 13. Januar 1961) 
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Uber die Lésung des Anfangswertproblems 
yy t+Ay +eBy=0, yO) yy Oe 
Von H. Seholz, Wien 


Zusammenfassung. Unter Zugrundelegung eines Satzes von Hermite wird ein Kriterium bewie- 
sen, mit dem man aus der Kenntnis der Koeffizienten der im Titel genannten reellen Differential- 
gleichung und ihrer Anfangswerte entscheiden kann, ob die Lésungsfunktion eine Jacobische ellip- 
tische Funktion, eine zusammengesetzte Kreis- oder Hyperbelfunktion oder eine rationale Funktion 
ist. Ferner wird eine Methode entwickelt, die es erméglicht, das Anfangswertproblem der obigen 
Differentialgleichung fiir beliebige Anfangswerte zu lésen. Die SchluBformeln sind so gebaut, 
da® man, ohne ein elliptisches Integral auf Normalform transformieren zu miissen, in ganz einfacher 
Weise die reellen Bestimmungsstiicke der Loésungsfunktion berechnen kann. Die Methode eignet 
sich auch zur Berechnung von Integralen, die auf elliptische Integrale erster Gattung fihren bzw. 
zur Bestimmung ihrer Umkehrfunktion. 


Einleitung 


Hermite hat folgenden Satz bewiesen’: Besitzt die Gleichung P (w’, w) = 0,w = 
= w (z) keine beweglichen kritischen Punkte, so ist ihr Geschlecht 0 oder 1. In diesem 
Fall ist das Integral der Gleichung entweder eine rationale Funktion, oder es 1aBt sich 
rational durch Exponentialfunktionen oder elliptische Funktionen ausdriicken. Dieser 
Satz lat sich nun anwenden auf die Integration der Gleichung w’" = P (w), wo P (w) 
genau den Grad 2 m hat. Zerlegt man P (w) in Linearfaktoren und beriicksichtigt die 
Fuchsschen Bedingungen, so lassen sich, weil die obige Gleichung nur feste kritische 
Punkte haben darf, genau acht Normalformen angeben. Fiir m = 2 ergeben sich 
dann die folgenden fiinf Normalformen (a; + a; fiir 7 + 7). 


1. w’? = by (w — aj)? (w — a) 


2.0 =05-(0.s= ay)? 

3. w’ = by (w — a,) (w — a.) 

4. w? = by (w — a)? (w — ag) (w — as) 

5. w’® = by (w — a,) (w — a.) (w — as) (Ww — ay). 


Die Gleichungen (1) bis (4) haben das Geschlecht 0, waihrend die Gleichung (5) das 
Geschlecht 1 hat. Daher sind auf Grund des Satzes von Hermite die Integrale von 1. 
bis 4. entweder rationale Funktionen oder rational darstellbar durch Exponential. 
funktionen, wahrend sich das Integral von 5. rational durch elliptische Funktionen 
ausdriicken 1aBt. 

Diese Erkenntnisse kann man dazu beniitzen, die explizite Lésung des komplexen 
Anfangswertproblems 


w' + Aw+ 2 Bw? = 0, w (0) = wo, w' (0) = wo’, (Wo, wo’) ¥ (0, 0) 


mit beliebigen komplexen konstanten Koeffizienten A, B + 0 anzugeben. Einmalige 
Integration ergibt 
w= — B (ws + sw? — @) C = Integrationskonstante 

womit der obigen Gleichung eine der Formen 2. bis 5. zugeordnet werden kann. Fiihrt 
man die GréBe D? = (A + 2Bw,?)? + 4 Bw,” ein, so lat sich, wie in I bewiesen 
wird, ein Kriterium angeben, mit dessen Hilfe es moglich ist, allein aus der Kenntnis 
von A, B, wo, wy’ zu entscheiden, welcher der drei moglichen Funktionsklassen die 
Loésungsfunktion angehért. AnschlieBend wird die explizite Lésung des Anfangswert- 
problems im Komplexen angegeben. Nun tritt aber in mannigfachen technischen und 
physikalischen Problemen 1 % 1° 18, 14 22, 23 das Anfangswertproblem in der reellen 
Form y" + Ay + 2By? = 0, y (0) = yo, ¥' (0) = ¥'o, (Yo: Yo’) + (0, 0) mit y = y (z) 


Uber die Loésung eines Anfangswertproblems 129 


und beliebigen reellen konstanten Koeffizienten 4, B + 0 auf. In der Praxis ist man 
meist an reellen Lésungsfunktionen mit reellen Bestimmungsstiicken interessiert, so 
daf die weiteren Untersuchungen im Reellen gefiihrt werden. Dazu wird gezeigt, dab 
man sich auf D? = 0, also D reell, beschranken kann, da namlich durch eine Trans- 
formation der Fall D? < 0, also D imaginar, auf den vorigen zuriickgefiihrt werden 
kann. Das im Komplexen bewiesene Kriterium wird ins Reelle iibertragen und hat 
dann folgenden Wortlaut: 


Die Integration des Anfangswertproblems y’’ + Ay + 2By3 = 0, y (0) = yp, 
y' (0) = Yo’, (Yo: Yo ) + (0, 0) mit beliebig reellen konstanten Koeffizienten A, B + 0 
und D? = (A+ 2By,?)? + 4By, = 0 fihrt dann und nur dann 1. auf nicht aus- 
geartete einfache und zusammengesetzte Jacobische elliptische Funktionen, wenn 
0 <D+/A\|, 2. auf ausgeartete Jacobische elliptische Funktionen, namlich auf 
zusammengesetzte Kreis- und Hyperbelfunktionen, wenn D = 0, A +0 oder 
D =|A\| + 0, 3. auf eine rationale Funktion, wenn D = A = 0 ist. 


Damit sind die Klassen der Lésungsfunktionen festgelegt, und es dreht sich nun darum, 
diese Funktionen auch tatsichlich explizit anzugeben, und zwar in reeller Form. In II 
werden die drei Hauptfalle getrennt untersucht. Kernproblem ist die Lésung des Haupt- 
falles 1, also die Bestimmung derjenigen Jacobischen Funktionen, die das konkret 
vorgegebene Anfangswertproblem lésen. Es wird sich zeigen, dafs dies wesentlich davon 
abhangt, welche Vorzeichen A und B haben und ob D ¢ | A_| ist. Die Bestimmung 
der Jacobischen Funktionen, welche der Gleichung y’2? = — B (y' oy — 4 
= — B(y* — b,?) (y? — b,?) geniigen, wurde bereits von Briot und Bouquet®, aller- 
dings nur fiir ganz spezielle Anfangswerte, durchgefiihrt. Dazu wurde die rechte Seite 
der obigen Gleichung auf die den Jacobischen Funktionen eigentiimliche Form (1 — wu?) 
(1 — k? u?) mit 0 < k? < 1 transformiert. Wenn man aber allgemeinere Anfangswerte 
beriicksichtigen will, insbesondere den Fall y,.’ + 0, so ergeben sich bei der obigen 
Transformationsmethode Schwierigkeiten. Daher wird hier ein anderer Weg eingeschla- 
gen. In ihrer einfachsten Form wurde nachfolgend erlauterte Methode vom Verfasser 
auf ein ganz spezielles Beispiel aus der Theorie der Fahrzeugfederung in einer friiheren 
Arbeit!* angewendet. Die Methode besteht in folgendem: Bedeutet allgemein y = p q(x)” 
eine der zwolf Jacobischen Funktionen, so geniigt sie der Gleichung 


y’ + f(y + 29 (kh) y = 0 


mit speziellen Anfangswerten. Will man allgemeinere Anfangsbedingungen beriicksich- 
tigen, so verallgemeinert man y = pq (x) zuy=apq(wx-+y), woa, w, y Konstante 
sind. Diese Funktion geniigt dann der verallgemeinerten Gleichung 


2 w? 
y+ of (h*) y +g (ke)? = 0 


mit allgemeinen Anfangswerten. Bemerkenswerterweise tritt die konstante Amplitude 
a bereits in der Differentialgleichung auf, was bei linearen homogenen Gleichungen 
nicht der Fall ist. Soll nun y = apq(w a+ y) Lésung von 
y’ +Ay+2By=90 
sein, so mu gelten: 
2 
A= otf (2B ——, 9 (lc) 


; d 
Yo=4PT(Y), Yo = ala, Pg (ox piv )| Sele: 


Wie sich spater in II zeigen wird, sind von den zwolf Jacobischen Funktionen in 
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diesem Zusammenhang nur sieben wesentlich. Fiir diese sind f (k?) und g (k’) linear 
in k?, so daB 


A OP Gye ee k?), B Gs ee Og k*), 


wo die Zahlen «,, %», %, %, Charakteristisch sind, gesetzt werden kann. Aus diesen 
beiden Gleichungen lassen sich nun w? und k? durch a? ausdriicken. Da aber anderer- 
es AY ag Fa, he 
seits a =F ya, kB 
Schranken: s ( 5) ava (5). Damit hat man die allgemeine Lésung der gegebenen 
Differentialgleichung gefunden: 


und 0 < k? < 1sein soll, so ergeben sich fiir a? im allgemeinen 


y=apglo(?)a+y), P=’), ssas 8. 


Die noch freien Konstanten a und y werden durch die Anfangswerte festgelegt. Zu 
diesem Zweck wird die gegebene Differentialgleichung einmal integriert. Mit der GroBe 
D? ergibt sich fiir « = 0 die Identitat 


z ; D—A , D+A 
yo? = — By —sa) (ye + IB ), 
: : = D— A.) Da A ; : eee 
wo wegen D? > 0 die beiden GréBen >=, gp Stets reell sind. Andererseits gilt fiir 


die sieben Funktionen allgemein die Identitat 


Yo? = — By? — € 4”) (yo? + 4 (k?) a?), 
wo ¢ und 4 (k?) bezeichnend sind. Beide Identitaéten ergeben zusammen die weitere 


(yo? — €4*) (yo? - 207) = (y.’ is 7) (y.? si Fr): 

aus der wegen der Schranken s, S fiir festes « und / a? durch D, A, B bestimmt wird. 
Ferner wird dadurch festgelegt, ob D < | A|oder > | A| ist, woraus wieder Schranken 
fiir y, und y,’ folgen. Mit diesem so errechneten a? bestimmt man w? und k? durch 
D, A, B. Aus yy) = apq (y) ermittelt man | sn y| = sin y, setzt « = are sin k und erhalt 
y =F («,g) mit 0 sy < K. Je nach den Vorzeichen von y, und y,’ hat man fiir die 
einzelnen Funktionen y, y + 2 K,2 K — y, —y zu nehmen. In einem speziellen Fall 
mu man auch a <0 zulassen, wahrend sonst stets a > 0 ist. Damit erhalt man die 
gesuchte Losung des Anfangswertproblems fiir den Hauptfall 1: 


y =a(D, A, B) pq [wm (D, A, B) «4+ y), kb? = k* (D, A, B), y=F (a, go), « = are sink, 
y = arc sin (sn y). 


Hauptfall 2 erledigt sich durch Spezialisierung der vorigen Ergebnisse, wobei bemerkens- 
wert ist, daB in diesem Fall a? nicht mehr in gewissen Schranken liegt, sondern durch 
A und B festgelegt ist, wodurch y,’ bereits durch Vorgabe von y, bestimmt ist. Dadurch 
tritt der Hauptfall 2 und ebenso auch Hauptfall 3 nur fiir ganz bestimmte Konstel- 
lationen der Anfangswerte ein. 

Die Tabellen 7 und 9 fassen die Ergebnisse iibersichtlich zusammen und kénnen als 
praktische Berechnungstabellen beniitzt werden. 

III. bringt emige Anwendungsbeispiele. 


1. Beweis des Kriteriums. Die drei Hauptfalle 
In der komplexen Differentialgleichung 
w+ Aw + 2Bu® = 0, w= w (z) 


mit beliebigen konstanten komplexen Koeffizienten tritt die Verinderliche z nicht 
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explizit auf, so da man bekanntlich? ohne Beschrankung der Allgemeinheit die 
Anfangswerte an der Stelle z, = 0 vorgeben darf. 
Ks ist dann das Anfangswertproblem 


wr Awe 2bw = 0, w (0) = wo, wv’ (0) wy’, (Wy, W5') #1 (0; 0) (dea) 


zu losen. Die Spezialfalle A = B = 0 (Polynome), A + 0, B = 0 (einfache Exponen- 
tialfunktionen) bleiben auBer Betracht, so daf fiir das Folgende stets B+ 0 voraus- 
gesetzt ist. (1,1) ergibt einmal integriert die aquivalente Differentialgleichung 


Z C 
yee (ws ee 3) (1,2) 


(C Integrationskonstante), deren Integral auf Grund des Satzes von Hermite eine ratio- 
nale Funktion von elliptischen Funktionen oder Exponentialfunktionen oder der 
Veranderlichen z allein ist. Welcher der drei Funktionsklassen das Integral angehort, 
hangt vom Geschlecht’ der Gleichung (1,2) ab. Ist das Geschlecht 1, hat also (1,2) nur 
einfache Nullstellen, so tritt der erste Fall ein, ist das Geschlecht 0, hat also (1,2) mehr- 
fache Nullsteilen, so tritt der zweite und dritte Fall ein. Es handelt sich nun darum, 
aus der Kenntnis von A, B, wy, wy’ heraus zu entscheiden, wann (1,2) einfache bzw. 
mehrfache Nullstellen hat. Diese Entscheidung liefert das folgende 


Kritervum 

Die Integration des Anfangswertproblems w’’ + Aw + 2 Bw? = 0, w (0) = wW», 

w’ (0) = wy’, (Wo, Wo) + (0, 0) mit beliebigen konstanten komplexen Koeffizienten 

A,B+0 und D? = (A + 2 Bw,?)? + 4 Bw,? fiihrt und dann nur dann 

1. auf elliptische Funktionen, wenn 0 + D+ 4A, 

2. auf Exponentialfunktionen, wenn D = 0, A + 0 oder D= A # 0, 

3. auf eine rationale Funktion, wenn D = A = 0 ist. 
Zum Beweis driicke man zunachst die Integrationskonstante C' aus (1,2) durch die 
Anfangswerte aus. Es ist dann 

D? — A? 


Cc =w?tB Wy + Aw? = er (1,3) 


Setzt man (1,3) in (1,2) ein und zerlegt die rechte Seite in Linearfaktoren, so gilt mit 


bP = _ , 62 = 2B (1,4) 
w= — B(w* — 6,?) (w? —b,") = Bm w). (1,5) 
Fiir die Diskriminante //? des biquadratischen Polynoms p, (w) erhalt man 
4 Dt (A? — D? 
pe ee 41,8) 


Daher hat p,(w) dann und nur dann einfache Nullstellen, wenn 0 + D+ A_ baw. 
dann und nur dann mehrfache Nullstellen, wenn D = 0, 4+ 0 oder D= A+0 
oder D = A = 0 ist. (1,5) nimmt dann eine der drei folgenden Formen an: 


I w’? = — B(w? — b,?) (w? — b,?) 0 +b, +b, +0 
— A 
IT a) w' = //— B (w* — 6,3) NE anes ets ae Xe 
A 
b) w’? = — B (w? — b,?) w? 510, b= 5, +0 


b, =O, 
Ill w = /—-Bw bt Din, =; 0: 


Q* 
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Durch geeignete Transformationen’ erhalt man die folgenden Losungen: 


I Daa a 6D | 
el Cee 3D—2A+ 12) (2 + »). 
mit den Invarianten . 
2(A— Dt (D — A)|= 
pe = (G, 28, m+ 3m); g=4| 2B (m3 — s,m? + 58. m — 83), 
_24—3Dy,/ 2B 2B Bi 2 Bs 8, 


a op - DLA’ 2 ep? 8 ep y eee 3 
so) oJ feno|—e/3 +] 40 le [-2 V4 + ol] 


b) w= eee 2/lexp (: \—A+ ») 4 exp i= 2\— 4—,)\ 


1 


Mie (27 Su, (eS 


y ist jeweils durch wy, zu bestimmen. Damit ist das Kriterium bewiesen. 
Die weiteren Untersuchungen werden nun im Reellen gefiihrt. 
Gesucht ist die reelle Losung des Anfangswertproblems 


y' +Ay+2BY=0,y = y (x), y (0) = Yo ¥' (0) = Yo> (Yo Yo) + (9, 9) (1,7) 


mit beliebigen konstanten reellen Koeffizienten A, B+ 0. Die Bestimmungsstiicke 
der Lésungsfunktionen sollen ebenfalls reell sein. Mit 


D? = (A + 2 By3)? + 4 B yy? (1,8) 
gilt 
pee YG eae (1,9) 
D—A D+A 
6)? = 9B? 6.2 = — 2B (1,10) 
Zunachst soll bewiesen werden, dafi man sich auf reelle b,? und b,?, also reelles D, 
d.h. auf D? = 0, beschranken kann. Angenommen, es ist D? = — A? < 0, also D=iA, 
dann sind b,? und b,? konjugiert komplex. Dieser Fall kann nach (1,8) nur dann eintreten, 
2 
wenn B <0, also y)?> oS ist. Setzt man 
A—iaA ; A+iaA = : 
bf? = SS’, b.? = 402) oO. =e O ae te (1,11) 
so fiihrt die Transformation 
— l+4u y? — 00 
Y= 00 jy Y= Tes (1,12) 
die Gleichungen (1,9) und (1,7) in 
Ay ee G2UE 2 ( raga 
= — 4 Bo*(1—w) (w+) (1,13) 
wu’ + 4B (o? — 77)u — 8 Bo? w= 0 (1,14) 


iiber. In (1,18) sind die den 6,” und 6,? entsprechenden GréBen wieder reell. Es geniigt 
also, im folgenden stets D? = 0 vorauszusetzen. Dann aber hat (1,9) fiir 0 + b, + 6, +0 
die den einfachen und zusammengesetzten Jacobischen Funktionen eigentiimliche 
Form‘. Ubertragt man nun das im Komplexen bewiesene Kriterium ins Reelle, 
so hat es mit den iiblichen Modifikationen den folgenden Wortlaut: 
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Kritervcum 
Die Integration des Anfangswertproblems 
y+ Ay + 2 By? = 0, y (0) = Yo. y' (0) = Yo's (Yor Yo) (0, 9), 
A,B+0 beliebig reell, D? = (A +2 By,2)??14By,"%= 0 


fiihrt dann, und nur dann 


1. auf nicht ausgeartete einfache und zusammengesetzte Jacobische elliptische 
Funktionen, wenn 0< D +/ A), 
2. auf ausgeartete Jacobische elliptische Funktionen, namlich auf zusammenge- 
setzte Kreis- und Hyperbelfunktionen, wenn D = 0, A + 0 oder D =| A! + 0, 
3. auf eine rationale Funktion, wenn D = A =: 0 ist. 
Die im obigen Kriterium genannten Bedingungen legen zunachst nur die Klassen der 
Lésungsfunktionen fest. Um zu erreichen, daB diese und ihre Bestimmungsstiicke reell 
werden, sind weitere Bedingungen zu erfiillen. Diese werden in II fiir die obigen drei 
Klassen getrennt entwickelt. 


2. Untersuchung der drei Hauptfalle 
Haupijfall 1: 0< D+ | A} 


Ks gilt zunachst nach (1,9) die Identitat 
i psSA D+A 
Siar 8 GED EY) (pol g 2) 
Als Integrale treten Jacobische elliptische Funktionen auf. Sie werden hier als reelle 
Funktionen mit der reellen Periode 4 K, wo K das vollstandige elliptische Normal- 
integral erster Gattung bedeutet, betrachtet und sollen folgendermaBen bezeichnet 
werden* 15 24 


(2,1) 


sn x= sinam x Cn x = cosam «x dn x = 1 — k? sn? x 
1 1 Fl 1 
ae 2 SS SS nd x= -— 
Uae Snax bec nar dna 
j Sux Cnet Ti dna (2.2) 
LE ey OE EE Ser aa pe 
sn x Cne dnx 
i SSeS ; — dcx = — 
sd x dna cd x dn x CN 


mit 0<k?<1lundx«—F (4,¢),¢ =am2,« = aresink. Bezeichnet man allgemein 
mit y = pq (x)}* eine der obigen 12 Funktionen, so geniigt sie einer Gleichung vom 
Typ (1,7), namlich 

y +h )y +29 () y= 90 (2,3) 


mit ganz speziellen Anfangswerten. Zum Beispiel gentigt y= sna der Gleichung 
Yo t By — 2 yt = 0, Yo. = 0,49 = 1. 


Um nun (1,7) fiir allgemeinere Anfangswerte lésen zu kénnen, ist es zweckmaBig, die 
12 Funktionen durch Hinzunahme freiverfiigbarer Konstanten zu verallgemeinern. 
Naheliegenderweise wird man y = pq (x) zu y=apq (wa + y), wo jetzt a, w, 7 
konstant sind, verallgemeinern. Diese Funktion geniigt dann der Gleichung 


20? Bas ) 
yf" =a w? f (k?) y+ ag (k ) y* = 0, (2,4) 
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wie sich durch geeignete Transformation von (2,3) leicht ergibt. Da nun (2,4) den 
Modul &2 als Parameter enthalt und iiberdies die konstante Amplitude a bereits in 
der Differentialgleichung auftritt, was ja bei linearen homogenen Gleichungen nicht 
der Fall ist, werden sich fiir a? im allgemeinen Schranken ergeben. Das bedeutet 
andererseits, da auch den Anfangswerten im allgemeinen Schranken auferlegt 
werden, innerhalb deren sie dann allerdings willkiirlich vorgegeben werden kénnen. 
Spezialisiert man nun (2,4) fiir die 12 Funktionen (2,2) und vergleicht mit der 
Grundgleichung (1,7), so ergeben sich mit 


A a w* f (k*), B =" (k?), Ox ao y == TEE 


die folgenden Normaldifferentialgleichungen 


Tabelle 1 
Nr. y= Y) (a) Aw? a? Bia 
I acu 1 2h = 6 hee 
II adnu —(2—k*?) < 0 ] 
III asn u 1 + k2=> 0 = ie 
IV ans U l+k>0 = 
V anew 12 2E2 = 6 (lb ==) 
VI asc u —(2 — k?) < 0 —(1 — k*) 
VII acs u —(2— k*?) <0 — ] 
VIII adeu 1+ kh? >0 -- J 
TX asd u p= 22.0 k? (1 — k?) 
X and u —(2 — k?) < 0 1] — k? 
XI ads u l= 2k 20 — ] 
XII acd u 1+kh>0 — kh? 


Von diesen 12 Funktionen kommen zur Lésung von (1,7) aber nur die ersten sechs als 
voneinander wesentlich verschieden in Frage, da sich namlich die Funktionen VII 
bis XII nur durch Vermehrung von u um K, 2 K, 3 K aus den vorhergehenden Funk- 
tionen ergeben, wie die folgende Tabelle 2 }° zeigt. 


Tabelle 2 


bat oe Ke 6 VAD ON wah clea ie 
sn cd wu — snu —cdu 
cn — k’ sd u —cenu k’ sd u 
dn k’ ndu dn wu k’ ndu 


Aus Kontinuitatsgriinden, die sich spater ergeben werden, wird im folgenden auch noch 
VII mitgenommen. Ordnet man die Funktionen I bis VII der Tabelle 1 nach den Vor- 
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zeichen von A und B, so zeigt sich, daB bestimmten Vorzeichenkombinationen von A 
und B nur ganz bestimmte Lésungsfunktionen zugeordnet werden kénnen (Tabelle 3). 


Tabelle 3 
A B y = y (2) 
> 0 =O acnu 
== 10) > 0 acnu,adnw 
= 0 <a asnu,ancu,ansu 
== (0) <0) anew 
<0 = 0 asCU,ane u, acs WU 


Diese Zuordnung ist aber im allgemeinen zunachst nicht eindeutig. Es sind also zusatz- 
liche Bedingungen notwendig, die eine eindeutige Zuordnung erméglichen. Solche Be- 
dingungen ergeben sich dadurch, da8 man a? unter der Voraussetzung 0 < k? < 1 
abschatzt. Dazu betrachtet man die den Funktionen I bis VII der Tabelle 1 entsprechen- 


den Koeffizienten A und B. Man erkennt, daB 22 bzw. a linear in k? sind. Man kann 


also setzen 


2 
A = ow? (%, + «2 k?), B= a (a + o%4 &?), (2,5) 
woraus sofort folgt 
A +a,k 
Pe 20) 
Aa, — «a Ba, Oy Oy 2 a Ba,—A as 27) 
eo = d ea Og Oly > — A oy — a? Bay (2a 


Die Zahlen «,, «5, %3, %, haben fiir die Funktionen von Tabelle 3 die folgenden charak- 
teristischen Werte: 


Tabelle 4 
Nr y = y (2) oy Oy ag Ly d 
I acn uw : 1 =—2 0 1 1 
LE adnu = 1 1 0 —1] 
III a sn U 1 1 ) ==] —] 
IV ans u 1 1 al 0 1 
V ane u 1 =) —| 1 al 
VI asc u —2 ] =] ] =) 
VII acs wW —2 ] —] 0 Ud 


Setzt man nun in (2,5) die speziellen Werte aus der obigen Tabelle ein, so kann man a? 
abschatzen und erhalt folgende Schranken: 
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Tabelle 5 
a ee 
Nr. (cain Gan Chay CA) a Schranken fiir a? 
A>0O | pe et 
3 : A k? = Hl 
1 (1; — 2, 0, ) ip [= Oe < 0 a2 Sof os 
| B 
| 1 —A : 
1 (—25 1, 1, 0) Taig © 2 — fe 9B <a Sa 
UI 1 108 9) A = Vea = 
ee a = Boe 2B 
A 1 A ope! 
IV (1, 1,1, 0) ae ae Se ee 
2 : | 7 
A yee Asc Dey) Raeimarogy 
V (1, —2, —1, 1) mas 12 | r 
: | di 0) plhig] gee een 
| 4 1—K , 
VI (—2, 1, —2, 1) ae eae Oraa <oR 
| { 1 A A 
] | SO =o 10 eee ee eee 
VII ( y 4 sooo oh. an 


Bei I und V sind die Falle A> 0, also 0 < k?<4, und A <0, also $< k? <1, zu 

ne: : ah ke 5 - L— 
unterscheiden, wobei sich zeigt, daB bei I fiir A > 0 |— gq baw. bei V fir A < 0 542-4 
beliebige positive Werte annimmt, so daf also a? in diesen beiden Fallen nicht beschrankt 
ist. Ferner ist auch noch der Fall A = 0, also k? = } zu untersuchen, doch ist dies 


in der obigen Form nicht unmittelbar méglich, da sich fiir a? der sinnlose Ausdruck 


ergibt, dem erst im Zuge der weiteren Untersuchungen ein Sinn beigelegt werden kann. 

Berechnet man nun nach (2,6) mit den Zahlen «,, «., «3, «, der Tabelle 4 die den 
Funktionen I bis VII entsprechenden Werte von w? und k? und nimmt noch die Schran- 
ken von a? aus Tabelle 5 hinzu, so hat man damit die allgemeine Lésung der Grund- 
gleichung vor sich. ,,Allgemein“ bedeutet hier, da die Konstanten a und y durch in 
gewissen Schranken beliebig vorgegebene Anfangswerte bestimmt werden kénnen. 
Daf man die Anfangswerte nicht ganz beliebig vorgeben kann, wie das bei der allge- 
meinen Loésung linearer homogener Differentialgleichungen mit konstanten Koeffi- 
zienten moglich ist, liegt einmal an der Nichtlinearitat der Gleichung, wodurch die 
Amplitude a bereits in der Gleichung auftritt, zum anderen daran, daB der Parameter 
k? zwischen 0 und 1 liegen soll, weswegen a dann eben im allgemeinen beschrankt ist. 

Fiir I aus Tabelle 5 soll nun der Vorgang zur Bestimmung der allgemeinen Lésung 
erlautert werden. 

Hier ist «, = 1, «, = — 2, ag = 0, a, = 1,d = 1, nach (2,7) also @2 = 4 4 2a 8 

2 

es ser Fir A > 0, also 0 < k? < }, ergeben sich keine Schranken fiir a2, also 
lautet fir A > 0, B> 0 die allgemeine Lésung y=acn (wx+y), w= A+ 2a? B, 


2b 
i? = 45 eae Fir A = — | A |< 0, also } <k? <1, muB a? > be sein. Daher erhalt 
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: Aika : 
man fir 4<0, B>0, a> =! die allgemeine Lésung y =acn (wx + y), 


2 
ot = —|A| + 2a*B,k? = — {7S 5. Num lat sich auch der Fall A =0 be- 
2 2 
handeln. A = 0 bedeutet k? = } und es folgt daher w? = 2a? B,4 =k? = oe = sa : 


also mu a?+ 0 und endlich sein. Man erhalt daher als allgemeine Lésung fiir 
A=0, B>0 y=acn(wux+y), w? = 2a? B, ae 

Ahnlich geht man bei II bis VII vor. 

Zur Lésung des Anfangswertproblems driickt man a? durch D, A und Baus. Hier geht 
die Voraussetzung 0< D+) A) und die Beschranktheit von a? ein. 

Fir die Jacobischen elliptischen Funktionen gilt in einfacher Verallgemeinerung 
bekannter Resultate * 2 


Gy ae UY) 6") (Ye A(R?) a"), (2,8) 
wo ¢ und 4 fiir die Funktionen I bis VII die nachstehenden Werte haben: 
Tabelle 6 
Nr. I | ii Ill HW V | VI VII 
€ 1 I 1 1 1 S| —] 
1 a k2 ae ow } ~. i he2 1 ] od ke 
ke ta Ke 1_ Fk 1 — # 
(2,1) und (2,8) ergeben die weitere Identitat 
DA Dawa : 
(yo? — © a2) (yo? + Aa’) =(yo?— sa) (ye + a) (2,9) 


Da nun entweder 0< D<|A| oder D> | Aj gilt, dies aber von den Anfangs- 
werten und den gegenseitigen Vorzeichen von A und B abhangt, ergeben sich Schran- 


A : 
ken fiir die Anfangswerte. Es ist klar, daB fiir y.>+ —,,, B> 0 stets D? > 0 ist. 
Ist aber B < 0, dann mu8, damit D? > 0 ist, 


A-+2B y,?)? A 
i — aL Ys y+ IR (2,10) 
sein. Wie leicht zu beweisen, ist fiir 
1A>0,B>0 stets D>A (2,11) 


2ew = OB = 0 
a) D>|A|, falls y.?> — y,.*(A + By,*) 
by0< D=|A), falls 0< yo? =<.— 9° (A + B yo’) 
—A 
und gleichzeitig yo? < —- 


} (A +2 B y)?)? 
a A> 0: Bie 0202 Map <p sae as 


a) Ds A, falls 0< Yous = oe (A a B y?) 
A 
und gleichzeitig y.” > —; 
b).0< D= A fiir 


A , (A +2 By’)? 
a) 0< Yor = Sop falls O< iy Mir a 
Aye tesa ae 7 falls 
; (A + 2 B yy’)? 
ye (AB fe) a = ap 
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4,.4=0,B<0 D> 0, falls’ y,2=< — B y,* 
2B 2)\2 
oF Ax ®) Bay, Oe 
a) D>|A|, falls O< yo? < — yo? (A + B yp’) 
Dp Oe Dies A|, falls 
1g a Ee Baye 
= Ye (A BY) = Ye" = =a 


Nun 1a8t sich a? fiir die einzelnen Funktionen I bis VII durch D, A und B ausdriicken. 
Da « und 4 charakteristisch sind fiir diese Funktionen, so entspricht jeder eine ganz 
bestimmte Form der Identitat (2,9), die nun entsprechend den Schranken von a? ent- 
weder durch 0< D</|A| oder D>j|A| befriedigt wird. Beachtet man, da8 zur 
Berechnung von y fiir I, II, III y,?< a* und fiir IV, V y,?= a? sein mu, so erhalt 
man mit (2,11) die entsprechenden Schranken fiir die Anfangswerte. Mit dem so be- 
stimmten a? lassen sich auch w?, k? und y durch die Anfangswerte ausdriicken. y wird 
so berechnet, dai man jeweils | sny|—=singy und hierauf y = F (a, p) mit «= 


=are sink, 0s as 5 bestimmt und nun die 4 K-Periodizitaét ausniitzt, um den ge- 


gebenen Vorzeichen von y, und y,’ gerecht zu werden. 
Bei IV etwa geht man so vor: Hier ist d>0, B<0, y=ans(wx+y), «4 =1, 
A 


y. 6s 
te, lL, a= — Lo, = 0, dl, also 02> —o7 bir = ag ; 97 = Se 
bie Da A D+A 
é=1, A= — k*., (2,9) lautet dann (y,? — a?) (y,? — k? a?) = (yo? IB ) (yo? +ee 
DAs 
also kann nur a? = =s mit 0< D<A gelten. Da y,?=a?> —jq, so folgt 
A+2 B y?)? E +A 
aus (2,11) 3.b) 8) — y,2 (A + B yo?) < yo? < ( “Be 4 Be . Mit ot _ 
3 D+A A—D j 
ist dann w? = 5 : ee ay | sn y | =sing=p0- Aus ¥, = ans y, 
Yo =—aons*ycenydny folgt die Zuordnung 
] / 
eee > 0 z= 0 
>0 2K —y y 
: ” A . A+2By,?)? 5. 
Somit autos turA>0, B=0rg= ep — Yo (A+ By*)<yo?< fee die Losung 
_ VRtA4 (e /D+A Su. d= Ds 1 1/D+A ee ea! 
y= ap” x| O) ey Ah — SEeE De ego seiag a Fi =o? a = are sin |/ 7, , 


Se (ase): 
Ahnlich geht man bei den anderen Funktionen vor. Bei II, also y=adn(wx+y), 


mu man, um auch yy < 0 erfassen zu kénnen, a < 0 zulassen, da ja bekanntlich dn u 
fiir reelles w stets positiv ist. Bei I und V tritt der Fall A = 0 auf, der jetzt ohne Schwie- 


rigkeiten erledigt werden kann. Fiir A + 0 ist bei I a? = Z ee baw. bei Va? = 2t4 


Wie man sieht, gelten beide Formeln auch noch fiir A = 0. aos 

FaBt man jetzt alles iibersichtlich zusammen, so erkennt man, daB sich die fiir I 
bis VII aufgestellten Bedingungen untereinander ausschlieBen, also eine eindeutige 
Zuordnung zwischen den Koeffizienten der Gleichung, den gegebenen Anfangswerten 
und der diesen entsprechenden speziellen Lésung herstellen. Die nachfolgende Tabelle 7 
gibt diese Zusammenfassung wieder. Dabei wurden die Schranken fiir a2 durch die 
Schranken fiir die Anfangswerte ersetzt. 


5 0>»D 
= 4-—VWZ a&—VWG ge &4—VWG ad — % a A — a Z zt = 
( MG “ MG Gy < cob ie | we+e M e+e | Tees were | Se MG | MSt 0 
‘ a0 : 
a a ies | ad i | a VY G 4— J G | a ZL 0 | Wf z4d a — o< 
. | “= | | 7 : 0 
o> 0< OS sei. 0s Of. ei a a= G2 she | 0% | OF ah 90S o> | o< eis a 
wn) | awn | Ce Gn) (en) | (q (g*2) | (e (gL) | (q (zt) | Wegwersy\s | 
a}OMSSULTUY Joep UseyolozIoA WOA JlOysisavyqy ud uoA sunMUIYSOg 
aah 4 | Coy oO ee I T= 
ae e zs D— |" 80D > > Ph > (Pkg +7)°4— ‘o> PFRI0 > 0/00 
a D V-a| a% |\y—alaesty—| 7° if eAgst+vV) * ore , 
o 
= oh + 2? | g.|\¥—-d| az_| aw—v aioe ? pb Saas 
£ —S sd Vena eo p+y—|\n os 0 e > Ph > (PAT+V) LA — 0< RMA IO >\0 >](4 
& [°A | F-@| az ivt+a|qwesy |” ve Avtget+y) * : 
= 4 ; 
o Oy ee m4) ai r =. : 
= ([)- A dé _|dz% FT PErV vg-y-|n ou vip—< (hg +) Ph — > 2h S00 + hlO>/0 > [ie Fz) 
Sp z\ » I d r—div+al| geur Zz V V we aie Vo)kc Z 
I ate ae 
Se Of wae ! es 
S| i dé_ | dé Si lad .VZ-V—|n ou vo <a LG Fr) gh — > 2h DT = 4K 0 >|0 <](o 
3 te) u/| UD aly ol eee L.Po-V Ver (°Ag+PV)s a) eG ( 
Z ; - z X 
= [4 | CLV tV i eo=2| | Pes av = oR (.0 5 apse PIR Sek ig 
v— |n sun >i AI VY) Ch = eh I 0. = 
i oo Sl peta iy Pola ee Pies er eee a eae : y= (q 
q 
= ee = Bs 
DB D é GHAV AS hel Ne 0 n us D > pis Re los Cs KR sy Lh Sh S019 S10= betes 
4 fe|. |\ytala—-vViv-a|. re— | 2° Hae Kee eh ge +). V ; : wet 
oO 
3 ( ee G Vee GA ae Die? , ne ae of (yaar of 2s op 
oe) po Sh eh gd + Vg he eh SO > eh > 010 <|0>|(q 
g io ‘nT weaias |\y~dal geser_|. 7" e : i Van 
» Dp 
=>)]—T1/ dé Le 1 PETV Zz nw vip—< eal °f — <, °K “BIqotoq “A10 <|0 >|(e(z° 
ae} ; | Oe te Giy=a anh ;LVETV V (49 + Ve cy [qoroq ‘: 0 >} (ZL) 
wy | de dG TMEV Z of of = < 
aes Ly TMEV |" wWdIV +T< Siqoteq ,°A 0<|0=| (T‘L) 
dee) = ty ¢ \y-a\y—a|l a# ee 
Ree ee a a ee ne i 
d uis He Veo OH 4 | o x)h=h 
UIs ee a 3¥ al | z (2) ae arora? oie ee 
sunso7y ojerzedg Sunso'y 9UIOULOST[V 
a A Sen en Se i en San 
At+ro=n (db) g= 4 ‘buys = dus ‘y= ous 1 >t >0 ‘hgh + clehg 6 + ¥) = I ad Gea Oe | 


fh = (0) f A = (0) 4h 0 = figs +hp +h woe bunsoT apjarzads pun aurawabyp ° 4 9112q0L 


ad 
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Hauptfall 2: D=0, A +0 oder D=|A\ +0 
Hier lautet die Identitat (2,1) so 
A\ 
Yor= —B (0? + aa oder (2,12) 


A 
Yo — By? (ye ae +): (2,13) 
Als Integrale treten, da B+ 0 vorausgesetzt ist, nur zusammengesetzte Kreis- und 
Hyperbelfunktionen, also ausgeartete Jacobische elliptische Funktionen auf. Diese 
entarten bekanntlich fiir 2 = 0 in Kreisfunktionen und fiir k? = 1 in Hyperbel- 


funktionen. 
Speziell ist fiir 


eS okies (2,14) 
snu = sin u, cnu = cosu, dnu = | 
bzw. fiir ior", LK “egos (2,15) 


snu = Tgu,cnu =dnu = Sec u. 
Die Bestimmungsstiicke der ausgearteten Funktionen I bis VII, jeweils unter AuBer- 
achtlassung des Falles B = 0, ergeben sich einfach durch Spezialisierung der bisherigen 
Ergebnisse fiir k? = 0 baw. k? = 1. Fiir diese festen Werte ergeben sich aus Tabelle 5 
fiir a? keine Schranken, sondern ebenfalls feste Werte, desgleichen fiir wm? aus (2,7); 
A aus Tabelle 6 nimmt, abgesehen von 4 = + oo, was B = 0 entspricht, die Werte 
0, 1, —1 an, wodurch sich im Verein mit (2,12) und (2,13) den ausgearteten Funktionen 
I bis VII die speziellen Werte D = 0 oder D = | A| eindeutig zuordnen lassen. FaBt 
man alles zusammen, so erhalt man Tabelle 8. 
Tabelle 8 


Vi (oe) 


TIIA1 <1, Os 419 Ze 
cA} 
B 


A 


LV 
4 Ap ee 


Ai) a ey es 


= 

| 

bo 

| 

| 

by) be 
a | | 
bal | yl | So) BS | ba] 

(=) 

——— 


VII} (—2, 1, —1, 0) - A 1; 1 |}|—1] 0 0 —A 


2B 2B 


Nun ordnet man nach D = 0 bzw. D =| A |, bestimmt die Vorzeichen von A und B 
nach Tabelle 1 und berechnet sin y nach Tabelle 7. Die dort stehenden Formeln behalten 
fiir k? = 0 bzw. k? = 1 ihre Giiltigkeit, da ja sn? u + cn? u = 1 auch fiir smu = sin w, 
cnu=cosu bzw. snu=Tgu, cnu = Secw richtig bleibt. Die entsprechenden 
Schranken fiir y,” und die Berechnungsvorschrift fiir y ergeben sich unmittelbar aus 
den elementaren Eigenschaften der entsprechenden Funktionen. 
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Hauptfall 3: D= A =0 
Hier hat die Grundgleichung (1,7) die Form 
y’ +2By=0. 


Also ist 
gira! Bys 
und 
bri ‘ i! 
sea x aa Se 
MS Ong Oy! 4 = oy 
Yo 


Die nachfolgende Tabelle 9 enthalt die Ergebnisse des Hauptfalles 2 und 3, wobei die 
Bedingungen fiir y,’ fehlen, da dieses durch A, B und y, bereits bestimmt ist, ferner die 
expliziten Transformationsformeln (1,11) bis (1,14). 


Allgemeine und spezielle Lésung 
von y” + Ay + 2 By® = 0, y (0) = yo, y’ (0) = yo’ 
ae Or) Ae OOden y= Al | TP OUk 0 ks 1 
Ketan) yon —— Ovi res inane erp) 


[EOS ep prsing = Tey. y= F(S., “| = In (tg + sec g). 
RB Anfanus- D Allgemeine Losung Spezielle Losung 
Were k? Ol == 0) (G3) a w? sin y Yo >9| YQ <9 
Pe i A A |Yo!| | 
<0 05 yo" < —oR 0/1 jasnu=aTgu ery hea - y>0ly<0 
A A A a 
= aCt a 7 > Ol pa0 
<0 Yo ae FF 0 1 Jans u=aCtgu op! 2 oF y> y< 
EA : —A | 
Sa == — >0 0 
>0 Yo =ZR Yo Oo} 0 | 0 jadnu=a ap s cea a a< 
y : i A \— Al: |Yol ax 
<0 Yo 0 | 0 jaseu=ateu 3B |? lVar+y. y |\n-—y 
A : A : a Mae 
2 > ——| A | 0 | = a cosec wu | ——| A —.- ytr 
<0 Yo = —zR ansu=a =p eal y 9 
etd. acnu=adnu=|—A Ya VEN ox 
>0 0<y¥"°S yz |-4 ] =a Secu Re = 4) te) a 0|a<0 
te z i ; 
<0 Yo + 9 —A| 1 jacs wu =aCosec u BR ja tye y>0\yv<90 
3 DA) 
Baz0 + 0 ! Y ! 
< U; ’ = ——————— > aa 
Yo y aV—Bt+y Yo 
, (A +2B y,°)? 
Transformation fir D? = — 42, D=idA, B<0, yo? > ee 
= ee ae Te eas 
,_ VAP + A? 14a eas a 
Sea TS 1— wu’ ° 2 By? —VA?+ A? 
age AB SAR As 
; / y zs 2 - 
ieee (+ ay al 


Mo PAU + 2(A 4+ VA2 + A?) u3 = 0, D’ = 2VA2 + A?. 
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3. Anwendungsbeispiele 


Praktisch geht man zur Lésung des Anfangswertproblems y” + Ay + 2 By’ = 0, 
Yo. Yo so vor, daB man zuerst D? bestimmt. Ist D? < 0, so fiihrt man die in Ta- 
belle 9 angegebene Transformation durch. Ist D? = 0, so tritt genau einer der drei 
Hauptfalle ein, wofiir die Tabellen 7 und 9 die den gegebenen 4, B, yo, yo’ entsprechen- 
den Lésungen liefern. Einige Beispiele aus der Praxis mégen dies erlautern. Dabei 
ergeben sich bereits bekannte Resultate, doch ist im Einzelfall der Rechenaufwand 
weitaus geringer, da sich die Lésung einer algebraischen Gleichung und die nachfolgende 
Transformation des elliptischen Integrals erster Gattung auf Normalform eriibrigt: 
andererseits gelingt es mit den entwickelten Formeln, viele ahnlich gelagerte Probleme 
einheitlich und allgemein zu behandeln, wie das erste Beispiel zeigen soll. 

1. In } 3, 28, 14, 20, 13 tritt ein Anfangswertproblem der folgenden Bauart auf: 


y" + 2b y + 2b ey? = 0, y (0) = Yo ¥ (0) = Yo. & S 9. 
In den oben genannten Arbeiten wird fiir spezielle 6? und « bzw. fiir yy) + 0, yo’ = 0 
oder yy = 0, yy + 0 die Gleichung untersucht. Hier soll die Gleichung unabhangig 
von speziellen Annahmen rein mathematisch gelést werden. Die spezielle physikali- 
sche Deutung wird nicht durchgefiihrt. Es ist fiir y” + 26%? y+ 26? ¢y? =0 und 
beliebigen y, und y,’ zunachst A = 26?, B= b?«, D=4b*[(1 4+ «€ y,?)? + « B? y,?] 
mit £?b?=1 oder wenn man d = ,|/(1 + € yo2)? + & B? yo? setzt, D = 267d. 


Nun sind folgende Falle zu unterscheiden: 
le> 0, Ara? O40, Be aa DIS ae 
Wie=—4<40, A=]20s 0, B= — bs 20, B= 2h aie 


d= ya £1 Yo) ep pe 


(te 22 ; 
eek Rn dann ist d reell 
€,B 
1 A 
phen 
a) Yo eee ae! 
1 
b)} je > = 
& 
eon A ¥ A ) 
OE ee Oo ee 
i Vom hed 
— 2 
Bp) ihn te ee SRLS, 


; : 2)\2 
B) Yo 2 : ae 
a), b), «), B) wie oben und 


i A 
2s = 
) == —R 


\ 12 (1 — & Yo")? * 6 = 
C) a6 > — fe > also d=16, D=2b261, A=2676 
mit 6= Ve, B? yy —(1—2, tg eel 


Entsprechend diesen Fallunterscheidungen ergeben sich aus den Tabellen 7 und 9 
folgende Lésungen: 


I ily y=|/ cn (bx /2d+y) 


Pe Sot oh hee 
= — sing = ]/1 Pee 


A) Yo”? < 


, also d= D=0 


(a 
& 
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IT A)a) (7,3a) y =) —4 eel eave Se) 


kh? = Lad sing = yo Vets 
ApR)a) (3b) = yes ns (ba Vd +1 i 

mise sme = [A 
A)b)p) (7,80)y=|/*" ne (bx |2d-+y) 

aie — , sin g = yi se 


B) a) (9,1) y=j-Te be +>), sing =| ¥5| Ve, 


1: 
B)b)«) 6) (9,2) y= =Ctg (bx+y), sing =— Lica 
Ve, | yo! Vey 
B) c) Tabelle 9. Hier ergibt sich im mathematischen Sinne keine endliche Lésung, da 


namlich sing = 1 und daher y > co. Physikalisch dagegen laBt sich dies entsprechend 
deuten. 


C) Transformation aus Tabelle 9 


‘1+ 8 oe ——= 
= 8 (oa fi + B+ Led (oxi +845 
14+J1+e8 seo e,—V1 + & 
| a ————= ; sing = |/1 — u,2; uw) = ——= 
211+ 8 p= YI —ug?s Ho é, +11 + & 


Setzt man speziell 2 6? = a (Yo 0, ¥o = 0,80 ist fire > 0 d=1-+ ey,? und nach I 


/ ¢ 
y=yoonx|/ S(+eye), are 
bzw. fir ¢= — ¢,, €,: yo>?< 1, d= 1 — & y,’,, also nach II A) a) 
Snes . | &y PrN sr WN Sag 
y = yosn(v | aia ) + K aha el Seemann 


Das sind die Resultate aus **. 

Als weiteres Beispiel sei die Lésung der Pendelgleichung g + x? sing = 0 mit 
den Anfangswerten g (0) = 0, q (Q) = 0 (Q = Hochstausschlag) gesucht. Einmalige 
Integration ergibt 

q? = 2x? (cos q — cos Q). 
Setzt man jetzt }4 


sin = y; sin 5 = 7, late I, 
so folgt 
yl +P (L+ ay — 2x8 y= 0, yo = 0, yo! = 41. 
Mit 2b? = x?(1+ 7’), 2e=— Wp 5 ist d=, also nach II A) a) 
Y = HU, ke? = 4? 
und daher 


soa ge fa ) leh AA yin 
q = 2 are sin (sin sn xa), k® = sin mii 


Das ist das Resultat aus ™. 
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2. Bestimmung der Umkehrfunktion gewisser elliptischer Integrale. 


a 


i dt é 
—— —,.a> 8, OS y<a; y=y(ax) ist gesucht. 
a) 2 iWeeay a. B y yy g 
J ee ee 
Es ist ys =a, y = — Vy t+ 2P yr +at, Yo = — \/2 a? (a2 + B?), also 
y’ —2RPy—2y¥=0,%=4,Y) = — 202 (a? + B?), a> Bp. 


Es ist D? = 4 (B4 — «4) < 0, daher Transformation Tabelle 9. 


At = 4 (a4 — BA), y? = a2 T—™, thy = 0, Ug’? = 2 (a? + BY), 


=p? 

Al = 4°82, BUS 2a" 1B!) MO Aes fiat) 
a? —. B 
2? 
u=—cn(2aa+ K), also 


giten (2a%-+ K) 
CO oR (Qaa+ K)° 


Opie ove Oa oF, ARID OS hey 


2 tne 
tres ET 
J 


b) x ;a>B>y>0; y=y(z) ist gesucht. Substitution y = * : 


ap dt = 
ya — a2 t?) (1 — B? #?) 


» ho = 0, w= Y(1 — a2 u®) (1 — BP ut), wo’ = 1, 
wu’ + (a? + B?) wu — 2 a? 2 u® = 0, D? = (a? — B?)?, (7,3a) 
1 fe 


i 
G= a kh? = a: (Oo? ==07 ay ==); U =| Sn an, also 


Y= 418 OX. 
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(Hingegangen am 16. Januar 1961) 


Umstrémung einer Kugel in einem Zylinder* 
Von Mouaffak Daaboul, Damaskus 
Mit 9 Textabbildungen 


Zusammentassung. Das Problem ist gelost, wenn das Geschwindigkeitspotential in einem belie- 
bigen Punkt des Gebietes gefunden ist. Mit Hilfe der Greenschen Formel wird zunachst eine Integral- 
gleichung fiir die Werte des Potentials auf dem Zylinder aufgestellt, die durch eine Laplace-Trans- 
formation auf ein unendliches Gleichungssystem fiihrt. Die Werte des Potentials auf der Kugel 
und anschlieSend in einem beliebigen Aufpunkt werden durch weitere Anwendung der Greenschen 
Formel gefunden. 


I. Einleitung 


Das Ziel der vorliegenden Arbeit ist es, die Stromung einer inkompressiblen, idealen 
Fliissigkeit um eine Kugel in einem Zylinder zu finden, wobei der Mittelpunkt der Kugel 
auf der Zylinderachse, die ich als x-Achse nehme, liegt. Ich nehme zunachst an, die 
Kugel bewege sich mit einer konstanten Geschwindigkeit U in der Richtung der posi- 
tiven xz-Achse in einer sich in Ruhe befindlichen Fliissigkeit. Da aber die Stromung 
anfangs wirbelfrei ist, bleibt sie es immer und die Strémung ist eine Potentialstromung. 
Die Randbedingungen sind: 


oP = — Ucos®@ auf der Kugel 
2 — 0 auf dem Zylinder und im Unendlichen, 


wobei y das Geschwindigkeitspotential, n die Normalenrichtung im betreffenden Punkt 
und der Winkel zwischen der «-Achse und dem Fahrstrahl des Punktes ist, wenn wir 
den Koordinatenursprung in den Mittelpunkt der Kugel legen. 

Den Fall einer ruhenden Kugel in einer mit der Geschwindigkeit U str6menden 
Fliissigkeit bekommt man, wenn man der Fliissigkeit und der Kugel die Geschwindig- 


keit — U erteilt. 


* Auszug aus der von der Fakultat fiir Naturwissenschaften an der Technischen Hochschule 
Wien genehmigten Dissertation. Referenten: Prof. Dr. G. Heinrich und Prof. Dr. H. Parkus. 


Ingenieur-Archiv XVI/2 10 


146 M. Daaboul: 


Il. Harmonische Funktionen mit einer auf der Kugeloberflache 
verschwindenden Normalableitung 


Ich suche nach einer harmonischen Funktion mit einer auf der Kugelobertflache 
verschwindenden Normalableitung. Fiir diesen Zweck sei H, in der Abbildung } 
ein Punkt auBerhalb der Kugel mit OH, =c, und H, sein durch Spiegelung an der 
Kugel entstehendes Bild mit OH, =c., d. h. c,-c, =a? (a ist der Kugelradius). 
In A, befindet sich eine Doppelquelle in der Richtung H,O von der Starke uw, = 4% 
und in H, das Bild dieser Doppelquelle, d. h. eine 
Doppelquelle in der Richtung OH, mit der Starke 


a as 
Be bee = ae 


Das Geschwindigkeitspotential in einem Punkt , 

x, welches von einer in H gelegenen Doppelquelle von 
der Starke « und der Achsenrichtung (J, m,n) her- 
riihrt, ist 


ey oe 
Abb. il Pp ray 4n y2 ) 


worin # den Winkel bezeichnet, den die von H nach M gezogene Gerade 7 mit der 
Achse (/, m, n) einschlieBt. Daher ist das Geschwindigkeitspotential, das von den beiden 
Doppelquellen in H, und #H, herriihrt, gleich 


fy cost, Me COS By cos B, a cosd, 
TA Tea dine, 742 4H SP on ne : CATE ee 
Oi! CO ae Cent GON Gy 1 
Fl Ota tae! Oe eamt a OCs, Fic. ee a BO a Ma 


2 
a- 


esi de, 
Cy 


und integrieren beziiglich c, von c, bis oo, d. h. beziiglich c, von c, bis 0. Definieren wir 


Wir betrachten c, als Veranderliche, multiplizieren beide Seiten mit dc, = — 


co 


weiters @. = | gy, dc,, dann erhalten wir: 


1 
=a +\a awa & (1) 
Im letzten Integral haben wir wu statt c, geschrieben 


he a 2S aie, oe 
r,*=— OH, +OM — 20M -OH, cos x MOF, 
=r+ u?— 2rucosa 
(OF = 00M =r] = OR eee 


2 hat ebenso wie p, eine auf der Kugeloberfliche verschwindende Normalableitung, 
weil c, ing, nur als Parameter auftritt. Wenn wir H, auf der Zylinderoberflache nehmen, 
dann ist y, in jedem Punkt des Gebietes, das wie in der Abbildung 2 durch 
S,, S, und S, begrenzt ist, stetig und hat auBerdem stetige Ableitungen beliebiger 
Ordnung. SchlieBen wir das Gebiet durch zwei Kugelschalen um O mit dem Radius 
d co, dann kénnen wir auf das entstehende Gebiet den Greenschen Satz anwenden 


ae 2a 
| | vege as =i P n, 4S. (2) 


Ss 
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Hier ist m das gesuchte Potential, und fiir gy, ist H, ein fester und M ein veranderlicher 
Punkt. Die Glieder in (2), die aus den Kugelschalen entstehen, verschwinden, weil 


y (bis auf eine additive Konstante), of , Y, und (solange H, eine endliche Ent- 
fernung von 0 hat) gegen Null streben, wenn 
M ins Unendliche gebt. S Ky 


AuBerdem verschwindet das Integral auf e 
der linken Seite fiir S, und Ss und das auf 


53 
der rechten Seite fiir S, und es besitzt weiters 
fir S; den Wert — 27 qm, wenn 7, gegen 
Null geht. m 
(2) reduziert sich dann auf ite 

7 ay i Be 
foo: on O51 = — 22 Gm + | \o-tds, 
‘st ss 
oder 

aK ae ae ase, 

Su s, 


Ill. Integralgleichung fiir den Potentialwert auf dem Zylinder 
Wir berechnen nun den Wert von 
aah F) 
| lo: : an ds,. 
me 
Wir verwendeten die Kugelkoordinaten (0 x Polarachse) H, (c,, , x) N (a, 3, y) 
(vgl. Abb. 3) 


4 
<¥ = —U cosd 
on 
c. 
1 ay aoe 
ae rh a ar ar ry OU 
0 
idle ep 
1) Ce an ds, 
1 1 Ms 
(a ce? +a? —2¢,acosa 
Lite i re 
meee (=) = P,, (cos «) pth 
Siem i 


P,, sind Kugelfunktionen. Nach dem Additionstheorem fiir Kugelfunktionen hat man 
P,, (cos «) = P, (cos w) P, (cos #) + 2 (ox A z “i P™ (cos w) P™ (cos #) cos m (y — yp). (4) 
m=1 = 


P” (cos w) sind die zugeordneten Kugelfunktionen. 


] op 
ele ont 


0 


eed 


" (— U cos? 
| ‘7 
Q) 


a? sin 0d0 dy = 


Ua f | yi/e\ gt e: 
ee | Gas P,, (cosa) cos}: sin ddd dy 
0 a 
= —"% and (=)' P, (cos w) | Px (cos 8) - cos # sin Odd. 
Cy cy 


eo é 


10* 
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Setzen wir cos # = ¢ = P, ein und verwenden wir die Orthogonalitatsbeziehungen der 


Kugelfunktionen 
te 0 fiir m + n 
{ P.(0)- Pa (6) do = nite 
2, Ini fer m=n, 
dann bekommen wir 
1 
1 } ap 22 U a? Nel a \r 
[Jz a, =-B2* S(2) pean rwre 
= =! 2 Uae (5) 
E YF Ba) oP me 
*) {| T On ds. 
S; 
Hier ist ; : abe 
Uw \n 
=_ : 2 P, 
r= Va? + uw? — 2au cos « Beles) es) 


weil w< @ ist. 
Nach demselben Verfahren bekommt man 


Lee 4 4x 
| ee py ener pe eae ee tp Cen 
Oe eR eee, Bi wiih <G) 3 

S, 


i oy Bet le ee. iT ee ie eral: 
J | Pala, Gs =] | ay ae mah an (fs du Ty in) ds, 
S, S 


AY] 
2 
4 a 1 O1arieeg 
ak le aap poe : 
gu U o,2 COS @ ! fu | | runes ax) du 


Cy 
4 a Le 4n 
caches ted Say § pen Lee jee eee 
= 3 7 U (2 C08 @ 4 all 3 u U cos w) du 
0 


at ee U7 as 4 2U coso c,? 
Sie ys Oe © 3 = 9 


—2nU @ cosw 
(Ore 


oy 22 Ua cos w 
{fo “an ds, = +E (6) 


2 
Gy 
Si 


Wir durften die Reihenfolge der Integration und Differentiation vertauschen, weil 
der Integrand gleichmafig stetig ist und gleichmafig stetige Ableitungen in bezug auf 
seine Argumente hat. 


Setzen wir (6) in (3) ein, dann bekommen wir 


22 Ua cosa 


¢ ris é 
2 Ae a aes + | Jo: an 02 Us: (7) 
Sy 


Diese Gleichung ist eine Integralgleichung mit dem Kern = Y2, der nur bei Zu- 


sammenfallen der beiden Flachenpunkte H, und M singular ist, und zwar, weil die 
Begrenzungsflache tberall eine Tangentialebene besitzt, nur so schwach, da8 der 
dreifach iterierte Kern beschrankt bleibt, so da die Fredholmsche Theorie auf diese 
Integralgleichung anwendbar ist. 
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IV. Der Zusammenhang zwischen Fourierkoeffizienten des Potentials auf der 
Kugeloberflache und den Werten des Potentials auf dem Zylinder 


Fiir die weiteren Untersuchungen leiten wir nun eine Beziehung zwischen den Fou- 
rierkoeffizienten des Potentials auf der Kugeloberfliche beziiglich des Systems der 
Kugelfunktionen und den Werten des Potentials auf dem Zylinder her. 

Zu diesem Zweck wenden wir noch einmal die Greensche Formel an: 


J Jered] Jose 8 


Dieses Mal wahlen wir fiir y’ die tee sie sind durch die Beziehung 
Py =e | ape se (¢) 


definiert, wobei ¢ = cos? und P(¢) die Kugel- a 

flachenfunktionen y-ten Grades sind. Diese Kugel- 

funktionen sind harmonische Funktionen. x 
Das Gebiet bestehe aus dem Zylinder (mit dem 

Radius 6), der Kugel und aus zwei Kugelschalen ‘Sy 

(S, und 8,) um O mit dem Radius d= co (Abb. 4). 

Die Glieder in (8), die aus den Kugelschalen ent- Abb. 4 

stehen, streben gegen Null, wenn d- oo, weil der 


> 


Flacheninhalt S, und S, endlich ist und weil q,’, se und Se gegen Null streben, 


n 
wenn d-— co geht. 


Die Formel (8) reduziert sich also auf 


ts aay’ apy’ 
JJws + ds, = | fe: fas +] [e- - ds». (9) 
Si Se 
Ich rechne nun die Integrale in dieser Formel aus. 
ek 
| fer ae as 

Ss: 

OU seams eS ere - Ucos8 = — UC 

ds, = a? sin 3 dy dd = — a*® dy df. 

ey eel 
Py (é 
[for as, =f [2 (00) (etapa = 
“Si Ores 
u 
—22U( —4nxU 4 
ne VP. (O)tde = Ser Ot = — aU 5. 
= 4 
ite oe 
yO eee: 

s 
Dryas. apy == (ae 1) 

p=¢ (9) =97 (2), ww =r on gral an\s, av+2 P, (¢) 

Ub 
: Z 1 : 2n(v +1) 

[ [o- Bas, = + 8a |@ (0) Py (0) db = —— Ay, 
“s, ay 


wobei A, die Fourierkoeffizienten von y auf der Kugeloberflache in bezug auf die Kugel- 
funktion P (¢) sind (siehe °). 
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“0 ag’ 
3) | Jean dee 
Sy 
: py’ Py (¢) ae (oad) _ or | 
Pp 5, = YP (Cc) = YP (2); an Ss =< rv Ting . Wa & = yv +2 Ps (f) on / 5 
ae 2... : : ( ob a sin? & cos & 10 
ake ab) = b (10) 
er ON era eee b : 
Sa | ee ee 
: i a0 = cos? B® ~—s sin 8 cos B ee 
weil tg 0 =—- ‘9g Siitim ga Fe b , 1 =a ist. 
Daraus folgt 
op ¢ Py’ (0) + (vy +1) Pr (¢) Sorc _ Pv¢iis) star Wes BC 
rll a brs sin 0 = — 72 Sin d. 


Die letzte Umformung folgt nach der Rekursionsformel (siehe *) 


(voted) se (CEP Cul (Cy a eee a (Gye 
Also: 


2 

Sa Oy’ . Poo (6) : ‘ * 

| ly on ds, = fo: ER sin +t 0-bdxdy 
0 


= 20d? (pay, 


rv +3 


Setzen wir die gefundenen Werte in (9) ein, dann erhalten wir 


(oo 
0. 1 im coger 
et = — ee NA, b 2a bP (peeetOL 
oder 
aie terete (y +1) 2U 
Jo- pes 42 = Tag, Av — Spe byl. (11) 


V. Das Gleichungssystem fiir die Fourierkoeffizienten des Potentials 
auf der Kugeloberfliche 


Um die Integralgleichung (7) naéher zu untersuchen, fiihren wir die folgenden 
Umformungen aus (Abb. 5). 


i : : 1 : 
Wir machen wieder von der Entwicklung >” nach Kugelfunktionen Gebrauch 
2 
a 1 
TY Vu +r? —2Qurcosa 
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Setzt man in (7) ein, dann folgt 


@ i 
2% 9H = we coso + |g = Fi * ds, + 
S," 
a Ep @ Sy ett 
ih i ? en - eth) P,, (cos «) ds, 
—co 0 az 
_22Ue ie bP a ra 
= as ae Ag Sleek 
cos w + | ein (Jat ole oe 2 (ey vee (cos «) dx dy, 
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weil der Integrand stetig ist und stetige Ableitungen hat und weil die Normalableitung 


auf dem Zylinder von y unabhangig ist. 
Mit Hilfe des Additionstheorems (4) bekommt man 


9 22 Ue 
4 P Hi = — 7) 


cos w + \ fe ; — (-) ds, + 
Se 


sees 
226 Nev) (Owais @ Py (cos®) 
Yi She Gate. (cos w) paca eed 
Da aber 
é@ Py (cos #) Py +1 : b P’y + 1 (cos #) 
a are hy Ee sinv’t? 9= ea ws 
ees 
: eae Cy 
ist. hat man 
; 22%U a3 ; a7 
229m = ie COS @ +| \¢ oa (,) ds, + 
Sp 
onder & 
ub? 7 i Ges P’y +1 (cos 9) 
= a jo 2a ovat Py (cos Ort pore 
Uae it i 4) ff 
GH = cane Oe m+ 2a) | aaron (dss + 


se 7 vy azv+2 


eae 
C Sirradt Ae a evri Py (cos w) ees rv +3 dx (12) 
a Uae — ON ve 
pm = 3 ee cos w + es oy (=) A,* P, (cos w) + 
Gy seen 
tan) Jos on (;) dso, (13) 


und zwar nach (11). 


Ich bezeichne das in (13) auftretende Integral mit J und ferner die Zylinder- 
koordinaten von H, (c,, , x) mit &, 6, z, die von N (r, 3, y) mit x, b, y und die eines 


beliebigen Punktes M des Gebietes mit 2,, e, y,. Wir haben dann (bei festem H, und 
veranderlichem M) 


— £2 + 9? + b? — 206 cos (y — 1) 


‘ eg bee '(y =P) 
fle 1 : 
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Lassen wir M-—WN streben, dann haben wir 


3 2bsint4—* 7 

om _ _ __4_ y= (w — 88 + 46% sin? 45° 
oder 

é (1 h? 

a | = se spe t) 3K (ae 7), 


wenn wir 25 sin“ 5 mit h bezeichnen. 
Alle Glieder von (13) sind Funktionen von £; multipliziert man beide Seiten mit 


Aol 
\2% i 
integriert iiber € von — oo bis + oo und bezeichnet 
+ co 
1 : 
SS SS 7a 
we (a) = =| mn (6) FE, (14) 


dann ergibt sich (man darf die Reihe gliedweise integrieren, weil sie gleichmaBig 
konvergiert ist und aus stetigen Gliedern neuigt 


Sie 2) 
me) = oe |e ects + — ye avt1 A fe cos @) eesde + geld 


3V2n J % 20 —, evi 
. ars (14) 
wobei 
cos w = ae c2 = & 4. 6? 
cos } = SS vy? — g? + b? sind. 
x 
Wir rechnen die Integrale in (15) aus: 
ee 
1 —— {Teed 
) ral $ 


+ co +oo 2x 


aa Teta = aah ee ‘ it |* (x) K (x — &) try 


Schreiben wir € — x =u, €= 2+ u, dann haben wir 


2x [tet am ae [ear (x) ts) i{ fx (u) e'*" du |dy |, 
6 


== (c=) 


\ 


weil g unabhangig von yp ist. 
Beriicksichtigen wir (14), dann konnen wir die letzte Gleichung in der Form schreiben 


aree 2% +0o 

1 b aay 

ay {Fett ae = 35 Pa: | fea" K (wu) du dy. 
— Go 9 as 


In dem Doppelintegral der rechten Seite konnen wir entweder iiber u und dann iiber w 
integrieren oder erst tiber y und dann iiber w. 

Beriicksichtigen wir die folgende Integraldarstellung der Hankelschen Funktionen 
im Falle rein imaginarer Argumente (siehe 8) 
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-- CO + co 
eB x ? 


d : 
apa = [EPs 
(eS 


+ co 
ed i a dian | @) (| Bl) = 2/8! 
tc aye ln ‘a m|B|*H,® (|B|)=2|B|k, (8), (16) 
dann bekommen wir im ersten Fall, wenn wir erst iiber uw und dann iiber y integrieren 
+ oo 2% +00 
lepb sae: b ¢ h? eau 
ee rag SS Eee ME ER i Ss 
= Sea vs | IE but + raya, Oy 
—co 0 —co 


= G2 \2|ah|b, (lab) -d 


 \ab sin 25" ky (2 ab sin4 >) dy 


2 


uw 


= "28+ {asin y ky (2b asin y) dy. 
0 


In dem zweiten Fall, in dem wir erst iiber y und dann iiber w integrieren, bekommen 


wir, wenn wir die Transformation u = hv = 26 sin oF -v ausfiihren und die Inte- 


graldarstellung der Besselschen Funktion nullter Ordnung beriicksichtigen : 


[ensnvdy — = 2{ eb sinv dy ned 6 (B) 
9 0 


co 2a 


1 ie oe : 5 
one iat = 2iabv sin (p— x)/2 q do. 
nal EE =| (Vo? + 1)8 | ‘ a ee 
0 


Jo (2xbv) 
+1) 


Das Integral (es dv hat einen Sinn und ist beschrankt mit der Schranke 1, weil 


Jo (2a bv) d r 1 
hc ae “|v 7) ia 
0 

ist; wir bezeichnen es mit k,, (b«), d. h. 


Jy (2 0b) wi ; " 

iF dv = — {asinyk, (2 basin y) dy == ky (b«) (17) 
3 

A (Vo? + 1) - ‘ 


und 
wale ee =e (b a). (18) 
"f Pa (cos co) 
n(COS@) 4 > 
2) ose eg 8 dé = ©. 


Durch partielle Integration findet man 


ry dey (cos ) 
a Ore ee 
io fe dE or dé. 
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Es ist aber 
d P,, (cosa) b? Py’ (cos @) (n + 1) Py (cos w)-& 
at “opti, ea cnt en +3 


wegen 


Wir machen -von den beiden folgenden Rekursionsformeln der Kugelfunktionen 
Gebrauch (siehe 4) 
P,! (cos w) — cos w P’,+1 (cos w) = — (n + 1) P, +1 (cos ) 19) 
P’,+1 (cos w) — cos w P,! (cos w) = (n + 1) P, (cos @). 


Wenn wir die zweite mit cos = — multiplizieren und zur ersten addieren, erh alten wir 
3 


b2 REF . 
oe P,’ (cos w) — wae ee oie. (n + 1) Pr+1 (cos @). 
1 1 
So findet man: 
@ Pn(cosw) _ (vn + 1) Pn +1 (Cos @) 
CL Cte ie ent? : 


Ks ist also 


+co 
a poner’ Py +1 (08 w) 


erin ents dé 
_ Gert 
Cn La Cn +1 
oder 
phn Sve (i anti 
On SS ap NAS eel lg ek ee ee Co; 
wobei 
Paar : eres « 
Cy d — — Yk 
litte” ipa ae oaer 
Ol, lie 
T Pa (COS@) 54 yA 
a= | dee eae ei) (20) 


Aus (15), (18) und (20) findet man 


AU aiathy (ba) | 2ky(ix) S&S (ia)vavt! 


iin an 


Pa = Px Ky (b %) 4 A, (21) 


Wenn wir A, durch seinen Wert in (11) ersetzen und beriicksichtigen, daB die Integral- 
gleichung 


+ co 
1 : 
Px (2) = =o (6), oe dE 
die Lésung & 


+ co 
y (8) = ial (a) e~ 95 da (22) 
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hat (siehe 7), dann nimmt (21) die Form an: 


220 a 
Px (2) [1 — key (b a)] = Ekg (b a) + 
2h (ba) ieee P’y +1(cos 8) 
toe a wm, =D + fae a 
_2t1U da 2hy bo) = B(ta)varv+! Were cos #) pet 
eee ae Se oeah |e [me ye-'*8de| da, 


und durch Vertauschung der Integralzeichen, was wegen der gleichmaBigen Stetigkeit 
und Differenzierbarkeit des Integranden erlaubt ist, erhalt man 


2iU aa 
+ (2) (1 — by (6 a)] = = key (b a) + (23) 
a y (ta)va2v+1 2 Venint- Py +1 (cos #) 
+z helba) 3 7 (v — 1)! ( i] me al w+3 dx dé, 


Wir konnen das in Klammer stehende Integral mit Hilfe des Ausdrucks fiir c, leicht 
auswerten; differenziert man namlich (20) nach 6, dann erhalt man: 


On Lag 6 Pn(cosw) 
ab =+ fe By apni oF 


Hier durften wir unter dem Integralzeichen differenzieren, weil die erforderliche 
Stetigkeit und Differenzierbarkeit des Integranden gegeben ist. 
Wegen (19) gilt 


@ Pn(cos@) — bé P,/ (cosa) b(n + 1) Pn (cos o) 
ab en t1 ae ent en +3 
@ P,(cos a) ; cos w Py’ (cos @) + (n+ 1) Pn (cosw) _ b P’n +1 (C08 @) 
ab Cel et en +3 aaa en +3 
Es ist also 
+ co 
8Cn ef 25 P’, +1 (Cosa) oe ee } 
a | ae dé = = (ia) a ke’ (ba). 


Durch Vertauschung von « mit — « gilt 
+ co 


"e—ta P’y +1 (cos ie 
i anes Cane 10 elcmiel 0 a) 
und nach bekannten Formeln (siehe °) 
keg’ (y) = — hy (y) 
ky (— y) = — ky (y) 


bekommt man 


ap Ss 
CES tag ] n+1 cos @ lé == 2 (— 1) . n I 
as b n! ( ) : 


ese ED, 


156 M. Daaboul: 


Setzen wir in (23) ein, erhalten wir 


@) 1 —& bal = bya) 


+ co 


L 2 SoM Macs seas = cA (ae es) 
cial The ie (b aes el (»— 1)! Jo (&) vy! ~ (y+) Ek, (bé) dé 


=o) 


oder, wenn wir 1 — k, (b «) mit J (b «) bezeichnen: 


27U @ a hig (b & 9 (D a) fou | Waert} (a Sy Ek, (b "ae (25) 


Pi) et (ba) rab a @— Dl o+1! 


vy=1 


Multiplizieren wir beide Seiten dieser Integralgleichung in ¢y, mit 


ai+1 bai (i)J 


ia aD ees 
und integrieren iiber « von — co bis + co, dann bekommen wir, wenn man 
oe 
2 (i) bai ; Pee eo Die on 
Sega pit We fee (a) aft hy (0d) do = dj (26) 
schreibt, 
+ co 
2bU (—4)i C., hy (ab) ky (ba) 
= = ‘ G VE ati lait? 1 (ba) da 
+ co 


2 ms lacy GY TEE DIB ass met ky (b x) ky (b a) 
| 2 Gait Go ee 


Fiihren wir die folgenden Konstanten ein: 


+ co 
26: Ree ie eb OO) oe (OC) 
Mee iy J Teg oe 
oder wenn wir 7 = mw schreiben 
=O 
2(—1)) j+1h, (w) ky (u : 
Bs Be CP wf on 0 (u) IA He es 


— co 


dann nimmt die obige Gleichung mit B; = = die Form an 


co 
Mjti , Ow My +j 


G+pit —(— " Gepre—m 
Mit Hilfe der M@; kann man die B; berechnen. 
Nun kénnen wir leicht eine Darstellung fiir ¢,, («) finden, wenn wir (26) in (25) einsetzen 
21 U a ak, (b a) 2) ky (Oo) ee, (Nae ave 
— }. —— s 7 ‘ 
janl@o) | |an le) = @-t or (28) 


v= 


pee. Bo =e (27) 


Pu (%) = 


Ich versuche nun die Formel (11) umzuformen und einen Zusammenhang zwischen 
A, und B, zu finden. Es ist nach (22) und & 4) 


= 
Py td v+1 
[« (€), 5 df == aa 2x alm 6 nfe a5 at — dé da 


co 


2 Ee Gs 
Von: RR 


ee (a) te Whe) da. 


LS) 


Umstrémung einer Kugel in einem Zylinder 157 


Setzen wir in (11) ein, dann finden wir 


+ co 
es CNC el ae ey 2U 
Jan dv! fer pu () hy (bada= iit Ay — 353 Oy, 
d. h. wenn wir (26) beriicksichtigen, 
A, 0,=—U «a B,, wenn vi | 
Ay = 6 = UB, (29) 


Das bedeutet, dali die Fourierkoeffizienten nach Legendreschen Polynomen des 
Potentials auf der Kugeloberfliche uns jetzt mit Hilfe (27) bekannt sind. 

Da dieses Potential auf der Kugeloberflache stetig ist und nur von # abhangt, ist 
es durch seine Fourierkoeffizienten eindeutig bestimmt und ist gegeben durch die 
Entwicklung 
2” = 1 


Ge) = A, P,, (cos 9), (30) 


So 


wenn diese konvergiert (siehe °). Das ist aber immer nach dem folgenden bekannten 
Satz der Fall (siehe 4): Die Entwicklung (30) der Funktion g (#) nach Legendreschen 
Polynomen konvergiert in jedem Punkt des Intervalls — 1 <cos# <1, in dessen 
Umgebung die Funktion beschrankt ist, gegen den Funktionswert, und zwar gleich- 
mafig in jedem abgeschlossenen Intervall, in dem die Funktion stetig ist und im Inter- 
vall — 1 <cos# < + 1 liegt. 


VI. Das Potential in dem Gebiet 6<R 


Der Wert von ¢ in einem beliebigen Punkt des Gebietes ist durch die Formel 


4% 9p = ih : ds, + | Je 5 (+ -)dsy + ae males )ds, (31) 
gegeben. 
(5) gibt uns den Wert des ersten Integrals 
ie : a ds, = : a U S cos *, (32) 
Sa 


wobei (R, %, y) die Kugelkoordinaten von P sind (Abb. 6). 
Wir berechnen jetzt das zweite Integral der rechten Seite. Zu diesem Zweck nehme 
ich zuerst an, die Kugel befinde sich allein im unbegrenzten Raum mit dem Randwert 


y = A, P, (cos-?). 


Das Potential hiefiir ist durch 
an+1 


pP = saa An Pr (cos ¥) 
oder mit Hilfe des Greenschen Satzes durch 
"1 ay Pete Cor Gen 
4nge = —| [+52 ds,+| [o- an(>)es (33) 
Sa 


gegeben. 
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gp und die beiden Integrale der rechten Seite sind wegen der Axialsymmetrie nur 
von R und # abhangig. Ich setze 


[ [2=t ds, = F(R, 8). 
“ge 
Es ist dann 


5 
eye er — i 
say) 2Rr 
und 
) ( 1 ) ( 1 R? — 7? — a 
on =) ea 7 2a7° 
Abb. 6 Also ist 
oF a? — R? CAgiin 1 
aR 2R \| a 6—zRP 
Daraus folgt 
Ain Pr ok BOE ] 
| | 2 4 = 3 _ pon | gm: 
Si 
ra) 1 k?—a Ag Tz 1 
UA es \] pi tae 
i S, 
R oF 1 
anon ae Pf. 


Setzt man in (33) ein und beachtet, daB die Normalableitung des Potentials auf der 


Kugeloberflache — A, P, (cos #) ist, dann erhalten wir 
an+t n+1 Via) 1 n R92 
4n A, Papa = a “4 eoRe a a oR 


Diese Gleichung ist eine Differentialgleichung, deren Lésung sehr leicht zu bekommen 

ist; sie lautet 

42 An Pn (cos #) an +2 
(l+2n) Re+1 


Nun ist K (cos #) = 0, weil, wenn R — ce strebt, F(R) gegen Null geht. 
Ks ist also 


F= 


LK (cos @) R’. 


F 4:% An Pp (cos 3) an +? 
a (2n +1) Re+t 


ae OED i R (/—4a(n+ 1) An Pr\ (a\n+? 
| 4, P, (cos w) alae rar a ( DP yeieek )(¥) 
Sy 


1 42 An Pnatt+2 
a (2n+1) Rr+1 


| \4, P,, (cos w) - = (--) as, = oy (3) A, P, (cos 9), 
Si 


und von dieser letzten Relation leitet man leicht ab 


1 


“i Saat g S da 
| | 2 5 A, P, (cos 0) (+) Cs 2a din A, (=) P,, (cos #) 
oe . n= } 
oder 
, a/l eS n+l 
| loa (=) ds, = 20 Yn Ad) P,, (cos 9). (34) 


S n=O 
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Die Annahme, die Kugel befinde sich allein im unbegrenzten Raum, hat uns zur Be- 


rechnung des Integrals | lo am (=) ds* gedient. 
Ke 


Ich komme nun zum letzten Integral in (31) 


ee 


2 


Wenn R <b ist, hat man 


fe al a/R 
ae om (=) P; (cos «) 


aN 3 Ge. ee 
YP (cos 2) - P; (cos w) + 2’ eS ai Pm 


ele nee 


1 (RV C GOs 
ealG VP ,/ (cos w) P; (cos 9) —S®. 4 
7 


Gj on 


€ cos @ 


i 
9 Y1j—m)! 
sae — Gem) [P;" (cos w)]’ Pj (cos 9) cos m (xy — yp) —z, 


Nach (10), (13), (19) und (11) ist (Abb. 7) 


dc b 6 COS @ bcos @ 
Dik Cr on * (i 
eG 3 Te 
‘ : y > P41 (Cos @ 
0p (=) den = 200d’ BP, (008 8) [2 dn 


n 
—co 


e . 
ae 5 ees 

2 7h ° 
e . 


Tr 


ly — (—) ds, = 2am + 1)A,° (=) P,, (cos #) a4 U R (cos 9) 2. 18) 


Aus (31), (32), (34) und (35) findet man den Wert des Geschwindigkeitspotentials in 
dem Gebiet R<b durch die Konstanten A, 


gp = 5 US, 008 oO Se 
SS fm fa\n+t je Uh ff 1 NE 
21/3 (z) a Ge |A P, (cos #). (36) 
Der Wert von @ auf der Zylinderoberflache ist durch 
l +co 
Y (x) = Von | P (a) Cn = de: (37) 


gegeben, wobei nach (28) 
2% Ua « ky (b x) 


v+l 


2 Ky (ba) yr ia)” a" 


= ——— we —d,, 
Px (a) }221 (ba) V2a 1(ba) > (v— 1)! 
(38) 
Ei, (— ta x)J 
é — = 
Faw : 


Die zweite Reihe konvergiert absolut und fiir jeden Wert von « und x. Den Beweis der 
absoluten Konvergenz der vorletzten Reihe bringe ich spater. Wir konnen daher die 
beiden Reihen multiplizieren und gliedweise integrieren 
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+ CO 
1 * Ny wi(—ta)i+1 hy (ba) 
ee 3 
OE) Vas | os j! 1 (ba) d 
tae (é.a)¥ (— ici + hy (dx) + av +! 
“yy yas (tay (— 44% 9 (ba) - a 
ae iba ue eS 
Mit den Konstanten 
+ co 
.. ( a) ky (b x) 
pitig |e da =D; 
schreibt man 
Ud y(—1)5 xvi Coe . ai dy Dj+vavrl 
PAB) senha s Di 37 aN n= | L) (jr eae ge Cy 
Jj ,v= 


Die erforderliche Konvergenz dieser Reihe beweisen wir gleichfalls spater. Ferner kann 
man eine Entwicklung fiir das Potential auf der w-Achse folgendermaBen ableiten 


4 JJe ps () dsy, 1° = b? + (% — ¢)?, 
} oe Su (=) = ga (Abb. 8) 
: co 

7 (foc) ee 


Abb. 8 So es 


Nach den Gleichungen (16) und (28) schreiben wir: 


+ co } + co + co 
] an = ¢ ae b2 eiax 
ge aa Le. de = 2x Jpla)re ware ent 


I, = 2-229, (a)bak, (ba) 


+ co 


Ue) afer" Ty, (a) da 
ae A 
I (x) = 24 Vad ( ei er aaule 
eae: 
: Pep d. fete art pe (40) 


Mit Hilfe von (31), (32) und (34) bekommt man den Wert von » auf der Zylinderachse. 


Vil. Untersuchung des Gleichungssystems 


Bevor wir das Gleichungssystem (27) untersuchen, wollen wir die Konstanten M ; 
abschatzen : 
5s 
‘ap wee 2) “uj+1 ky (uw) ky (w) 
Mie yr ey ia 


ah) 


du, 
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wobei 
2 : é 
1(u) =1—k, (uv) =1— = }wsin y ky (2 usin y) dy, 
6 
ist 
1) Mee j= 1, 805,05, 


weil der Integrand fiir solche 7 ungerade ist. 
Fiir gerades 7 hat man 


4 bare r 2j+1k k 
M2; == (— Ly weit U ay o (%) du. 
0 
eo ee Se ee 


2) Das Integral hat fiir 7 => 1 einen Sinn, weil der Integrand gegen Null geht, wenn 
uw gegen Unendlich oder gegen Null strebt. 
3) Das Integral ist fiir beschranktes 7 beschrankt mit der Schranke 


[Maj < oe +1 GN 


wobei c; eine gegen Null strebende Konstante bedeutet, wenn 7 gegen Unendlich geht. 
Durch numerische Integration findet man leicht die Werte 


M, = — 0,834 ue, M, = 2,480 ee, Me = — 10,944 yi M; = 86,400 we 
M, M M M 
(1!) = — 0,834 pv, a 4 = + 0,620 5, GB ByF ae 0,304 u, aye = 0,150 p? 


mit einem Fehler von 2%. AuBerdem nehmen | G! = | bei wachsendem j ab und haben daher 


die Schranke 
| M,; | 
ne 


Ich komme nun zum Gleichungssystem (27) 


<cp?it+! mit c = 0,84. (41) 


(—1)*-! Ma; ee see 
2 (8+ Wiel tt 
Beataber 0.0 fir 7 = 1) 3, .. und. = O'fir 7= 0, 2, ..., nach (29) wegen 
gy (x) = — y(— a) auf der Zylinderoberflache. Das Gleichungssystem reduziert sich 
damit auf 
—] Me (n+j +1) Mej+2 ae , 
> |a+ @j+2)! al Jepsen = — (2742)! pO 2 aa 
Mit dem Ansatz 
Bajar =p yr + BO (42) 
bekommt man 
=e F ME pee a) lehs(2) — Man+e = Maj+e 
= [2 a i+ 2)1(2 m)!| [Be eet CU ej) ein 0-3 a) 
oder 
7 M3) G+ 74-7) (2) 1 oa M2 (n, +1) M2(n +540) ‘ 
2'[>/ + gyrni ga) Be =ayta ay Om +21 2m)! (43) 
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Die Reihe auf der rechten Seite konvergiert absolut, weil nach (41) 


Ss Ma(n+1) Me (n+j+1) 7 al Me (n+1) Manjit 
Qj7+2HlQn[Ant2)!| ~— | (G+nt+ Ii? Qn +2)! 


| Mejia] 0 | Meman| © | Mac ue +1) | ~ (n+1)!(n+1)! ors pe 
<(G4 DIP =, @nt2! Site a,° @n+gr % 
| Miaj+2)| 48 | Mej+2| che 


<¢Tg epi l—@ ~ TG +i 1— a 


Wiederholt man dasselbe Verfahren 0 Mal (9 > oo), dann bekommt man fir Bete das 
Gleichungssystem 


4 | M ; 
ay} Miers Ee (atiby oo 
On! 6 ai | B 


tng 
(= Lesh 1 ae Me (m+j+1) Ma(m+m+3) +++ Malng_1+ng +1) Mang th) Lge 
— — Oo- +) 
ce kth?) bee Boeke (2m1)!(2nm1 +2)!... (2me)!(2me + 2)! 
Panes 
wenn man 
(0) (= })e ‘e ne peat Rast Se eae eee 44) 
Bein OF +2)1, eek “en m)!(2nt + 2)!... (2my-1)! (2me—1 + 2)! ata 
one 
setzt, mit der Abschatzung 
| | Mo; +2 | ( CHI if 
Sos1l<Te4 pi T= a 


Anderseits ist (siehe *) 


g | Mej+2| et e+o-1 eo, Ao+3Q+3 
otal SG Sy he 
l my t++. +g =G 


Ses 


Hier bemerkt man leicht, dab 
kleiner als die entsprechenden Glioder von |S, sind, d. h. die Reihe 8S, konvergiert 


gleichmaBig, solange uw < | ist. 
Wegen der Eindeutigkeit von 5; hat man nach (42) und (44) den folgenden Wert 


<1, wenn 9~ co, weil die Glieder von |S, .; 


Mo; +2 i y M2(m +1) M2(n Seine) 
B : cme HW » 1 1 
7+1=— Bet BH, Om +7Qny! 7 
4 (— 1je+1 y" M2 (ny +j+1) Pe alee M2 (ng +2y +1) M2 (ng +1) (45) 
(27 + 2)! mh ade: $aeat (ni)! (221+ 2)) oo. Onegin 1)! ; 


Diese Reihe konvergiert gleichmapig. 
Da M2; das Vorzeichen (— 1)/ hat, hat B2;,; auch das Vorzeichen (— 1)/. 
Ferner gilt Bo;+1—> 0, wenn 7 > oo. 
Die Entwicklung (36) stellt eine alternierende Reihe dar, die in dem Gehiet R < ) 
sehr rasch konvergiert. 
Summieren wir die Doppelreihe (45) spaltenweise, so bekommen wir 


1 | Maj+e| Ma M2 (j +2) 
B role — 5 J J 
| 2+ (27 +2)! iz |M,| Zs onanth DUAR =f 
ot ; '( pee |: ) 


al | Me Me (j +3) | | Ma Me Ma (; +2) | rT | Me Ma(j+a)| 4 
| 2 = EE Eee 
4161(1— 3) (2!)2 (41+ (1 ms Sly 6181(1— [Maly 
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, [Me Ms Mo(j+2)|+|Ma Ms Mo(j +3) | | Ma M Ma(j +2) 4 
Me \ \3 M2z|\4 
2! (41)26! (1 “ Se) ane aye (1 = oe) 
| Mio Ms (j+5) | ae ety Mss) +\M, Mio Ma(j+4)| 
Mal 2 [Mal 2 
8110! ( | 214i6!i(1— 21) 
4 Me Mie Maceo] [2M Me Magis | +/Met Maso 
Ma “ [Me]\4 
(4!)2 (6!)2 vee “ ny (21)2 (41)3 6! (1- af ) 
4 |M,? Ma Me (j +2) Beis 
21)5 (4 (1 Me y 
(21)5 (41) — 3H) 


Man sieht leicht, daB diese Entwicklung sehr schnell konvergiert, weil der Ausdruck 
in der ersten eckigen Klammer von der GréSenordnung y2/+3, in der zweiten von 
uw7~+10, in der dritten von w?/+14, in der vierten von m2/+18... ist. 


1 : : 
Im Falle « = 5 haben wir (durch numerische Rechnungen) 


B, = + 0,443 7s Bz; = — 0,028 we B; = 0,001 we Bz, = — 0,00005 Le 
Fiir den Fall u = . == es hat man 
B, = 0,560 w3 Bs = — 0,090 3 Bs = + 0,010 v3 Bz = — 0,0008 py’. 


mit einem Fehler von 3%. 
GemaB (29) hat man 
Ua 


A= 3 aw As OU 8 Bye 1. 


Setzen wir die gefundenen Werte in (30) und (36) ein, dann finden wir 


1) Fir den Fall u =5 auf der Kugeloberflache 
gy = U a [0,583 cos # — 0,012 P, (cos #) + 0,001 Ps (cos #)... | 
2) Fiir den Fall u = es 


y = U a [0,880 cos & — 0,14 P; (cos #) + 0,022 P; (cos #) — + ...] 
Im Gebiet R <b hat das Potential den Wert 


1] 
1) Fir (aaa 9 


gp = Ua [0,528 (e) cos } — 0,005 (5) P, (cos 0) + 

0) ,0005 (4). P; (cos 9) +).. | 
+. Ua ol (5) cos # — 0,056 (3) P, (cos 8) + 0,012 (F) P, (cos 3) + 

24 003 (5) P; (cos) +. . .| 

2) fir «= 3 
gp = va 0,63 ( a cos ) — 0,06 (-F) Ps (cos #) + 0,012 (+) P; (cos 0) — 
+. J+ (029 (; ) cos 0 0,18 (5, ) Ps (cos #) + 
+ 0,05(7')° Ps (cos #) -—- 0,01 (=) P; (cos 9) ) 


Hes 
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Der Wert von ¢ auf der x-Achse ist nach Gl. (36) und (40) fir 2 > 0 


pret a S 4p (Ep + 4S) + tae 


12) = 2100 | gar Moby hay 


T(ba) 
1 ee b 
a i le ay “J see Oe ae tt da 
oder 
+0co 
a : * sx hole) hy (x) 
1(=) =1(X) = 21 prale Bee cae 


+ co : na 


— co 


4 . 6 ° . M; 4 ; 
Entwickeln wir e~‘** nach ihrem Argument und bezeichnen wir rears mit M; 


(M; = 0; 9 = 1,3; 5; ..., dy = Ohy 2, 3,-6)2....), “dans bekommen wi, 
I (X) Ua j Az; +1 lw Xt pevte 
4n 4 A M3342 005 0) 2 pep (2)! M's +541) dey +1- 


Man sieht leicht, daB alle hier eintretenden Reihen konvergent sind und daf die Reihe 
iiber 7 alternierend ist. Ferner findet man leicht, daB (x) gegen oo we 2(1 ++ 2 B,) 
strebt, wenn x gegen Unendlich geht. 
Der Wert von p auf der Zylinderober fliche 

Die Gleichungen (37) und (38) geben uns den Wert von @ auf der Zylinderoberflache 


x iia Tie Mas ju +t XI Dj+v dy 
eG) == Wha te 2 1 
j ES, 
mit 
+00 ke, (B 5 ony (<x) 
al oye | &) Kg (& 
D; = b+ (i [ CE da = (i)' | WCW 


Wie bei M2; kann man leicht beweisen, da8 D,; einen Cae hat und auBerdem wegen 
Ke (@)o< 2—— (a > 1) die Schranke | D2;' <c, r|(2 9 + 5) +690 besitzt: fe; -cund se, 


sind beschrankte Konstanten). Die Reihe (38) ist konvergent, weil d2,,, von der 
GroBenordnung wt! ist. 

Die Doppelreihe (46) konvergiert fiir alle Werte von x. Ich suche nun nach dem 
Grenzwert von g(x) auf der Zylinderoberfliche, wenn x > 0 geht. 

Aus (37) und (38) folgt 
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co 


20 y . ky (bx) 2v+1 g2v+2 
neler 1) + (sin oo tes i= an dey +1 de 
x U wa ee ky («) 
»(X)=9(F)= 45 = (sina X yrs ada + 


6 % k a 2 2 
ein [sin XP) ge PY? day a do 


in Ubereinstimmung mit dem Grenzwert von p auf der Zylinderachse. 

Ferner bedeutet dieser positive Grenzwert, daB sich die Fliissigkeit auf dem Zylin- 
der in der negativen x-Richtung (d. h. in der Gegenrichtung der Kugelbewegung) 
bewegt. 

In einem beliebigen Punkt P ist mp durch (31) gegeben, wobei das erste Integral auf 
der rechten Seite durch (32), das zweite durch (34) gegeben ist; in dem letzten 
Integral von (31) brauchen wir nur den Wert von g auf dem Zylinder einzusetzen und 
iiber die Zylinderoberflache zu integrieren. 


Vill. Geschwindigkeit und kinetische Energie 


Die absolute Geschwindigkeit ist durch die Komponenten veg vy gegeben 


op eee lg 
CE eR v9 R 2d 
Die kinetische Energie ist 
y 1 ae Y 
P= — se] | o-o— 
(Ss) 
(4) yi 1 
Ks gilt auf S,: a = — Ucosd —» — A, P,, (cos #) 


é 
auf S, = 0. 


Daraus folgt 


= e 0 u| 1p 5 An Px (cos #) - cos 3 ds 


= — Sn oU A, a*{ cos? od cos d 


0 
] 
=xoU aA, =20U*a3(% + By). 


1 3 
9. Beisprel: (p= + = 5 und 7. 
Um einen Uberblick iiber den Einflu8 des Zylinders zu gewinnen, bringe ich die 
folgenden Tabellen. 
Die erste Tabelle gibt die Geschwindigkeit im Fall b > oo (d. h. die Kugel 
befindet sich allein im unbegrenzten Raum), und zwar in den Richtungen # = 0°, 


4 : 
30°, 60° und 90° mit dem Abstand k =a, 3 a, 24 vom Kugelmittelpunkt. 
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6 : ates 1 : : 3 a 3 
Die zweite Tabelle gilt fiir den Fall - =, die dritte fiir den Fall ==]. Abb. 9 
zeigt die Verteilung der Absolutgeschwindigkeit im Falle w= 5. 
Tabelle I. = 0 
bi 1 So 2 
a 3 
vo/U 0 0 0 
10 
vr/U| 1 0,422 | 0,125 
bs/U | 0,25 0,105 | 0,031 
0 = 30 
vr/U | 0,866 | 0,365 | 0,108 
bo/U | 0,433 | 0,182 | 0,054 
& = 60 
Dr/U | 0,5 0,211 | 0,062 
ps/U | 0,5 0,211 | 0,062 
o = 90 ? 
br/U! 0 0 0 
Abb. 9 
: 1 
Tabelle 2. 4 = — 
a 1,2 1,5 1,8 2 
a 
be/U 0 (i «are 0 
G0) 
br/U 0,580 | 0,298 | 0,175 | 0,135 
pe/U 0,267 0,157 0,078 0,032 0,020 
@ = 30 
Dr/U 0,867 0,490 0,241 0,143 0,105 
be/U 0,391 0,234 0,116 0,061 0,033 
o = 45 / 
Dr/U 0,706 0,390 0,179 0,087 0,045 
ve/U 0,507 0,304 0,172 0,108 0,080 
o = 60 
Dr/U 0,269 0,111 0,044 0,002 
vo/U 0,601 0,415 0,192 0,185 0,179 
v3) 
vr/U 0 0 0 0 0 
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Tabelle 3.. un = 4 


us 1 1D) uel 4 
a | 3 
va/U 0 0 | 0 
o=0 ih 
vr/U 1 0,66 | 0,49 
bo/U 0,402 0,226 0,136 
o = 45 5: 2 
vr/U 0,707 O.842 7 6.839 
ve/U 1,30 104 | 0,92 
d= 90 - | 
Dr/U 0 0 | 0 
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Der Tragheitspol der Kurbelschleife* 


Von K. Federhofer +, Graz 
Mit 3 Textabbildungen 


Zusammenfassung. Die Trigheitspolkurven fiir die schwingende und die rotierende Kurbel- 
schleife werden angegeben. 


Einleitung 


Bei der kinetostatischen Untersuchung ebener Getriebe erweist sich die Benutzung 
der erstmals von H. Winter! eingefiihrten Tragheitspole der Getriebeglieder als sehr 
vorteilhaft, denn sie sind die Mittelpunkte von geradlinig begrenzten Kraftbiischeln, 


* Die vorliegende kleine Arbeit fand sich im NachlafB Prof. Federhofers. Sie wird hier zusammen 
mit einem Nachtrag zu den in dieser Zeitschrift Bd. 1X, 8. 75, und Bd. XIV, S. 244, erschienenen 
Verzeichnissen der Arbeiten des Verstorbenen veroffentlicht. Die Schriftleitung 

1H. Winter, Sitz. Ber. Akad. Wiss. Wien, Abt. Ila 139, 151 (1930). 

R. Beyer, Maschinenbau. RM-Afg. 8, 642 (1935). 

R. Kraus, Maschinenbau. RM-Afg. 4, 337 (1936). 
G. Gerber, Maschinenbau. RM-Afg. 8, 533 (1940). 
O. Tolle, Ing.-Arch. 22, 227 (1954). 
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welche die Systeme der Tragheitskrafte der Getriebeglieder kennzeichnen. Fiir die 
graphische Ermittlung des Tragheitspoles einer zwangslaufig gefiihrten ebenen Scheibe 
gibt es eine Reihe von Konstruktionen!. Der Ort aller Lagen des Tragheitspoles eines 
Getriebegliedes fiir den Gang des Getriebes ist die Tragheitspolkurve*. Fiir die Koppel 
des wohl am haufigsten verwendeten zentrischen Schubkurbelgetriebes la8t sie sich 
mit groBer Naherung durch eine Gerade eho: zur Koppelachse darstellen*, die vom 


Schwerpunkt der Koppel die Entfernung -“ == “ besitzt (i, Tragheitshalbmesser, 7 Kur- 


belhalbmesser, o ='s/r mit s als Entfernung Ae i Salnwerpuniiee der Koppel vom Kurbel- 
zapten). 

Besonders einfache Ergebnisse bestehen beim Kreuzschieber ; die Tragheitspolkurve 
ist im bewegten System ein Kreis, im festen System eine Ellipse‘. 

Fiir die schwingende und rotierende Kurbelschleife liegen bisher noch keine Ergeb- 
nisse vor; sie sollen im nachstehenden entwickelt werden. 


Die schwingende Kurbelschleife 


Dieses in Abb. 1 dargestellte Maschinengetriebe besteht aus der um den festen 
Punkt 2 rotierenden Kurbel 2.4 = r und der im Kurbelzapfen A gelenkig angeschlos- 


Abb. 1 


> K. Federhofer, Osterr. Ing.-Arch. 5, 240 (1951). 
° O. Tolle, Ing.-Arch. 21, 365 (1953). 

4H. Weirich, Osterr. Ing.-Arch. 9, 230 (1955). 

O. Bottema, Osterr. Ing.-Arch. 13 103 (1959). 
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senen Schwinge A D, die in einer um C drehbaren Hiilse gleitet. Die Kurbelschleife ist 


schwingend, wenn C auBerhalb des von A beschriebenen Kreises liegt (AC =a> rye 


fiir a<r kann die Schwinge volle Umdrehungen ausfiihren, dann liegt eine rotierende 
Kurbelschleife vor. 


Der Tragheitspol 7 einer ebenen Scheibe ist durch die Lage des augenblicklichen 
Drehpoles P, des Wendepoles J, des Schwerpunktes S 
und durch den Tragheitshalbmesser 7, vollkommen be- 
stimmt. 

Man findet ihn zeichnerisch aus diesen Elementen 
durch folgende einfache Konstruktion (Abb. 2): Ziehe 
JS bis zum Schnitt W (Wendepunkt) mit dem iiber 
J P geschlagenen Wendekreis; dann ist 

ed 
Say 


Die Konstruktion verwertet die von H. Winter® be- 
wiesene Kigenschaft des Tragheitspoles einer zwangs- 
laufigen Scheibe, daB er identisch ist mit dem Schwin- 
gungsmittelpunkt der um den Wendepunkt W des Abb. 2 
Strahles SJ drehbar gedachten Scheibe. 


Bezeichnen rp, rj, rr die vom Schwerpunkt S gemessenen Ortsvektoren der Punkte 
aes = == , 2 
P, J, T, so ergibt sich, da SW = SP cos 6 ist, ST = —“* 


SP cos o 


und hiemit der Orts- 


vektor des Tragheitspoles ® 
as? 
"J tp-ty 


tr = — (1) 


oder mit D — PJ als Durchmesser des Wendekreises 


: 2 ts! (2) 
UE ogee ar ORE (ry? + rp?) * “ 


a) Drehpol und Wendepol 


Der Drehpol P liegt im Schnitt von A mit der in C errichteten Normalen zur 
Schwinge; seine Koordinaten in bezug auf das durch S gelegte (&, 7) Achsenkreuz 
sind weiter unten angegeben. 

Den Wendepol J findet man durch folgende Konstruktion’: Bringe die in A zur 
Schwinge errichtete Normale mit CU zum Schnitt in G, verbinde G mit P und tiber- 
trage den Winkel GPA = CPt; dann ist Pt die Polbahntangente (nach Bobillier); 


die Polbahnnormale Pn | Pt enthalt den Wendepol J. Macht man ferner CP = PW,, 
so liegt W, auf dem Wendekreis, daher ist die in W, zu PW errichtete Senkrechte ein 
zweiter Ort fiir J. 

Die Konstruktion von J und damit auch von 7' versagt in den beiden Totlagen des 
Getriebes (Kurbelwinkel y = 0 und z) und auch dann, wenn die Pole der Bewegung 
auBerhalb der Zeichenflache liegen. Fiir diese Falle und auch fiir Kontrollzwecke ist 
die rechnerische Festlegung der Lage von J und des Tragheitspoles 7’ erforderlich. 


5 Siehe FuBnote 1, S. 167. 
6 Siehe FuBnote 2, S. 168. f 
7 Mit geringfiigiger Abanderung schon angegeben von F’. Wittenbauer, Graphische Dynamik, 


Berlin: J. Springer (1923), S. 115. 
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Sei y der Drehwinkel der Kurbel aus der linken Totlage und Winkel X2AC = /}, 
ferner r/a = A (gegeben), so ist nach Abb. 1. 


sin B ey 
sine)” 
woraus 
1+ Acos ‘ 
ctgP = “Tang (3) 
Die Lage des Drehpoles P ist bestimmt durch 
aD ean ge 
ates he (y — B) 
oder nach Eliminierung von f mit Hilfe von (3) 
_ _asing — 
Chis 1+ cosy ” 2) 
worin 


w (py) = V1 + 24 cosy + 2. 


Da AC = aw, so sind die Koordinaten des Drehpoles P fiir das (&, 7) Achsenkreuz 
E=SO=AC—s=uw—s, 


aw sin | (5) 


be RUE, = ereaspe 


Die Konstruktion des Wendepoles J nach Abb. 1 fiihrt zu folgenden geometrischen 
Beziehungen 


eo ACO aw 
u 3 
OS cosB 1+ Acosq’ 
er A (A + cos ¢) 
GU =a l+Acosg@ ’ 
CP = (op ee) 
sin « 


Wird letzterer Wert aus dem Dreieck AC P folgendem Wert 


a cos B 1+ Acos@ 
cos (y — B) Me A+ cosy 


gleichgesetzt, so ergibt sich eine Gleichung fiir den Winkel « = GPA mit der Lésung 


P= 


sin p 
8B O =a paet ceca (6) 


fa A+ cos \? 
4 (py) =A (rates! ; (7) 


Durch den Winkel « ist nun W, J bestimmt; es wird W, J = PW, cotg « oder wegen 
PW, = PC mit Benutzung von (4) und (6) 


WoT 1+y cos » 
Lene % (A+ cos 9) ° 


Hiemit berechnen sich die Koordinaten des Wendepoles J fiir das (é, 7) System zu 


he ee 2 oz w {1+ (A+ 2cos @)} 
Ey SC+W,J=AC s+ W,J=a x4 (A + cos ¢) —o|,| (8) 


= Rua; um _ sing 
ny =CW, 2 eo 1 mecnega 
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Hierin bedeutet s = AS und o = =. 


Fir die linke Totlage folgt hieraus mit gy = 0: 


fy =a[So* —ol, 


gy = —a|"5™ + ol, AG 


Fiir die Sonderlage y = + wird 


Lea Mata a5 
£, =a] + i+ # —ol, 
9, Vit # 
UN, == 2a i : 
Fir gp = — > bleibt €; ungedndert, wahrend 7», ein negatives Vorzeichen erhalt. 


6) Tréagheitspol 
Nach Gl. (1) ergeben sich die Koordinaten des Tragheitspoles Tim (é, 7) System mit 
5 2 
éy = — §, ee 


42 
‘Ss 
NT, itn elep ony * 


Aus 
tp ty = §p €; + yp ny 


ergibt sich nach etwas langerer Rechnung 


et = N =Q? [zx (A+ cosy) + 47 (2+ 3A cosy + A*)] —oQ[1+y7(24+ 3 cosy)] +0, 
worin 
Ww 


Q (0) = Fat con 9)’ 


Hiemit wird schlieBlich 


f= — == [0 {1 + 7 (4h +2 con )} — ol. 
1a 
2, 2w sing i 

YE. aN A+ 08.9 * | 


Die Lage des Tragheitspoles der Schwinge ist somit abhangig vom Kurbelwinkel 9, 

vom Verhaltnis 2 der Kurbellange r zur Festpunktentfernung 21C = a und von dem die 
AS 

Lage des Schwerpunktes S der Schwinge kennzeichnenden Parameter o = ~~ = — 

Fiir die beiden Totlagen vereinfachen sich die vorstehenden Ergebnisse sehr 


erheblich. 
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Mit y = 0 wird zunachst wi (0) 


K. 
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1 
1+, ~()=4, VO)= 7, 
(0) EAE een) aa eee 
Hiemit wird 
_ ie (+44)? 
Er (0) an 1 aN (0) i GO|, (12) 
nr (0) = 0. 
Mit p=x wird w(z)=1—A, z(a)=A, 
1 
@Q (x) = Ee eg 
May 2 ee 


Hiemit ergibt sich 


; jay (ase 
ST (x) = te a N (x) A I oo ’ (13) 


IT (zr) => 0. 


In beiden Totlagen liegt der Tragheitspol auf 
der Schwerachse AC der Schwinge. Durch die 
beiden Gleichungen (11) ist die Tragheitspol- 
kurve im bewegten System in Parameterdar- 
stellung (Parameter ~) gegeben. Sie ist eine zur 
Schwingachse symmetrische Kurve, deren Form 
fiir eine schwingende Kurbelschleife mit den 
Annahmen a =27 (d.h. A= 44) und 's—@ 
(d. h. o= 1) in Abb. 3 fiir den Drehwinkel- 
bereich y = 0 bis + 110° und g = + 120° bisa 
dargestellt ist. 


Im dazwischenliegenden Bereich gy = 110° 
bis 120° nehmen die &;-Werte au®erordentlich 
rasch zu, wie die Zahlenangaben fiir die redu- 


zierten Werte &; = ee und 97= a ray Abe 
belle 2 zeigen. 

Fiir jene Kurbelstellung g = g*, in der der 
Wendepunkt W des Strahles S.J mit dem Schwer- 


punkt S zusammenfallt, also fiir SW = 0, wird 


ST = oc. Der Strahl SJ ist dann Asymptote 
der Tragheitspolkurve. 


Die Bestimmungsgleichung fiir y* lautet 
gemaf Gl. (11) 
Nae, 


Da&B mit VN = 0 auch SW — 0 wird, demnach S 
mit W zusammenfallt, zeigt die direkte Be- 


rechnung von SW. 


Der Tragheitspol der Kurbelschleife iz} 


Tabelle 1 
ree eS oe 
pe 0 30 60 90 100 110 120 130 150 180 


sT| —2,00 —1,647 —1,197 —1,485 —2,2997 —6,784 +7,464 42,744 +1,713 + 2.00 


NT 0 —0,471 —0,536 —0,369 —0,292 —0,233 0 +0,094 +0,412 0 
Tabelle 2 

a tt As ll SS fg Ee ee, See nee 

po 113 114 115 115° 5’ 115° 830” 115° 6’ 115° 7’ =115° 10’ 116 118 


éT| —16,65 —31,983 —363,8 —2844,7 —6595,7 +9233,6 +1813,2 4+507,6 +245,4 +12,367 


nT| —0,284 —0,402 —3,180 —24,05 —55,58 +77,54 115,12 44,15 +1,376 +0,017 


Nach Abb. 1 ist 


SW =SJ— PJ cose: 
da 


2 


ry? 


1 
T= ae (B® + ny) = @ [1 + 27 (A+ 2008) + 72 (4 + B+ 4 A008 9)]— 
—20Q0(1+ 2yc0sy+4+ Ay) 4+ o?, 
PW, aw 
ae x (A + cos g) VA PZ O08 Ga 23, 


e=y— a, 


PJ= 


so folgt das einfache Ergebnis 


a? 
SJ 


SW = 


N, 


wonach die Bedingung SW = 0 gleichbedeutend ist mit N = 0. 

Aus N = 0 ergibt sich g* = + 115° 5’ 41”. 

Die beiden symmetrisch zur &-Achse liegenden Asymptoten schneiden sich auf der 
&-Achse im Punkt (& = 16, 7 = 0) und sind gegen die &-Achse geneigt unter + 0°28’ 54”. 

Abb. 3 zeigt, daB im Gegensatz zum zentrischen Schubkurbelgetriebe und zum 
Kreuzschieber der Tragheitspol bei der Kurbelschleife in einem gewissen Bereich des 
Kurbelwinkels sehr heftige Lageninderungen innerhalb des bewegten Systems erfahrt. 


(Hingegangen am 18. Febiuar 1961) 
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174 A. Huber: 


(Mitteilung aus der Physikalischen Forschungsabteilung der Osterreichischen Stickstoff- 
werke A. G. in Linz-Donau). 


Ein zusammengesetzter Warmeaustauscher 


Von A. Huber, Linz-Donau 
Mit 4 Textabbildungen 


Zusammenfassung. Der in einem Ammoniak-Syntheseofen an den Kontaktraum angeschlossene 
Warmeaustauscher ist gegeniiber den sonst tiblichen Typen meistens dadurch ausgezeichnet, da} 
das kalte Frischgas einen Teil der vom heiSen konvertierten Gas aufgenommenen Warme an das 
in einem Ringspalt herangefiihrte Frischgas abgibt. Die zur Berechnung der in einem solchen Warme- 
austauscher vorhandenen Temperaturen erforderlichen Formeln werden entwickelt und an einem 
Beispiel wird der Unterschied gegen einen gewéhnlichen Gegenstrom-Warmeaustauscher gezeigt. 


I. Beschreibung des Warmeaustauschers 


Aus dem Kontaktraum K (Abb. 1) strémt das heiSe konvertierte Gas durch die 
Rohre eines Biindels nach unten und gibt dabei einen Teil seiner Warme an das zwischen 
den Rohren und Schikanen nach oben strémende Frischgas ab, das durch das Zentral- 

rohr Z abflie8t und oben in den 
K Kontaktraum eintritt. Vorher 
YON aber gibt dieses erwarmte Frisch- 
gas einen Teil der vom konver- 
tierten Gas empfangenen Warme 
an das im Ringspalt R nach 
unten str6mende Frischgas ab, 
das von unten her in den 
eigentlichen Warmeaustauscher 
eintritt. Die durch den Druck- 
mantel an die Umgebung ver- 
lorengehende Warmemenge konnte 
zwar auch theoretisch beriicksich- 
tigt werden, sie ist aber so gering, 
da sie vernachlassigt werden 
darf. 

Wir bezeichnen den Innen- 
raum aller Rohre des Biindels 
mit 7, den Raum zwischen ihnen 
mit JZ, den Ringspalt mit J// 
und die Temperaturen der in 
diesen Raumen stro6menden Gase 
mit 7 (1), ¢ (1) und @ (1), wobei 7 
vom Austritt des Gases aus dem 
Kontaktraum nach unten ge- 
rechnet werden soll. Aus der 
Abb. 2 sind die Strémungs- 


ab 


Beli ay U, =L, Ty, = . 
und Warmeiibergangsverhaltnisse 
Abb. 2. Schema der Abb. 3. Schematischer zwischen diesen Raumen sowie 
Gasstrome VerlaufderGastemperaturen die Bezeichnungen der entspre- 


chenden GréBen zu entnehmen. 
Der Verlauf der drei Gastemperaturen ist in der Abb. 3 schematisch dargestellt. 


Obwohl die Raume JI und I/J im entgegengesetzten Sinn vom Frischgas durchstrémt 
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werden, zeigen t (1) und # (J) den fiir Gleichstrom typischen Verlauf, wobei insbesondere 
die # (l)-Kurve bei / = L eine horizontale Tangente besitzen muB. 


II. Die Differentialgleichungen fiir die Gastemperaturen 
keal 


Ks seien W; ) a 


str6menden Gase, die im allgemeinen Funktionen der gesuchten Gastemperaturen 
sind. Um aber zu geschlossenen Ausdriicken fiir 7’ (J), ¢ (1) und @ (1) zu kommen, ersetzen 
wir die W; durch konstante Mittelwerte, die aber von den noch unbekannten Austritts- 
temperaturen abhangen, so da man eventuell mit geinderten W; die Rechnungen 
wiederholen mu, wie es auch sonst bei der Berechnung von Warmeaustauschern not- 
wendig ist, wenn man genauere Ergebnisse haben will. Da in JJ und JJI gleich groBe 
Gasmengen von gleicher Zusammensetzung strémen, deren Molwarmen sich nur wegen 
der héheren Temperatur in // geringfiigig unterscheiden, so wollen wir zur Verein- 
fachung der folgenden Entwicklungen noch W, = W, setzen. 


Mit U und U’ seien die Umfange eines horizontalen Schnittes mit den Ubertragungs- 


flachen von J nach JI bzw. von IJ nach JIT bezeichnet und mit k bzw. k’ = ae die 


entsprechenden Warmedurchgangszahlen, von deren Temperaturabhangigkeit wir 
jedoch wegen der ihnen anhaftenden groBen Ungenauigkeit absehen wollen. 
Die Warmebilanzen fiir eine beliebige Schichte von der Dicke di lauten sodann fiir 
die drei Raéume: 
I: W,-dT +k-U-dl-(T —) =0 
Tl: W,-dt +k-U-dl-(T —1i)—k’-U'-dl- (¢ — 3) = 0 
ITI: W,:d6+ k’- U’-dl- (8 —t) = 0, 
wozu bemerkt sei, da dT’ < 0, dt < 0, jedoch dd > 0 fiir di > 0. 
Mit den folgenden Abkiirzungen fiir die positiven Konstanten: 


(1 = 1, 2, 3) die stiindlichen Wasserwerte der in den drei Raumen 


= 4: = ~y-=° (1) 
werden die obigen Gleichungen: 
T = 4-7), (2) 
’=b:(t— T)+c¢-(¢— 9%), (3) 
0’ =c-(t — 9). (4) 


Dieses System ist nun fiir die folgenden Randbedingungen zu integrieren : 
T(0)=T7; t(L)=8(L); 9 (0) =%, 


wobei aber nur T, und %, vorgegeben sind. 


Ill. Integration der Warmebilanz-Gleichungen 


Zunachst liefert die aus (2), (3) und (4) sofort ableitbare Beziehung 
1 


i! , , 
MW = 5 (t —9') 


wegen ¢(L) = 8 (L) mit T (L) = T, das Integral 


a 


Tl) —T, =F: tO —9), (5) 


L76 A. Huber: 


was wegen 7 = ul auch physikalisch einleuchtet. Eliminiert man damit ? aus (3), 
so erhalt man zunachst 
fal te (P— TY] 
und nach Einsetzen in die differentiierte Gl. (2): 
TT =a-t—T7)=6* (14 e-(T—.17,))—a-T, 


wofiir man mit ° 
T —T,=v/(l) und b—a=2p 
schreiben kann: 
vo’ —2p-v'—be-v=0. (6) 
Das allgemeine Integral dieser Gleichung lautet mit 


q=|/r+be; a2=pty: 
v (l) = C,- e@! + Cy: e&!, 

oder, da v (L) = 0: 

v(l) = @- [ee C-)) — eal), 
Die Integrationskonstante C wird wegen 0, — oe, = 2q und (2) bestimmt durch: 

vj 29 °C =a" (ty), 
so da’: 

Tee es Spy (1 — t,)* fea ¢ = ea Ce: (7) 


Aus (4), (5) und (7) folgt nach Integration mit # (0) = 0: 


e ek eukes 
pe eS eg! = 
Qo ( ) Q1 


- (et! —1)], (8) 


Ch =p Sag (ys aa 


Damit liefern (5) und (7): 


Mit den Abkiirzungen 


bre. (lee ee = 


Cy CE SA ae Sa ena aA Be ALS | 
sa (e : ee ee | Qs 0: 


(10) 
ergeben (7) und (9) fiir ] = 0 und (8) fiir / = L die folgenden Beziehungen: 

T, = 7T,+(T, —1,)- A; 

ig = 0, (Dy ea 


Ho Oy (1 1) ae, 


Die Auflosung dieser Gleichungen nach den gesuchten Austrittstemperaturen 7’ tp une 
1, = %, liefern die folgenden Darstellungen mittels der vorgegebenen Eintrittstempera- 
turen 7’, und %, wobei 


D=1+A+B (10’) 


Kin zusammengesetzter Warmeaustauscher U7 


b 
und —- = W,/W, = w das Verhaltnis der Wasserwerte: 


Mp Tees (Lig 9,) A/D; (11) 
to = Bo + (To — Ho) w: A/D, (12) 
t, = 3, = 8 + (T, — 0) - BID. (13) 


Aus (11) und (12) folgt insbesondere 
(TZ) — Ty) * Wy = (to — %) - W>, 
wie es ja sein mu. Aus (4) sieht man sofort, daB wegen ¢, = 0, die 0-Kurve bei] = L 


eine horizontale Tangente haben muB. 


IV. Der Gegenstrom-Wiarmeaustauscher als Sonderfall 


Fir k’=0 werden c = 0, also: 0, = 2p=b—a und e,=0, so da wegen 


pitts 7 
iti? ——$= = 
¢—>0 2 


1 


zunachst aus (8) folgt: 
AL =i 5 = Uys 


wahrend (7) und (9) wegen [ue nes > die folgenden Darstellungen liefern: 


c>0%2 


POAT + oe “(r, — 4%) ee ®) 
(1) = ita (Ps —£)- (1 — 24-0), 


Um daraus die gebrauchlichen Formeln fiir Gegenstrom zu erhalten, haben wir die 
unbekannte Austrittstemperatur 7’, durch die vorgeschriebene Eintrittstemperatur 7’, 
zu ersetzen. Nach einigen Umformungen erhalt man fiir die Temperaturen in einem 
beliebigen Querschnitt im Abstand / vom Eintritt des heiBen Gases mit der Abkiirzung 


Zl) = soe tat! 
T (bl) =T, — (TT, — 4) ° 2 Y) (7’) 
t() =t, + (7) — 7,):w- [4 (L) —Z ()], (9’) 


also fiir 1 = L bzw. 1 = 0 die bekannten Formeln fiir die Austrittstemperaturen, als 
deren Verallgemeinerungen die Formeln (11) und (12) gelten kénnen. 
Dabei wird: 


hn 2 
Bee — ay 1) 
und 2p L= = -(w — 1), wobei F = UL die ganze Heizflache ist. 
1 


V. Ein Beispiel 


Fiir einen Syntheseofen, in den nach dem Warmeaustauscher noch kaltes Frisch- 
gas eingespeist wird, so daB W, > W,, seien die folgenden Werte gegeben: 


1 / 
W, = 14500; W. = 9000 f= 4,3; Sizes U' = 2.2m; 
F keal gh ; ere 
k = 410; ee (ier, eat T, = 540; 8, = 50°C. 
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Nach (1) werden damit: 


a = 0,875; b= 1.416: 


400 


200 


100 


Abb. 4. Vergleich der Temperaturen 
im zusammengesetzten und im ein- 
fachen Gegenstrom-Warmeaus- 
tauscher 


b 
c = 0,089; Md ime 1,61 
und nach (10) und (10’): 
A=1,712; B= 0,288; D = 3,000; 


also: 
A/D = 0,571; B/D = 0,096; w- A/D = 0,920. 
Damit ergeben die Formeln (11), (12), (13): 


T, = 540 — 490 - 0,571 = 260 
ty = 50+ 490-0,92 = 500 
t, = 0, = 50 + 490+ 0,096 = 97. 


Da die ganze Heizflache F = L - U = 133 m?, 
also Z (1) = 0,595, liefern dagegen (7’) und (9’) 
fiir drei Werte von ¢, die folgenden Austritts- 
temperaturen : 


KC eG 80 100 
Ts: 249 266 278 
i 519 520 621 


Die ¢, sind also fast unabhangig von ¢, und 
um rund 20°C hoher als friiher. Dieser Unter- 
schied kann aber unter Umstianden bei der Be- 
rechnung der Ammoniakausbeute zu einem 
merklich geringeren Betrag fiihren als bei Ver- 
wendung der aus der Formel (12) berechneten 
Temperatur fp. 


In Abb. 4 sind der Temperaturverlauf im 
zusammengesetzten (ausgezogene Kurve) und 
der im Gegenstrom-Warmeaustauscher fiir 
t, = 100° C (strichlierte Kurve) dargestellt. Im 


ersten Fall gibt das heiBe Gas die Warme- 


keal F : 
menge 4,05 - 108 ab, im zweiten Fall da- 


gegen werden nur 3,8- 108 ae ausgetauscht. 

_ Ich danke der Unternehmungsleitung der 
Osterreichischen Stickstoffwerke A. G. fiir die 
Erlaubnis zur Veréffentlichung dieses Auf- 
satzes. 


(Hingegangen am 6. Mérz 1961 ) 
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Buchbesprechungen 


Hauptprobleme der Bodenmechanik. Von J. Brinch-Hansen und H. Lundgren. Mit 150 Abb., XII, 
282 S. Berlin-Géttingen-Heidelberg: Springer-Verlag. 1960. Geb. DM 36.—. 


Wie schon der Titel des Buches vermuten la8t, ist sein Inhalt auf die wichtigsten Fragen der 
Bodenmechanik beschrankt. Neben einem Abri8 witber die Bodeneigenschaften und heute wiblichen 
Untersuchungsmethoden werden die Stromungsvorginge im Boden, die Setzung von Fundamenten 
und die sogenannten Bruchprobleme behandelt, wobei in diesem letzten Kapitel der Erddruck auf 
Stiitzwande, die Tragfihigkeit von Fundamenten und Pfahlen sowie die Standsicherheit von 
Boschungen zusammengefaBt sind. Hervorzuheben ist das Bestreben, die Unstimmigkeiten zwi- 
schen Theorie und Praxis bei den konventionellen Berechnungsmethoden der Bodenmechanik 
durch die Anwendung neuer Verfahren zu eliminieren. Auf gewissen Teilgebieten ist hier zweifellos 
ein wesentlicher Fortschritt erzielt worden, der nur deshalb nicht noch umfassender ist, weil an 
den herké6mmlichen Annahmen iber die Spannungs-Dehnungs-Beziehungen der Boden und ihren 
Festigkeitseigenschaften festgehalten wurde. Vollstandig abweichend von den bei uns iblichen 
Methoden ist die Art der Ermittlung der Standsicherheit von Bauwerken, indem die einzelnen, auf 
das Bauwerk wirkenden Krafte mit gewissen Faktoren, die Partialkoeffizienten genannt sind, 
multipliziert und die Bodenkennwerte ,,Kohision und Reibung“ durch ebensolche dividiert wer- 
den. Die Standsicherheitsuntersuchung mu8 dann mit diesen fiktiven Kraften und Bodenkenn- 
werten eine ,,nominelle Sicherheit‘ gré8er als eins ergeben. Wenn auch manches fir diese Methode 
angefihrt werden kann, so diirften doch die Nachteile die Vorteile iberwiegen. H. Borowicka, Wien 


The Numerical Treatment of Differential Equations. Von L. Oollatz. Translated from a supplemented 
version of the second German edition by P. G. Williams. Dritte Auflage. Mit 118 Abb., 1 Portrat, 
568 8. (Die Grundlehren der mathematischen Wissenschaften: Band 60.) Berlin-Géttingen- 
Heidelberg: Springer-Verlag. 1960. DM 93.60, Gln. DM 98.—. 


Die ersten beiden Auflagen dieses in der Zwischenzeit bereits zum Standardwerk gewordenen 
Buches erschienen 1951 bzw. 1955 in deutscher Sprache, wobei die zweite Auflage gegeniiber der 
ersten wesentliche Erweiterungen aufwies, die insbesondere die in der Zwischenzeit aktuell 
gewordene Anwendung funktionalanalytischer Methoden in der praktischen Mathematik, an der 
der Verfasser und seine Schiiler wesentlichen Anteil haben, beriicksichtigten. Die Ubersetzung 
wurde von dieser zweiten Auflage hergestellt, wobei eine Reihe kleinerer Erginzungen und Ver- 
besserungen vorgenommen und zusitzliche Anwendungsbeispiele geboten wurden. Die Uber- 
setzung wurde sorgfaltig und sachkundig vorgenommen, so daB die englische Fassung den klaren 
und iibersichtlichen Stil, den das Original aufwies, beibehalten hat. AuBerdem wurde noch die 
in der Zwischenzeit erschienene Literatur erganzend beriicksichtigt. 


Jeder, der sich auf dem Gebiete der praktischen Behandlung gewéhnlicher und partieller 
Differentialgleichungen umfassend informieren oder in einem speziellen Falle die zur Lésung 
einer Aufgabe gangbaren und zweckmaéBigen Wege kennenlernen will, wird auf das Buch von 
Collatz nicht verzichten kénnen. Besonders fiir den, der sich auf diesem Gebiet forschend betatigen 
will, ist natirlich das Buch schlechthin unentbehrlich. Fiir den Anfanger ist es gerade bei der 
Fille der vorhandenen und weit verstreut publizierten Literatur ungeheuer schwer, sich zu orien- 
tieren und zu entscheiden, was fiir ihn wesentlich ist, ja iberhaupt die Literatur zu finden. 


Das erste Kapitel des Buches stellt die wichtigsten Grundlagen, die zum Verstadndnis der 
numerischen Verfahren benétigt werden, ibersichtlich zusammen und gibt die Grundprinzipien 
an, nach denen sich die zahlreichen Verfahren zweckmaBig einordnen lassen. Auferdem werden 
die fiir die moderne praktische Mathematik wichtigen Ergebnisse der Funktionalanalysis dargelegt. 


Damit ist die Grundlage und das Einordnungsprinzip fiir die in den restlichen finf Kapiteln 
je nach Wichtigkeit ausfihrlich entwickelten und durch Beispiele belegten oder erwahnten Ver- 
fahren zur numerischen Lisung von Anfangs-, Rand- (und Eigenwert-) Problemen bei gewohn- 
lichen und partiellen Differentialgleichungen sowie von Integral- und Funktionalgleichungen 
gegeben. Dabei werden auch die verschiedenen Ansitze fir Fehlerabschditzungen besprochen. 
Es ist an dieser Stelle ausgeschlossen, auf Einzelheiten einzugehen. Es sei nur noch erwahnt, daB 
die Verwendung des Buches fiir die praktische Anwendung zur Lésung konkreter Probleme durch 
die in den einzelnen Kapiteln und im Anhang angegebenen Tabellen wesentlich erleichtert wird. 


Es hieBe Eulen nach Athen tragen, dieses Buch dem Fachmann eigens noch empfehlen zu 
wollen. Jedoch soll betont werden, da® auch derjenige, der sich nur gelegentlich mit praktischen 
Aufgaben bei Differentialgleichungen beschaftigen muB8, viel unnitze Versuche ersparen wird, 
wenn er sich zundchst in diesem Buch tiber die vorhandenen Moglichkeiten orientiert. 

E. Bukovics, Wien 
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Theorie und Praxis der Schmierung. Von D. D. Fuller. Ubersetzt und bearbeitet von W. Maier 
unter Mitarbeit von H. Paul. Mit 239 Abb., 442 S. Stuttgart: Berliner Union. 1960. DM 64.—. 


Fur den angehenden wie auch fir den in der Praixs stehenden Ingenieur ist die Schmierung 
ein wichtiges Teilgebiet seines Studiums bzw. seiner technischen Tatigkeit. Aus solchen Erwagungen 
heraus hat der Verfasser das vorliegende Buch geschrieben, das heute zu den Standardwerken auf 
diesem Gebiet gezihlt wird. 

Das Hauptziel des Verfassers war, nach seinen eigenen Worten, dem Ingenieur die Theorie der 
Schmierung in moglichst verstandlicher und brauchbarer Form naherzubringen. Zu diesem Zwecke 
leitet er die hydrodynamischen und hydrostatischen Theorien aus einfachen und leicht verstand- 
lichen Grundsiitzen ab. Durch klare Entwicklung der physikalischen Gedankengange wird trotzdem 
wissenschaftliche Strenge bei den einzelnen Ableitungen erreicht und durch die Vermittlung be- 
stimmter fester Vorstellungen wird der Ingenieur befihigt, die Schmiertheorie auf verschiedene 
Schmierprobleme anzuwenden, wobei er auch die Grenzen der Theorie kennenlernt. An Beispielen 
wird diese Beschrankung gezeigt. 

In zwilf Kapiteln, beginnend mit den Grundlagen der Viskositat und des FlieBens, werden die 
hydrostatische Schmierung, die in der letzten Zeit entwickelt wurde und zahllose, noch unerforschte 
Anwendungsmodglichkeiten in sich birgt, und die hydrodynamische Schmierung, Verluste im Lager, 
Sonderausfiihrungen, luftgeschmierte Lager, trockene Reibung und Lagerwerkstoffe eingehend, 
klar und itibersichtlich besprochen. 

Als Anhang hat der Ubersetzer, dem man fiir diese wertvolle Arbeit dankbar sein mu, ein 
Kapitel tiber Walzlagerschmierung hinzugefigt. 

Jedem Kapitel sind eine groBe Anzahl Literaturangaben beigegeben, so daf man mehr als 
400 Titel findet. Dieses Werk enthalt wirklich Theorie und Praxis der Schmierung, wie es die mo- 
derne Technik verlangt. Der Ubersetzer hat sich die Miihe genommen, ein Handbuch zu schaffen, 
das jedem Ingenieur zum Studium und Gebrauch wérmstens empfohlen werden kann. 


R. Bruniak, Wien 


Magnetverstirkerschaltungen. Von W.A.Geyger. Nach der zweiten Auflage des amerikanischen 
Originals tibersetzt und bearbeitet vom Verfasser. Mit 183 Abb., 374 8. Stuttgart: Berliner 
Union. Berlin: VEB Verlag Technik. 1959. Geb. DM 64.—. 


Das vorliegende Buch ist die vom Verfasser selbst tibersetzte und bearbeitete 2. Auflage eines 
amerikanischen Originals, dessen 1. Auflage 1954 erschienen ist und das als eines der wichtigsten 
Nachschlagewerke tiber die Magnetverstarker (Transduktoren) anzusehen ist. Das Buch behandelt 
Grundlagen, Eigenschaften und Anwendungen in einem ausgeglichenen Verhaltnis. Nach der 
Einleitung mit der geschichtlichen Entwicklung wird das ganze Gebiet in vier Hauptgruppen 
unterteilt, nadmlich: Magnetverstarker ohne Riickkopplung, Magnetverstirker mit 4uRerer Rick- 
kopplung, Magnetverstairker mit innerer Riickkopplung und Magnetverstirkerdifferenzenschal- 
tungen. Diese Hauptabschnitte zerfallen logisch je nach der Anzahl der vorhandenen Drosseln 
in solche mit 1, 2 oder 4 vormagnetisierten Drosseln. Fiir jede dieser Typen werden die Grund- 
schaltungen, die Arbeitsbedingungen und die Ausfiihrungen in besonderen Kapiteln behandelt. 
SchlieBlich werden in vier Abschlu8kapiteln Spezialausfihrungen der Magnetverstarker, die auf 
Frequenzverdopplung beruhen, weiters die technischen Eigenschaften unter besonderer Beriick- 
sichtigung des dynamischen Verhaltens der Magnetverstirker, die typischen Anwendungen und 
die Verwendung von Magnetverstarkerschaltungen zu Aufzeichnungen von dynamischen Hysteresis- 
schleifen behandelt. 


Die der Ubersetzung zugrunde liegende zweite Auflage unterscheidet sich von der ersten dadurch, 
da die in den letzten Jahren in den Vereinigten Staaten erfolgten bedeutsamen Entwicklungs- 
arbeiten zur Schaffung auferordentlich stabiler Magnetverstirker mit besonderen dynamischen 
Kigenschaften (sehr geringe Ansprechverzégerung) durchwegs beriicksichtigt sind. Dadurch, da8 
bei solchen modernen Magnetverstirkern die Ansprechzeit derart gekiirzt ist, daB — bezogen 
auf das Trigerfrequenzsystem — der Ausgangsstrom innerhalb einer Halbperiode der den Magnet- 
verstarker speisenden Wechselstromquelle seinen vollen, dem Eingangsstrom — bzw. der Eingangs- 
spannung — entsprechenden Wert annimmt, ist es méglich geworden, den Magnetverstirker mit 
besten Erfolgen in Regelungssystemen zu verwenden, die ihm noch vor wenigen Jahren grund- 
sdtzlich verschlossen waren. 


Durch die klare Darlegung der gesamten Magnetverstirkerproblematik, die gute Darstellung 
und das auf den neuesten Stand (1959) erginzte Literaturverzeichnis einschlieBlich der Patent- 
literatur wird das Buch allen Ingenieuren, Patentanwalten, Lehrern, Studenten und Neulingen 
eine wertvolle Einfihrung und Weiterbildung fiir dieses moderne Sondergebiet der Elektrotechnik 
bieten, H.von Bertele, Wien 
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Getriebelehre. Von P. Grodzinski +. Teil 1: Geometrische Grundlagen. Dritte, neubearbeitete Aut- 
lage. Von G. Lechner. Mit 131 Figuren, 164 S. (Sammlung Goschen: Bd. 1061.) Berlin: W. de 
Gruyter. 1960. DM 3.60. : 


Das vorliegende, den geometrischen Grundlagen der Getriebelehre gewidmete Bindchen, ist 
erstmals 1933 erschienen, geht aber im Grunde auf die 1912 in der gleichen Sammlung herausge- 
brachte ,,Kinematik‘‘ von H. Polster zuriick. Die jetzige Neubearbeitung unterscheidet sich von 
der Erstauflage vor allem durch den ganz neu gestalteten Abschnitt tiber die Konstruktion von 
Kurventrieben, sowie durch gelegentliche kleinere Einschaltungen. Ansonsten ist der bisherige 
Text mit allen seinen sprachlichen Nachliassigkeiten und sachlichen Mangeln leider kritiklos tiber- 
nommen worden. Dies betrifft besonders die Einfithrung in die wichtigsten Begriffsbildungen der 
Bewegungsgeometrie. Man vermifst hier saubere Definitionen, dafiir stéBt man neben zahllosen 
Ungenauigkeiten und Unklarheiten immer wieder auf Scheinbeweise, die vielfach in unqualifizier- 
barer Art mit ,,unendlich kleinen“ Differentialen operieren, und auf eine heillose Verwirrung bei 
vektoriellen GréBen, da dieselben nicht von Skalaren unterschieden werden! Als Beispiel die Be- 
schleunigungsformel auf 8. 29: p= pi+ pn=JVp2+ pn?! Manche VerstéBe sind wohl bloB als 
Versehen zu werten ; so wird etwa auf S. 15 behauptet, eine Nullstelle der Geschwindigkeit bedinge 
einen Extremwert der Beschleunigung, auf S.39 werden die Arten der Polkurven beim gleich- und 
gegenlaufigen Antiparallelogramm verwechselt u. a. m. Was soll man aber von solchen Ungereimt- 
heiten halten, wie daf die Tangenten benachbarter Bahnpunkte sich nicht blo® der Richtung, son- 
dern auch der GréBe nach unterscheiden (S. 23), oder da8 Winkelgeschwindigkeit und -beschleu- 
nigung reine Zeitmafe sind (S. 28)? 

Der anschlieBende, mehr getriebetechnische Teil kann eher gelten gelassen werden; er gewdhrt 
mit der Betrachtung von Gelenkvierecken, Nocken- und WAalzhebelgetrieben und einem Abrif8 der 
Verzahnungstheorie einen gewissen Ausblick auf den praktischen Einsatz der Kinematik. 


W. Wunderlich. Wien 


Problémes de théorie générale des oscillations et de chronométrie. Von J. Haag ¢ und R. Chaléat. 
Mit 66 Abb., 412 S. Paris: Gauthier-Villars, Eyrolles. 1960. 38.— NF. 


In diesem Buch werden schwingungsfahige mechanische Systeme mit einem und zwei Freiheits- 
graden, wie sie zur Zeitmessung bei mechanischen Uhren verwendet werden, sehr erschopfend 
dargestellt. Nach einem Kapitel, das der astronomischen Zeitmessung gewidmet ist, folgt die Theorie 
oben erwahnter Systeme, wobei auch verschiedene Dampfungsarten beriicksichtigt werden. Der 
Theorie der Storungen sind zwei weitere Kapitel gewidmet (Ankerprobleme), wo auch an Hand 
verschiedener Beispiele gezeigt wird, wie mit Hilfe von beabsichtigten periodischen ,,Storungen‘‘ 
das ebene Schwerependel naherungsweise isochron gemacht werden kann. Sehr ausftthrlich wird 
die Kompensation der Temperaturdehnung, und in einem eigenen Kapitel weiters die Theorie der 
Synchronisation behandelt. Aus diesen angefiihrten Problemen, die allerdings bei weitem nicht den 
Inhalt des Buches erschépfen, dirfte wohl ersichtlich sein, dafi dieses Werk einen gelungenen 
Beitrag zur Theorie der Uhren liefert und Ingenieuren, die sich mit derartigen Problemen aus- 
einandersetzen miissen, sehr wohl empfohlen werden kann. Auch die in grofer Zahl gebrachten 
Beispiele sowie die lakonische Kiirze des Begleittextes, der sich auf das Wesentliche beschrankt 
und andeutet, daf beim Manuskript dieses Buches ein oder mehrere Vorlesungsmanuskripte Pate 
gestanden haben, wird einen Ingenieur, der ,,in der Praxis steht‘, sehr wohl ansprechen. 


W. Peschka, Wien 


Hypersonic Flow Theory. Von W.D. Hayes und R.F. Probstein. Mit 92 Abb., 4648S. New York— 
London: Academic Press. 1959. $ 11.50. 


Wiahrend die Begriffe ,,Unterschallstromung“, ,,Uberschallstr6mung“‘ sowie das schwierige 
Gebiet des Uberganges der beiden Strémungsarten, die ,,schallnahe Strémung“ im Deutschen 
wie auch im Englischen verhaltnismafig kurz und pragnant bezeichnet sind, fehlt bis jetzt ein 
kurzer deutscher Fachausdruck fiir ,,hypersonic flow’. Man ist also gezwungen, sehr unsch6n 
von ,,hypersonischer Stroémung‘“‘, oder sehr langatmig von ,,Gasdynamik bei groBen Mach-Zahlen™ 
zu sprechen. Die moderne Flug- und Raketentechnik hat dieses Gebiet der Gasdynamik zunehmend 
zum Gegenstand intensiver Forschungsarbeit werden lassen. Die beiden Autoren haben sich zum 
Ziel gesetzt, die Ergebnisse dieser Forschung, soweit sie zuginglich sind, in ihrem Buch darzu- 
stellen, wobei sie durch eigene Arbeiten einen nicht unwesentlichen Anteil zum Inhalt des Buches 
beigesteuert haben. 

Das Buch beginnt mit Betrachtung der Voraussetzungen, die der Theorie gro8er Mach-Zahlen 
zugrunde liegen. Es folgen die Behandlung der Ahnlichkeitsgesetze, Konturprobleme, wie z. B. das 
Problem des giinstigsten Fligels bei vorgegebenem Auftrieb, die Analyse der Struktur von Ver- 
dichtungsstéBen, die Umstrémung abgestumpfter Korper, um nur einiges aus der Fille des 
Gebotenen herauszugreifen. Einen verhiltnismaig breiten Raum nehmen die Abschnitte ein, 
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in denen Zihigkeitseffekte besprochen werden. Am Ende des Buches finden wir ein Kapitel tuber 
die Strémung verdiinnter Gase und tber den Extremfall der freien Molekularstromung. Die 
Probleme des Fluges in groBen Héhen fallen in dieses Gebiet, das schon der kinetischen Gastheorie 
angehort. 

“Die Autoren vermeiden geschickt den leider sehr zur Mode gewordenen trocken-wissenschaft- 
lichen Stil. Sie haben die Miihe nicht gescheut, den Ergebnissen der Theorie, wo dies méglich war, 
Ergebnisse des Experiments an die Seite zu stellen; und die legendaren , gleichungsfreien Seiten“, 
in denen als die vielleicht schwierigste Arbeit eines Autors das Wesentliche eines Vorganges oder 
einer Methode herausgefeilt wird — in diesem Buch sind sie wieder zu finden. 

H. Bednarczyk, Wien 


Technik des Messens radioaktiver Strahlung. Von V. Kment und A. Kuhn. Mit 382 Abbildungen 
und 43 Tabellen, XIV, 602 S. Leipzig: Akademische Verlagsgesellschaft Geest & Portig K.-G. 
1960. DM 52.—. 

Das Buch enthalt eine ziemlich vollstandige Darstellung der heute gebrauchlichen MeBtechnik 
der radioaktiven Strahlung, wobei auch die physikalischen Grunderscheinungen, auf deren An- 
wendung diese MeStechnik beruht, in ausreichendem Mafe besprochen werden. Im ersten Teil 
des Buches werden auf 275 Seiten die theoretischen Grundlagen der Detektoren und ihre Anwen- 
dung behandelt und auch die Technologie der beim Bau der Detektoren erforderlichen Materialien 
beriicksichtigt. Ionisationskammer, Proportionalzihler, Geiger-Miiller-Zahlrohr, Funkenzahler, 
Szintillationszihler, Cerenkov-Zaihler, Kristallzihler, chemische und Thermodetektoren bzw. 
Dosimeter erfahren eine sorgfiltige Beschreibung, die durch zahlreiche Abbildungen, teilweise 
Konstruktionszeichnungen, veranschaulicht wird. Auch die gewohnliche und die Diffusionsnebel- 
kammer ist bericksichtigt. Im zweiten Teil werden auf 220 Seiten die elektrischen Grundlagen und 
deren Anwendung auf die in der nuklearen MefBtechnik verwendeten Schaltungen entwickelt, 
wobei auf die mathematische Beschreibung der Vorgiinge groBes Gewicht gelegt wird. Auf die 
Eroérterung der physikalischen Grundlagen der Zihlung von Strom- und Spannungsimpulsen und 
der Messung der Impulsamplituden folgt eine ausfihrliche Behandlung der klassischen Faden- 
elektrometer, der physikalischen Eigenschaften der Elektronenréhrenverstarker und ihrer Anwen- 
dung als Rohrenvoltmeter und Réhrenelektrometer, als Impulsverstarker mit zahlreichen Schalt- 
beispielen, die GroSenangaben der einzelnen Schaltelemente enthalten. Ferner enthalt dieser Teil 
des Buches die Multivibratorschaltungen mit ihren Anwendungen in den Untersetzern, Ampli- 
tudendiskriminatoren usf., Impulsdichtenmesser (rate meter), Koinzidenz- und Antikoinzidenz- 
kreise mit Réhren und Kristalldioden, Ein- und Mehrkanalimpulsspektrographen. Auf die Ver- 
‘wendung von Transistoren wird nur wenig eingegangen und auch die Nanosekundentechnik wird 
nicht eingehender behandelt, aber durch ausreichende Literaturnachweise hervorgehoben. Der 
dritte Teil des Buches, betitelt: MeBgerite und ihre Anwendungen, bringt auf ca. 50 Seiten den 
Zusammenbau der im ersten Teil beschriebenen Detektoren mit den im zweiten Teil gegebenen 
Schaltungen in den verschiedenen praktisch verwendeten Mefgeriten. Die Bilder beziehen sich 
auf Geraite des VEB Vakutronik, Dresden, und Frieseke und Hopfner, Erlangen. 

Zusammenfassend kann gesagt werden, da die Autoren sich mit diesem Buch ein groBes Ver- 
dienst erworben haben, da es in deutscher Sprache bisher kein Werk gibt, das in dieser Ausfiithrlich- 
keit unter Beriicksichtigung alterer und neuerer Entwicklungen auf diesem Gebiet ohne Bevor- 
zugung des einen Beitrages zuungunsten anderer das gesamte Gebiet der radioaktiven MeStechnik 
erfaBt. Besonders ansprechen wird den adlteren Leser, zu denen auch der Referent zaihlt, wenn er 
da und dort auf vertraute Abbildungen sto8t, die beweisen, da es auch schon vor dem Anfang 
des ,,Atomzeitalters eine ziemlich weit entwickelte radioaktive MeBtechnik gegeben hat, was in 
den meisten modernen Darstellungen nicht zum Ausdruck kommt. 

Das vorliegende Buch ist jedem fortgeschrittenen Studenten und auch jedem Wissenschaftler, 
zu dessen Aufgaben ein tieferes Eindringen in diese MeBtechnik gehért, sehr zu empfehlen. 


G. Ortner, Wien 


Elektronische Analogierechenmaschinen (Gleichstromanalogrechner). Von G. A. Korn und 
Th. M. Korn. Ubersetzt und bearbeitet von H. Katz. Mit 235 Abb., 466 S. Stuttgart: Berliner 
Union. 1960. DM 98.—. 


Das amerikanische Original war eines der ersten fiihrenden Standardwerke auf dem Gebiet der 
Analogierechentechnik. Inzwischen erschien eine stark erweiterte zweite amerikanische Auflage. 
Diese bildete die Grundlage fiir die nun vorliegende deutsche Ubersetzung, die von H. Katz in 
vorbildlicher und sachkundiger Weise durchgefithrt wurde. Die Bearbeitung erstreckt sich darauf, 
die amerikanischen Fachausdriicke moéglichst der deutschen Fachterminologie anzupassen oder 
zumindest sinngemai zu tibertragen. Das englische Mafsystem wurde dabei durchwegs auf das 
metrische Technische Mafisystem (m — s — kp) umgestellt. Andererseits wurden Formelzeichen 
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und Indizes nur dort geindert, wo im Deutschen andere Buchstaben allgemein tblich oder aus 
sprachlichen Griinden notwendig sind. 

Die Darstellung ist so gehalten, da8 das Buch sowohl der Anfinger als auch der fachkundige 
Ingenieur und Techniker mit Gewinn zu Rate ziehen wird. So zeigen die Kapitel 1 bis 3, welche 
Probleme und wie diese auf einem etwa schon vorhandenen Analogiegerat gerechnet werden kon- 
nen, wahrend Kapitel 4 bis 8 sich mit dem Entwurf von Rechenelementen bzw. von vollstandigen 
Analogierechneranlagen beschaftigen. 

Im Einzelnen bringt das erste Kapitel einen Uberblick iiber die mit Analogmethoden realisier- 
baren mathematischen Operationen. An Hand konkreter Beispiele wird im 2. und 3. Kapitel die 
praktische Durchfihrung behandelt. Hier zeigt sich, da8 die praktische Realisierung einer Schal- 
tung, obwohl sie im Prinzip recht einfach zu sein scheint, auf einige Schwierigkeiten stdBt, insofern 
ndmlich, als eine solche Rechnerschaltung stabil sein mu& und die Rechenelemente nicht tiber- 
steuert werden diirfen. Das erreicht man durch Anbringen entsprechender Ma&stabsfaktoren. 
Dieses Grundproblem der richtigen Dimensionierung wird am Beispiel linearer, gewohnlicher 
Differentialgleichungen mit konstanten bzw. verinderlichen Koeffizienten und nichtlinearer ge- 
wohnlicher Differentialgleichungen ausfihrlich behandelt. Dabei ergibt sich die Gelegenheit, 
die Wirkungsweise der einzelnen Rechenelemente genau zu besprechen. 

In den weiteren Kapiteln 4 bis 6 wird die Theorie und der Entwurf von linearen Rechenelemen- 
ten, Potentiometern, Gleichstromverstirkern, Multiplikatoren und Funktionsgeneratoren behan- 
delt. Kapitel 7 bespricht Hilfsschaltungen (Steuerschaltungen, Registriereinrichtungen) und den 
Betrieb der Rechenmaschine. Das 8. und letzte Kapitel zeigt, wie der Entwurf vollstandiger Analog- 
rechneranlagen durchzufiihren ist. 

Der Anhang enthalt Tabellen fir spezielle Analogrechnerschaltungen und ein ausfthrliches 
Literaturverzeichnis, durch das die am Schluf8 jedes Kapitels angegebenen Literaturhinweise er- 
ganzt werden. 

So ist ein Buch entstanden, das auch im deutschen Sprachraum Studenten, Ingenieuren und 
Wissenschaftlern gleicherweise groBe Dienste leisten wird. H. Scholz, Wien 


Transistoren. Von L. M. Krugman. 189 8S. Ubersetzt und bearbeitet von H. Michaelis (Elektroni- 
sche Reihe: Band 7 bis 9. Herausgegeben von A. Schure). Stuttgart: Berliner Union. 1960. 
DM 14.40. 


Das Buch wendet sich an den Techniker, Fachschiler und Studenten und ist als Einfiihrung 

in praktischer Form gedacht. Es beginnt mit wenigen Seiten tiber die Grundlagen der Halbleiter- 
physik und die Wirkungsweise der Transistoren. Nach Besprechung verschiedener Abarten der 
Transistoren wird ausfihrlich auf die Vierpolrechnung unter Zugrundelegung der Widerstands- 
matrix eingegangen. Das Kapitel ,,Transistorverstarker handelt unter anderem vom A- und B- 
Verstarker, Kaskadenschaltungen und komplementiar-symmetrischen Schaltungen. Es folgt die 
Besprechung von Oszillatorschaltungen. Den Abschlu8 bilden einige Seiten tiber das Verhalten 
des Transistors bei héheren Frequenzen. Die Frequenzabhangigkeit der Stromverstarkung wird 
angefiihrt und einige Neutralisationsschaltungen angegeben. Erst in diesem Kapitel werden die 
h-Parameter definiert und mit den z-Parametern in Verbindung gebracht. 
Im Vorwort vertritt der Verfasser die Ansicht, daB alle zum Verstandnis der Arbeitsweise und des 
Betriebes von Transistoren notwendigen Grundlagen gebracht werden. Es erscheint jedoch zweifel- 
haft, ob der Band geeignet ist, jemandem die physikalischen Vorgange im Transistor klarzumachen, 
da diese nur gestreift werden. Der Wert der ausfihrlich behandelten formalen Vierpolrechnung 
leidet erheblich an der Zugrundelegung der Widerstandsparamter, die im Falle des Transistors 
weder anschaulich noch praktisch gebrauchlich sind, da sie sich zur direkten meBtechnischen 
Bestimmung nicht eignen. Die sehr zahlreichen Formeln, die die charakteristischen GréBen (Ein- 
gangs- und Ausgangswiderstand, Leistungen usw.) durch die Widerstandsparameter ausdricken, 
dirften daher fiir die Anwendung wenig Bedeutung haben. 

Da die Schaltungsbeispiele und deren Besprechung mit einiger Ausfihrlichkeit gebracht werden, 
sind die spiteren Kapitel des Buches als gelungen zu bezeichnen. Zur Einfihrung in das Gebiet 
der Transistoren kann es jedoch auf Grund seines formalistischen Aufbaues nicht empfohlen 
werden. H. Potzl, Wien 


Impulstechnik. Von L. A. Mejerowitsch und L. G. Selitschenko. Aus dem Russischen tbersetzt von 
W. Glaser und W. Rohde. Mit 701 Abb., 680 S. Stuttgart: Berliner Union. 1960. Leinen DM 74.—. 


Der Band ist die Ubersetzung eines in dem diesbeziiglichen Prospekt als russisches Standardwerk 
bezeichneten Buches, das sich zur Aufgabe gestellt hat, die in den letzten zwei Jahrzehnten auBer- 
ordentlich angestiegene Verwendung von Strom- und Spannungsimpulsen (im Gegensatz zu den 
stationaren Formen von Gleich- und Wechselstrémen) zusammenfassend darzustellen. 

Das Originalvorwort definiert als ,,.Impulstechnik“ die Gesamtheit der Verfahren zur 
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Erzeugung, Mischung, Verstirkung und Messung elektrischer Impulse und ihrer praktischen 
Anwendung. 

Wenn man nach dieser Definition das Inhalts- und Kapitelverzeichnis der 665 Seiten Text 
ansieht, findet man, da® den Grundlagen und Anwendungen der harmonischen Analyse zusammen 
mit der Berechnung linearer Netzwerke 85 Seiten gewidmet sind, denen 533 Seiten mit Spezial- 
schaltungen folgen. Fiir ,,Messungen an Impulsen‘‘ verbleiben die restlichen 24 Seiten, wabrend 
den Anwendungen auer der Nennung von Telegraphie, hoher Geschwindigkeit, Richtfunkstrecken 
mit Impulsmodulation und Funkpeilung sowie Funkortung, Fernsehen und schlieBlich Messungen 
der Kerntechnik nebst den mit Impulsen arbeitenden Verfahren zur Steuerung und automatischen 
Regelung und elektronische Rechenmaschinen im Originalvorwort kein Abschnitt im Buch 
gewidmet ist. 

Als grundsitzlicher Mangel erscheint auch das Fehlen jedes Hinweises tiber die Erzeugung 
und Anwendung von Kondensatorentladungen fiir Réntgenblitze, Warmeerzeugung, fir akustische 
Impulse und vor allem fir Lichtblitze. 

Die Ausfiithrungen iiber die verschiedenen Schaltkreise zur Erzeugung sind eine sehr umfang- 
reiche Kompilation. Fiir den Schaltungsspezialisten wird dieselbe wahrscheinlich ein angenehmes 
Nachschlagewerk vorstellen, fir jeden allgemein an den Eigenheiten der Impulstechnik Inter- 
essierten wird aber das so einseitig orientierte Buch eine Enttduschung bedeuten. Der mathe- 
matische Teil findet sich im wesentlichen in jeder Einfiihrung in die moderne Laplace-Transfor- 
mation, der elektrische Teil weitgehend in jeder Vierpoltheorie. 

Der ungeheure Fortschritt der Elektrotechnik verlangt es, da dem nach Weiterbildung 
Beflissenen der wesentliche Sinn und der volle Umfang des fraglichen Tributs leicht zuganglich 
dargeboten wird, wenn der Titel grundsitzliche Informationen verspricht. Im vorliegenden Fall 
ist der gewahlte Titel aber vollig irrefiihrend. Der richtige Titel hatte nur ,,Berechnungs- 
methodik von impulsverarbeitenden Systemen“ hei®en dirfen. Bei dieser Beschrankung 
hatte ich mich vollig mit dem Prospekt dahingehend identifiziert, da es sich um ein systematisch 
aufgebautes Lehrbuch der Schwingungsanalyse der Impulstechnik zusammen mit den Schaltungen 
zur Erzeugung, Formierung bis zur Modulation und Demodulation handelt. 

H.von Bertele, Wien 


Principles of Feedback Control. Von OC. H. Wilts. Mit zahlreichen Textabb., X, 271 S. (Addison- 
Wesley Series in the Engineering Sciences). Reading Mass. — London: Addison-Wesley- 
Publishing Co., Inc. 1960. $ 9.75. 


Die vorliegende Darstellung ist fiir fortgeschrittene oder graduierte Studenten und fir bereits 
in der Praxis stehende Ingenieure geschrieben, die sich hinreichende Kenntnisse tiber die Theorie 
selbsttatiger Regelungen verschaffen wollen. Voraussetzung fir das Verstandnis sind Kenntnisse 
uber gewohnliche Differentialgleichungen und tber Funktionentheorie, fiir einige Abschnitte auch 
die Grundbegriffe der Laplace-Transformation. 

Es ist dem Verfasser gelungen, auf verhaltnismafig kleinem Raum einen sehr sch6nen und weit- 
gehenden Einblick in das heute schon allgemein anerkannte wichtige Gebiet zu geben. Gerade hier 
hat sich die stark vereinfachende und eine einheitliche Betrachtung scheinbar verschiedenartiger 
technischer Einrichtungen erméglichende Kraft der Verwendung mathematischer Ideen und Be- 
griffsbildungen gezeigt und es dirfte heute wohl keinen Regelungstechniker geben, der nicht die 
Wichtigkeit der Theorie fur sein Fachgebiet zumindestens einsehen wiirde. Wenn er nun vor die Not- 
wendigkeit gestellt wird, sich diese Kenntnisse anzueignen, wird ihm eine pragnante, alles Uber- 
flissige vermeidende und doch leicht verstandliche Darstellung wie die vorliegende sehr willkom- 
men sein, zumal diese auch den mathematischen Apparat so gering wie méglich halt. 

Nach einer kurzen Einleitung werden zunichst die wichtigsten Eigenschaften der Lésungen 
linearer Differentialgleichungen mitgeteilt, sodann wird auf das Ubertragungsverhalten linearer 
Regelsysteme eingegangen. Drei weitere Abschnitte sind der Stabilitiitsuntersuchung mit Hilfe 
der algebraischen Stabilitatskriterien, des Wurzelortverfahrens und des Nyquist-Kriteriums 
gewidmet. Dann wird die Stabilisierung von Regelsystemen besprochen. Weiters wird auf vermaschte 
Regelkreise eingegangen. Nach der Besprechung des Ubertragungsverhaltens geschlossener Regel- 
kreise wird noch Einiges tiber die statistischen Methoden zur Untersuchung von Regelsystemen mit- 
geteilt. Ein Schlu8kapitel zeigt die Behandlung nichtlinearer Systeme mit Hilfe der Beschreibungs- 
funktion. Im Anhang werden einige wichtige Erginzungen gebracht. 

Besonders wertvoll sind auch die zahlreichen instruktiven Abbildungen und die den einzelnen 
Abschnitten beigefiigten Aufgaben, die teils theoretische Erginzungen, teils reinen Ubungsstoff 
enthalten. Das Buch kann fiir den ins Auge gefa8ten Leserkreis sehr empfohlen werden. 


H. Bukovies, Wien 


Higentiimer, Herausgeber und Verleger: Springer-Verlag in Wien I, Mélkerbastei 5. — Fir d 
wortlich: Prof. Dr. H. Parkus, Wien IV, Technische Hochschule, Karlsplatz 13. — Druck: sbesameetny Wace oy 
GumpendorferstraBe 40 —44. ‘ j 
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besserte Auflage. Mit 169 Textabbildungen. X, 440 Seiten. Gr.-8°. 1960 
S$ 290.—, DM 46.—, sfr. 49.50, $ 11.50 
Ganzleinen 8 312.—, DM 49.50, sfr. 53.20, $ 12.40 
Zweiter Band: Integration und Differentiation der Funktionen von mehreren Veriinderlichen. 
Lineare Algebra. Tensorfelder. Differentialgeometrie. Zweite, neu bearbeitete Auflage. Mit 
136 Textabbildungen. VIII, 401 Seiten. Gr.-8°. 1958. 
S 290.—, DM 46.—, sfr. 49.50, $ 11.50 
Ganzleinen § 312.—, DM 49.50, sfr. 53.20, $ 12.40 
Dritter Band: Gewdéhnliche und partielle Differentialgleichungen. Variationsrechnung. 
Funktionen einer komplexen Veriinderlichen. Zweite, mit Berichtigungen versehene Auflage. 
Mit 107 Textabbildungen. XII, 512 Seiten. Gr.-8°. 1960. 
S$ 290.—, DM 46.—, sfr. 49.50, $ 11.50 
Ganzleinen § 312.—, DM 49.50, sfr. 53.20, $ 12.40 


Soeben erschien: 


Vierter Band: Integralgleichungen. Laplacetransformation. Randwertprobleme bei gewéhn- 
lichen Differentialgleichungen. Grundziige und Randwertaufgaben der Potentialtheorie. Mit 
49 Textabbildungen. VI, 335 Seiten. Gr.-8°. 1961. 
S 290.—, DM 46.—, sfr. 49.50, $ 11.50 
Ganzleinen § 312.—, DM 49.50, sfr. 53.20, $ 12.40 
Mit diesem Band liegt das Werk vollstindig vor. 


Grundziige der Tensorrechnung in analytischer Darstellung. Von Dr. phil. 


Adalbert Duschek, weiland o. Professor der Mathematik an der Technischen Hochschule 
Wien und Dr. techn. August Hochrainer, a. o. Professor an der Technischen Hochschule Wien, 
Direktor des Hochspannungsinstitutes und der Hochspannungsschaltgeratefabrik der AEG, 
Kassel. Vollstandig in drei Teilen. 


Teil I: Tensoraigebra. Vierte, erganzte Auflage. Mit 34 Textabbildungen. VIII, 171 Seiten. 
8°. 1960. 

Steif geheftet S 144.—, DM 24.—, sfr. 24.60, $ 5.70 
Teil II: Tensoranalysis. Zweite, erginzte Auflage. Mit 61 Textabbildungen. VII, 334 Seiten. 
8°. 1961. 

Steif geheftet S 246.—, DM 39.—, sfr. 41.90, $ 9.75 
Teil III: Anwendungen in Physik und Technik. Mit 25 Textabbildungen. VI, 250 Seiten. 8°. 1955. 

Steif geheftet S 144.—, DM 24.—, sfr. 24.50, $ 5.70 


Integ raltafei. Herausgegeben von Wolfgang Grébner, o. Professor an der Universitat Innsbruck, 
und Nikolaus Hofreiter, o. Professor an der Universitat Wien. 


Teil I: Unbestimmte Integrale. Dritte, verbesserte Auflage. Mit 2 Textabbildungen. VIII, 


SS ol i le Steif geheftet S 136.—, DM 22.70, sfr. 23.20, $ 5.40 


il I : i iten. 4°. 1961. 
Teil IL: Bestimmte Integrale. Dritte, verbesserte Auflage. VI, 204 Seiten. 
Steif geheftet S 170.—, DM 27.—, sir. 29.—, $ 6.75 


————E—EEEE ne 
Zu beziehen durch Ihre Buchhandlung 


SPRINGER-VERLAG IN WIEN 


Soeben erschien: 


Programmieren von Ziffernrechenanlagen. Von Dr. phil. Walter Knédel, Professor 
der Mathematik an der Technischen Hochschule Wien. Mit 18 Textabbildungen. VIII, 202 
Seiten. Gr.-8°. 1961. r 

Steif geheftet S 126.—, DM 20.—, sfr. 21.50, $ 5.— 
Ganzleinen S 145.—, DM 23.—, sfr. 24.70, $ 5.75 


Die Beniitzer von Ziffernrechenautomaten stehen vor der Tatsache, da diese Automaten weder 
die Formelsprache des Mathematikers, Technikers oder Naturwissenschaftlers noch die Fach- 
sprache des Kaufmanns unmittelbar verarbeiten kénnen. Das erfordert fur jedes Problem die 
Erstellung eines maschinengerechten Programms. 

Das Buch befa8t sich mit den hierbei auftretenden Aufgaben, wobei sich der Verfasser — gestiitzt 


auf seine nunmehr siebenjahrige Vorlesungserfahrung und praktische Arbeit auf diesem Gebiet 


an der Technischen Hochschule Wien — bemiht hat, Wesentliches klar herauszustellen, ohne 
sich in zu viele Einzelheiten zu verlieren. Der Leser soll in die Lage versetzt werden, einfache 
Maschinenprogramme selbst herzustellen oder von anderen geschriebene Programme zu verstehen. 


Die im Text verwendeten Beispiele sind dabei so gehalten, daB besondere Kenntnisse aus der 


mathematischen oder kaufmadnnischen Verfahrenstechnik nicht bendtigt werden. 


Electromagnetic Separation of Radioactive Isotopes. Proceedings of the Inter- 
national Symposium held in Vienna, May 23—25, 1960. Edited by M. J. Higatsberger and 
F. P. Viehbéck, Osterreichische Studiengesellschaft fir Atomenergie Ges.m.b.H. Vienna. 
With 168 figures. VIII, 318 pages. Gr.-8°. 1961. 5; ; 
S 340.—, DM 54.—, sfr. 58.10, $ 13.50 
Full cloth S 366.—, DM 58.—, sfr. 62.40, $ 14.50 


Im Mai 1960 veranstaltete die ,, Osterreichische Studiengesellschaft fir Atomenergie Ges. m. b. H.“ in 


Wien ein Symposium iiber ,,Elektromagnetische Trennung radioaktiver Isotope‘‘, an dem Fach- — ; 


leute aus zwélf Landern einschlieBlich der USA, GroSbritanniens, Frankreichs und der skandinavi- 
schen Lander teilnahmen, die das Thema in mehr als 30 speziellen Vortragen und Diskussionen 
behandelten. Der angekindigte Bericht bringt die Ergebnisse dieses Symposiums, die fir Wissen- 
schaftler vieler Forschungsgebiete von Bedeutung sind. — rid, 


Fortschritte im Bau und Betrieb von ElektrolichtbogenGfen. ,,Berg- und 
Hiuttenmadnnische Monatshefte‘ der Montanistischen Hochschule in Leoben. Schriftleiter: 
Prof. Dr.-Ing. K. Zeppelzauer, Leoben. (106. Jahrgang, Heft 9/1961.) . 


S 76.—, DM 12.—, sfr. 12.90, $ os 
Inhalt: Trenkler, H., Fortschritte im Bau und Betrieb von Elektrolichtbogenofen. — Original- 


arbeiten: Forst, H.-J., Lichtbogenschmelzofen groSer Leistung. — Liinig, H., Lichtbogenstahl- 
schmelzéfen als elektrische Verbraucher und Méglichkeiten zur Steigerung der Wirtschaftlichkeit. — 
Erni, E., Konstruktive Besonderheiten der 40-t-Lichtbogenéfen des Stahlwerkes Gerlafingen der 


Gesellschaft der Ludw. von Roll’schen Eisenwerke A. G. — Pungartnik, K., Erfahrungen mit den 


120-t-Lichtbogenéfen des Hiittenwerkes Rheinhausen. — Schéberl, A., und E. Horst, Der Einflu8 
einer induktiven Badbewegung auf den Schmelzablauf im 30-t-Lichtbogenofen. — Honig, F., 


Lichtbogenéfen — Zustellungsfragen, betrachtet von der Seite des Steinerzeugers. — Leitner, K.,. 


Zustellungsfragen von Elektrolichtbogenéfen. — Schéberl, A., Zustellungsfragen an Lichtbogen- 
dfen. — Schmidt, F., Zustellungsfragen an Lichtbogenéfen. — Schlager, W., Zustellungsfragen 


von Lichtbogenéfen. — Dobrowsky, F., Erfahrungen beim Betrieb eines sauer zugestellten 7,5-t- 


Elektrolichtbogenofens. — Petz, M., Sauerstoffanwendung im Lichtbogenofen. — Plessing, R., 
Sauerstofffrischen unlegierter und legierter Stahle im basischen Lichtbogenofen. — Diskussionen. — 


Schlu8wort und Zusammenfassung. — Umschau. — Buchbesprechungen. — Mitteilungen der 


»,Eisenhiitte Osterreich. — Hochschul- und Personalnachrichten. — Patente. 
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